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PREFACE 


The idea of expanding an arbitrary function in terms of the solutions 
of a second-order differential equation goes back to the time of Sturm 
and Liouville, more than a hundred years ago. The first satisfactory 
proofs were constructed by various authors early in the twentieth 
century. Later, a general theory of the ‘singular’ cases was given 
by Weyl, who- based it on the theory of integral equations. An 
alternative method, proceeding via the general theory of linear 
operators in Hilbert space, is to be found in the treatise by Stone 
on this subject. 

Here I have adopted still another method. Proofs of these expan- 
sions by means of contour integration and the calculus of residues 
were given by Cauchy, and this method has been used by several 
authors in the ordinary Sturm-Liouville case. It is applied here to 
the general singular case. It is thus possible to avoid both the theory 
of integral equations and the general theory of linear operators, 
though of course we are sometimes doing no more than adapt the 
latter theory to the particular case considered. 

The ordinary Sturm-Liouville expansion is now well known. I 
therefore dismiss it as rapidly as possible, and concentrate on the 
‘singular’ cases, a class which seems to include all the most interesting 
examples. In order to present a clear-cut theory in a reasonable 
space, I have had to reject firmly all generalizations. Many of the 
arguments used extend quite easily to other cases, such as that of 
two simultaneous first-order equations. 

It seems that physicists are interested in some aspects of these 
questions. If any physicist finds here anything that he wishes to 
know, I shall indeed be delighted; but it is to mathematicians that 
the book is addressed. I believe in the future of ‘mathematics for 
physicists’, but it seems desirable that a writer on this subject should 
understand physics as well as mathematics. 


NEW COLLEGE, OXFORD, 

1946. 


E. C. T. 
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THE STURM-LIOUVILLE EXPANSION 


1.1. Introduction. Let L denote a linear operator operating on 
a function y ^ y{x). Consider the equation 

Ly - Xy, (1.1.1) 


where A is a number. A function which satisfies this equation and 
also certain boundary conditions (e.g, which vanishes x = a and 
X b) is called an eigenfunction. The corresponding value of A is 
called an eigenvalue. Thus if is an eigenfunction corresponding 
to an eigenvalue A,^, 

LU^)^XMx}. ( 1 . 1 . 2 ) 


The object of this book is to study the operator 


L 


q{x)- 




(1.1.3) 


where q{x) is a given function of x defined over some given interval 
(a, fe). In this case y satisfies the second-order dilfercntial equation 

^_j_{A-g'(a;)}y =0, (1.1.4) 


and satisfies 

= 0 . ( 1 . 1 . 5 ) 

If we take this and the corresponding equation with m instead of n, 
multiply by ^Jix), ^J(x) respectively, and subtract, we obtain 

Hence 

b 

j = [^Jx)lli'^(x)->f)^(x)lp'^{x)]l 

a 

^ 0 


if and if^nix) both vanish at x ~ a and x ~ b {or satisfy a more 
general condition of the same kind). If A^ ^ A^^ it follows that 

b 

J dx = 0 . ( 1 . 1 . 6 ) 
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By multiplying if necessary by a constant we can arrange that 

b 

j^l{x)dx = l. ( 1 . 1 . 7 ) 

a 

The functions then form a normal orthogonal set. 

Our main problem is to determine under what conditions an 
arbitrary function f(x) can be expanded in terms of such functions, 
in the manner of an ordinary Fourier series. If this is possible, and 
the expansion is 

oo 

f{x)= 2 (1-1 -8) 

then on multiplying by and integrating over (a, 6), we obtain 

formally 

ft 

Cm = / dx. ( 1 . 1 . 9 ) 

a 

In some cases the eigenvalues are not discrete points, but form 
a continuous range, say, for example, over (0, oo). The expansion 
then takes the form 

00 

fix) = J ciX)U^) dX. (1.1.10) 

0 

All this has its simplest illustration in the case of ordinary Fourier 
series. Suppose, for example, that q(x) = 0, and that the interval 
considered is (0, tt). The solution of ( 1 . 1 . 4 ) which vanishes at x == 0 
is then y — sin(a:VA). This vanishes at a: — 77 if and only if A ™ n^, 
where n is an integer. These then are the eigenvalues, and the corre- 
sponding eigenfunctions are the functions sinno;. That an arbitrary 
function can be expanded in terms of these functions is the familiar 
theorem on Fourier’s sine series. 

1 . 2 . An argument which has sometimes been used to suggest the 
validity of the above expansions runs as follows. Consider the partial 
differential equation 

L/-4 (1.2.1) 

where / = /(rK,<). If f(x,t) is given for one value of t, it is fixed by 
this equation for a slightly greater value of t. Thus we should expect 
to have one solution, and only one, for any given initial value of 
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f{xj), i.e. for f{xj) equal to an arbitrary /(a:) when t ~ 0. Suppose 
now that the solution can bo expressed as a Fourier integral, 


00 



MO = J F{x,X)e-^dX, 

—00 

(1.2.2) 

where 

00 

F{x, X) = J f{x, /)e»^ dt. 

(1.2.3) 

Substituting in 

(1.2.1), we obtain 



oo oo 

J er^LF{x,\) dX === J Xe-i^‘F{x,X) dX. 

— 00 —00 


Since this holds for all values of we can equate the coefficients of 
which gives 

LF(x,X) — XF{x,X). 

Thus F{x,X) is an eigenfunction of the operator L belonging to the 
eigenvalue A. Now putting ^ = 0 in (1.2.2), we obtain 

00 

" jh) J 

— 00 

which gives an expression for the arbitrary function f{x) in terms of 
eigenfunctions. If/(^, were expressible in a Fourier series instead 
of an integral, we should obtain similarly a series expansion. 

The difficulty of justifying directly an argument such as the above 
is obvious. 

1.3. The argument assumes, for one thing, that f{x, t) is small as 
t -> oo, since this is required for the Fourier integral formula (1.2.2), 
(1.2.3) to hold. However, a more general form of the formula is as 
follows. Let 

oo 

J f(x,t)e}^dl (1(A) > 0), (1.3.1) 

0 

0 

F_{x,X) = ^^^f{x,t)eX^dt (I(A)<0). (1.3.2) 
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The inverse formula is then 


= +V(fc) / ■ 

tc-00 ic'-O) ( 1 . 3 . 3 ) 

where c > 0, c' < 0. Using this and (1.2.1), we obtain formally 

00 

LF.,(x,\) = J Lf{x,ty^dt 


V(277)J dt 


00 


if f{x, t) reduces to f{x) when i = 0. Similarly 


(1.3.4) 


LF_{x,X)= ^f(x)+XF_{x,X). (1.3.5) 

The method to be employed is therefore as follows. We construct 
the solutions and of (1.3.4) and (1.3.5) which satisfy given 
boundary conditions. Then (1.3.3), with / = 0, gives 

iV+oo it*' 4- CO 

ic — oo ic' — CO 


In the simplest cases F_(x, A) is found to be minus the analytic con- 
tinuation of J^^_(:r,A) across the real axis. Writing 



d)(a;,A) = -i^{2rr)F4x,X) 

(1.3.7) 

(1.3.4) becomes 

(L^Xmx.X) = -fix). 

(1.3.8) 

and (1.3.6) becomes 

ic' "h CO ic — oo « 


fix) 

=2^ J + / 

(1.3.9) 


ic' — CO ic + oo 


The expansion is then obtained from the calculus of residues, the 
terms of the series being the residues at the poles of 0(a:, A). 



1.3, 1.4 


THE STURM-LIOUVILLE EXPANSION 


5 


In any case, since (1.3.4) and (1.3.5) differ only in having the sign 
of i changed, the two terms in (1.3.6) are conjugate. Hence we 
also have ,v^ oo 

/(:r) = R y|) J F,.(x,\)dX 

ic — oo 

( ic+Qo V 

- J <D(a;,A)(/A . (1.3.10) 

iC — 00 

The expansion formula is obtained from this by making c -> 0. 

The above argument indicates in particular that A must be treated 
as a complex variable; but the analysis is of course still })urely formal. 


1 . 4 . In the particular case in which the operator L is given by 
(1.1.3), (1.3.8) is the second-order differential equation 

^f(x), (1.4.1) 

In this case the function (I)(a;,A) can be expressed in terms of the 
solutions of (1.1.4). Let Wy(^,i/i) or IF(^,0) denote the Wronskian 
TF(^, i/r) — (f){x)t/j'(x)~^'{x)ip{x) of two functions ^ and 0. Now let 
^(a:,A) and i/j{x,X) be two solutions of (1.1.4) such that lF((;t,j/f) — 1. 
Then a solution of (1.4.1) is 


X b 

0(a;,A) = #r,A) J <t>(y,X)S(y) dy A) J ^y,X)f{.y) dy, (1.4.2) 

a X 


as is easily verified by differentiating twice. 

Another starting-point for the theory, which (in the case of 
(1.1.3)) avoids an appeal to Fourier integrals, is as follows. Suppose 
the theory already established, and consider the properties of the 
function 

(^{x) 0(.t,A) = 

71 = 0 

This gives 


00 

2 


c„,'/'„(a^) 
A-A„ ■ 


(1.4.3) 


00 



-c„A„</-„(x) 


2 


a-aJ 
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This is (1.4.1) again, so that we are again led to consider solutions 
of this differential equation. If we can solve it, (1.4.3) then indicates 
that the terms of the expi^ion of f{x) will be the residues at the 
poles of ^)(a:;, A). 

Our general method consists of defining <I)(a:, A) by (1.4.2); integra- 
tion round a large contour in the complex A-plane then gives the 
value f{x), and the singularities of (!>{x, A) on the real axis give a series 
or an integral expansion as the case may be. 

1.5. The Sturm-Liouville expansion. We shall suppose 
throughout that q{x) is a real function of x, continuous at all points 
interior to the interval (a, fe) considered. In the classical Sturm- 
Liouville case (a, b) is a finite interval, and q{x) also tends to finite 
limits as a:: -> a and x~^b. 

The general theorem on the existence of solutions of (1.1.4) is as 
follows. 

Theorem 1.5. If q{x) satisfies the above conditions, and oc is given, 
the equation (1.1.4) has a solution ^(x) (a ^x ^b) such that 
(j)(a) = sina, ^'(^) “ —cos a. 

For each x, ^{x) is an integral function of A. 

Let 2/o(^) ~ sina:— (:r— «)cosa:, 

and for n == 1, 2 ,..., 

O' 

y,A^) =" J {?(<)— 

a 

Let |g(x)| < M, |yo(^)l ^ a ^x ^b, and let |A| < W. Then 

X 

2/o(^)l < I {M-\-N)K{x~t)dt ^ l{M-\-N)K{x—af. 

a 

For n ^ I, 

X 

yni^)~yn-li^) = / dt, 

a 

X 

\yn{x)-y ^ (-^+-^)(6-«) J \yn-lii)~ynS)\dt. 

a 

\yi{x)~yi{x)\ < J {t~afdt 

a 

_ K(M+Nf(b-a){x-af 
~ 3! 


Hence 
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= 2/o(*)+ 2 

n=l 

converges, uniformly with respect to A if |A| ^ N, and with respect 
to X over a < x < b. Since for > 2 


a 

the first and second differentiated series also converge uniformly with 
respect to x. Hence 

^"(x) == f {yn{^)-yl-i{x)} 

= {?(a:)-A}[2/o(x)+ {y n-\{^)—y n-M)]\ 

^ {9(a;)-A}<^(a;), 

SO that ( l ){ x ) satisfies (1.1.4). It also clearly satisfies the boundary 
conditions. 


1.6. Now let (j}{x,X), y(:r,A) be the solutions of (1.1.4) such that 
<l>{a,X) — sin a, ^'(u, A) —cos a, \ 

X(6,A) = siii/S, — cos^. 1 ^ ^ 

Then 

== fe(a:)-A)^(x);^(x)-{«2'(a;)-A}x(a:)9i(«) 

= 0 . 

Hence W (^, y) is independent of x, and so is a function of A only 
say a>(A). It is clear from the above theorem that it is an integra 
function of A. 

Let 

X b 

A) = J <f>{yA)fiy) dy J xiy>^)Ay) dy. 

a X 

( 1 . 6 . 2 ) 
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It is at once verified by differentiation that <I)(a:,A) satisfies (1.4.1), 
and also the boundary conditions 

<[)(«, A)cos a:+0'(«, A)sin a ~ 0, \ 

<5(6, A)cos^+®'(ft, A)sin^ == 0, / ^ ^ 

for all values of A. 

Suppose that the only zeros of a>(A) are simple zeros Aq, A^,... on 
the real axis. Then the Wronskian of (f){x,X,^) and x(^, A„) is zero, 
so that xi^yKi) ^ constant multiple of (^(o;, A^J, say 

y(^*,AJ -- A:,,^(.r,AJ. (1.6.4) 

It follows from the boundary conditions that is neither 0 nor oo. 
Hence <5 (a:, A) has the residue 

b 

^1/ 

^ (\) J 


at A ™ A„. The above formalities therefore indicate that there should 
kf' an expansion of the form 


/(^) = 


I 

tj, =fi 


b 

-r7hM^>K) f Hy>.K)f{y)dy- 
(K) J 

a 


(1.6.5) 


This is the Sturm-Liouville expansion. 

If we start with any two indei^endent solutions <^o(.r, A) and 
of (1.1.4), and write a>o(A) = 

A) == «+Xo( «)si “ Q'+<^ fl( ff-)si n «} 

c<Jo(A) 

and y(:r. A) is obtained by replacing a, a: by 6, Then 


<f>a( b, A Jcos^ +^o(6 , A J sin^ ^ Xo(^. A„)cos^+Xo(^. A J siii^ ^ ^ 
9So(a,A„)cffsa-f^;(a,AJsina xo(«.'^JsinaH-xo(a, A„)sma; 


and the analysis proceeds as before. 

The case of ordinary Fourier series is obtained by taking q{x) = 0. 
The equation (1.1.4) is then 


dx^ 


+ A«/= 0. 


Solutions are <f>Q{x,X) = cosajVA, Xo(^>^) — sina:VA, and a>o(A) = VA. 
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Take first the case a: — 0, == 0. Then 

sin((6— a)\/A) 

“ — a— • 

The zeros of w(A) are A,j = {mTl(b—a)Y {n ^ 1, 2 ,...). Wc have 
(K) = ^cos[(h-a)VA„} = — 2A ’ 

and h COH{bn7r/(b—a)} 

<^„(«,A„) cos{an7r/(b~a)}-^ 

Hence (1.6.5) gives Fourier’s sine series 






Similarly, by taking a ~ ^ we obtain the cosine series 


0 

/(^) = J /(y) <iy- 


3|f 


cosmTrf— -]/(?/) (Zy. 
b~ar ^ 


1.7. We have now to make the above analysis rigorous. We begin 
by proving some lemmas. 

Lemma 1.7 (i). Let </>(a%A) be the solution of (1.1.4) defined 

in § 1.5, and let A — Then 
(f){x) ~ cos{5(a;--a)}sina:-— 

X 

_sinsix (^}(.osa+- 

S S J 

a 

The last term is equal to 


X 

Z J mi{s{x~y)}{s-f{ij)+f{y)} dy 

a 

by the differential equation. On integrating by parts twice we obtain 

X 

J sin{,s(a;— y)}^"(y) dy — ,s^(:r)+sin{5(a;— a)}cosa:— 

a 

X 

---5COs{,s(a;— a)}sina— J sin{5(a;— y)}^(y) dy. 

This proves (1.7.1). 

4957 n 
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Lemma 1.7 (ii). For s — |5| > sin a ^ 0, 

^{x) (1.7.2) 

= cos{5(a;--a)}sina:4-0(|5|“%l^‘(‘*^“®^), (1.7.3) 

while if sinoL ~ 0, (l>{x) ™ 0(e'^i(^“®y |5|), (1.7.4) 

,, , sin5(x— a) , 

^{x) = >cos<x+0\-^j. (1.7.5) 

Each result holds uniformly for a ^ x ^ b. 

Putting (l){x) — (1.7.1) gives 


F{x) — lco8{s{x-~a)]sina~^^^^"'^^^^ — cosa:je”i^'^‘^“®^+ 

X 

+ i J 8\n{s{x—y)}e~^^^^^-y\{y)F{y) dy. 

a 

Let /X = max \F{x)\. Then it follows that 

a<,x<,h 

^ < |sina| + i^j^ + jt J \q(y)\,j,dy, 

a 

i.e. |sina| + |cosa|/|5| 

y. ^ , 

1— jlj J \<i(y)\ dij 

a 

provided that the denominator is positive. This is true if l^] is large 
enough, and (1.7.2) and (1.7,4) follow. Also (1.7.3) follows on sub- 
stituting (1.7.2) in the integral on the right of (1.7.1); and similarly 
for (1.7.5). 

Of course a general asymptotic expansion of ^(o;) as a function of 
s can be obtained by repeating the process. 

Lemma 1.7 (iii). //sin a 0, 

^'(x) — 5sin{^(a:--a)}sina:-f (1.7.6) 

and if sin a — 0, 

<j>'{x) — — cos{<s(a:— a)}cosa:+0(e'^l(®-®Vl^|). (1-7.7) 

This follows at once on differentiating (1.7.1) and using Lemma 
1.7 (ii). 

Similarly we obtain 

xix) ~ cos{5(ft— a:)}sin^, x(^) s &in{s(b—x)}8in ^ (1.7.8) 
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(or similar formulae if sirijS = 0); and 

co(A) ~ s[coH{s(x—a)}sin{s{h~x)]~\-sm{s(x—a)}cos{s{b—x)}]sin(xm\^ 
™ ^sin{5(6~-a)}sinasin^, (1.7.9) 

or similar formulae in the other cases. It follows in particular that 
CO (A) is not identically zero. Actually co is an integral function of s 
of order 1, and so an integral function of A of order 

1.8. Orthogonal property of the expansion. If L ~ q{x)~-d^ldx^, 
and F and G are any functions of x with continuous second 
derivatives, 

b b 

J F.LGdx - J G.LFdx 

a a 

b 

= - J (FG"-GF") dx = W,,(F, G)-Wf,(F, G). (1.8.1) 

a 

If F and G satisfy the same boundary conditions at a and at 6, so 
that F'jF ~ G'jG at these points, it follows that 

b b 

J F.LGdx^ J G.LFdx. (1.8.2) 

a a 

Now let F ~ (f>{x,X)y G “ ^(:r,A') be solutions of (1.1.4) satisfying 
(1.6.1). Then 11^(J^, (r) = 0. If A is a zero of <^(A), the Wronskian 
of ^{x,X) and vanishes, so that 

' x(*> 

where i is a constant, which is neither 0 nor oo, by the boundary 
conditions. Hence by (1.6.1) 

... .. sin)8 .. coS)8 

<)i(^,A) == jfy ^ ^ 

Similarly, if A' is a zero of co(A), </>(.r,A') satisfies the same conditions 
at X — b. Hence 

^ b b 

J (f){x,X)L(f>{x,X') dx = J (j){x,X)L^[x,X) dx. 

a a 

Since L(f>(x,X) — A^(a;,A), L^{x,X') ~ X'<f>(x,X') it follows that 

b 

(A-A') J 4>(x,X)<l>{x,X') dx = 0. 


(1.8.3) 
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b 

Hence, if A A', J (l){x,X)(f){x,y) dx — 0, (1.8.4) 

a 

the orthogonal property. 

If A “ u~\-w were a complex zero of cl)(A), then so w'. 

A' ™ A u—iv 

be. Since ^(x,X) is the conjugate of ^(x,X), (1.8.4) gives 

b 

J |<^(:r, A)|“ dx ™ 0, 

a 

which is impossible if cj) is not identically zero. Hence all the zeros 
of aj(A) are real. 

Again, taking F ^(a;. A), G = </>(a%A') in ^_8.1), we obtain 

b 

(A-A') [ cl>(x,X)cl>{x,X') dx = -<i^{h,X)<i>'(b,X')+4>'{b,X)<t,{b,X'). 

a 

Also aj(A) = —^(6, A)cos/3— ^'(6, A)sin/3. 

Hence, if sin^ 7 ^ 0, 

b 

(A-A')sin 4,{x,X)4,{x,X') dx w{X')<f>{b,X)-w(Xmb,X'). 

a 

If Xq is a double zero of a)(A), then (^(A^ii^^) 0(v-) as v 0. Hence, 

taking A AQ+^^^ A' — Aq—/^;, the right-hand side is 0{v^). But the 
left-hand side ^ 

2ivsinp j \(l){x,X)\^ dx, 

a 

SO that we obtain a contradiction. Hence all the zeros of co(A) are 
simple. If sin ^ — 0, 

b 

(A-A') J 4>{x,Xmx,X') dx = rt{w(A)rA'(/>,A')-a.(A')^'(6,A)}, 

a 

and the result again follows. 

It also follows from (1.7.9) that oj{X) ^ 0 for s = it {t > Iq), i.e. for 
A negative and sufficiently large. 

1.9. Theorem 1.9. Let f(y) he integrable over {a, b). Then if 
a <, X a b the Sturm-Liouville expansion (1.6.5) behaves as regards 
convergence in the saine way as an ordinary Fourier series. In parti- 
cular, it converges to i{f{x-\-0)-\-f{x~i))} iff{x) is of bounded variation 
in the neighbourhood of x. 
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Consider the integral 

27TI J 


(1.9.1) 


where 0(a;, A) is defined by (1.6.2), taken round a large closed contour 
in the A-plane. Let A = and lot the upper half of the A-contour 
correspond to the quarter-square in the .9-plane made up by the lines 

b~a \ " b~a ] 

t _ (o < a < 

b~a \ b—a / 

where s ™ or+^7; and let the A-contour be symmetrical about the 
real axis. Then (1.9.1) is equal to a finite sum of the Sturm-Liouville 
series. 

Consider the case sin a 0, sin^ 7 ^ 0. Then by Lemma 1.7 (ii) 

~ (y— «)}sinQ'-f-0{|.9|-V^?^~")} 

and “ cos^^ (6— a;)}siii/3+0(|.9|''V^^^"^^} 

on the above quarter-square; and by Lemmas 1.7 (ii) and (iii) 
co(A) ~ 5sin{.9(6— a)}sinQ:sin/3-j-0{f7^^~^*^}. 

On the quarter-square |sin{.9(/j~a)}| > and hence 

1 


1 


.9 sin{.9(6— a)]sin a sin /3 



Hence 


and 

a.(A) 


x (x,X )<f> {y ,X) ^ eos{s(/j ~x )}eoH{s(y—a)} 

ai(A) «sin{s(6— «)} ' 


/e<(!/-x)\ 

[AVJ’ 


X X 

J <f>{y,'^)f{y) dy == I 


COS [;9 (6 — .t) j cos{.9 {y — a)} 


.9sin{.9(6— a)} 


f{y) dy + 




J s""' ''^1/(2/) I ‘^y\. 


Let 0 < S < x—a\ then the last term is 


0 


+ 0 


X 

J I/to) I 

x—8 


dy\. 
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Since dXjds — 2^, this contributes to (1.9.1) 



\M\ 




,-ht 


+0 


X 

J \ny)\ 


dy 


The second term can be made arbitrarily small by choice of S; and, 
having fixed h, the first term tends to zero as n -> oo; for it is 


0 


I / J' 


{n } \)TT}(b~-a) 


^~-B(n+i)TTKb-a) 

^ 0|-^g| + 0(e-S<«+iW(6-a)). 


A similar argument applies to the other term in ^{x,X) in which 
^ ^ y ^ Altogether (1.9.1) is equal to 


f f % + 

2m J I j ssin{a(6— a)} “ 

^ a 

J 5sin{,s(o— a)} 


The first term is precisely what we should obtain in the corresponding 
problem relating to Fourier’s cosine series, and consequently it is 
equal to a finite sum of the Fourier cosine series of f{x). This, of 
course, is easily verified directly by the calculus of residues. The 
theorem on the relation between the Sturm-Lionville expansion and 
the Fourier series therefore follows. If sina =: 0 or sin^ — 0, the 
theorem can be proved in a similar way. 

In the case in which /is of bounded variation in the neighbourhood 
of the point x, it is also easy to obtain the result directly, without 
appeal to the theory of Fourier series. In the first place 


cos{s(b~x)}co8{s(y—a)} 

sm{s{b—a)} 


O 


^t{h~x)\t{]j~-a) 






whence the part of (1.9.2) with a ^ y < x—S tends to zero, as before. 
For x—S ^ y X, 


cos{5(6—a;)}cos{5(2/-“a)} 

sin{5(6-~a)} 

= — + 0(e-28<)], 
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tf 8 < b—x. The contriDOTlOH'ot tne O-term tends to zero, as before. 
The main term contributes^xo xhe ^/-integral in (1.9.2): 

X 

■?’>f 

x-h 

X X 

x~h x~h 

1 pisS n r 

X- h 

The first term is/(x — Q)/(2A), which contributes l f{x—0)y as required. 
The second term is 0(e~^^/|A|), which gives a zero limit. Also, since 
/ is of bounded variation, we can write f{x-~0)—f(y) — 9{y)~My), 
where g{y) and h(y) are positive and steadily decreasing, and tend 
to zero as y ~>x. By the second mean- value theorem 
X i 

J dy ~ g{x — S) J dy 


</(a:— 8)R J gtsCt-i/) dy 


g(x-h) 


This contributes to (1.9.1) 


J o |rW| = %(^-§)j, 


which can be made as small as Ave please by choice of 8. Similarly 
for the term involving h{y), and for the terms with y > Xy and this 
proves the theorem. 

/ k 

I / \ I ^ it \ I / \ \ 


Writing 


(l>{Xy A^^), 


i f/\ 

W iK)) 

the formulae (1.1.6)-(1.1.9) are now valid. 

If a: = a, similar results hold, and can be proved in the same way, 
provided that sin 0. If sin a — 0, 0(u, A) vanishes identically, an4 
the result fails unless /(a+0) — 0. The situation at x~b is.sii^ilar^ 
We also have the following result: 

If f{x) is any function integrable aver («, 6), and a ^ h 


^{x] A) » 2 


ST 
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This easily follows on integrating 

J Z—A 

round the same contour as before. The series is the Sturm-Liouville 
expansion of A), and it is easy to verify that 

h 

= J ^)>Pni^) dx. (1.9.4) 


l.IO. Example. Let q(x) ^ 0, a ^ 0, a o, h > 0. Then 


(f){x,X) 


sin X six 




sin{(6— .r)VA} 
~ VA 


cos y8 -L cos{(6 — x) VAjsin 


and 


6o(A) ~ - - sm(6VA)cosj3-(-cos(6VA)sinj8. 
VA 


The A„ are the roots of 


Now 


sin(6VA) 

VA 

tan bs 
66 * 


-tan^cos(6VA). 
tan ^ 

~~l7~ 


has all its roots real if — tan^ < 0 or —tan/S > 6, and two equal 
and opj)osite purely imaginary roots if 0 < —hin^ < h. Also 


AK) 


^ cos(6 VA, Jcos sin(6 VA,^)sin ^ ^ sin (6 VA,Jcos /S 
* 6VA„ 


2A„ 

6 co.s ^ cos(6VA„ 

2A,^ 




and 


— cos^cos(6VA^J— VA,jSin^sin(6VA,,) 
™ cos^cos(6VA,^)(l-l'Ay^tan2/3}. 
Hence the expansion is 


where VA^ may be purely imaginary. 



1 . 10 , 1.11 
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In the case tanjS ~ —6, we have 

co{X) = sin 6 VA — — — "J62Asin^+,.. 

so that there is a zero at the origin. In the neighbourliood of A = 0 

cf>(x,X) = — x+..., 

A) — (6— x)cosj3 f sinj8+... = ~a;co8j8-f ... . 

Hence J ? 

xcoB^ J yf(y)dy 2x J yf(y)dy 


^)(t, A) — Y R A21 — 

— ^smp.o'^A 


b‘^X 


Hence the first normalized eigenfunction is ^^{x) =■ Z^b~^x, 


1.11. An expansion involving Bessel functions. Bessel’s 
equation of order v is 


. 1 . 



( 1 . 11 . 1 ) 


and the standard solutions are J^{xs)y F^(a;5). Putting y ^ x~~^y^y we 
obtain the equation 


dx^ 




Vi ^ 0* 


( 1 . 11 . 2 ) 


Solutions of this are therefore x^J^{xs)y x^Y^(;xs). 
the form considered, with 


q{x) == 



Tliis equation is of 

(1.11.3) 


If we consider an interval a ^ 6, where a > 0, the conditions 

of the above theorem are satisfied. Let 

^o(a;,A) = xKT^ixs), = x%{xs), 

where s = VA. Then 

cooW = sx{J^{xs)Y^(xs)—Y^(xs}Jl(xs)} = (1.11.4) 

7T 

Taking a == 0 and yS = 0, wo have 

<f>(x,X) = lTT{ax)i{J^(xs)Y^{as)—Y^{xs)J^{a3)}, 
x(a;,A) = ^TT{bx}i{J^(xs)Y,{bs)--Y^{xa)J^{b8}}, 
and <o(A) = ^7ra^bi{J^{aa)Y^{b3)—Yy(aa)J,{b3)}. 
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Hence 

a>'(A) = -\7Ta^bi~{J,{bs)Yl{as)-J:{as)Y,(b8)}- 

Substituting for the F' by means of (1.11.4), we obtain 


/(A) = 


Y^(as)J l{as ) 2^ 
Jj(as) ^ TTasJ^(as) 


-JUas)Y,{ba)\- 


-iTTuib^— Jl{bs)Y^(as)~J^,(as) 
2s L 


Y,(bs)JUbs) 


Jy{bs) TTbsJ^(bs)j \ 


-11 


a Jl(as) 6 Jlibs) 

I TtZ 


n,m ^^bi Mbs) J,(m) } 

^ ’ "2s^ J,(bs)j- 


(25 J^{as) 2s Jy(bs) 

Hence if the zeros of co(X) are and VA^ = 5 ^, 

^ 2sl JAasMbSn) * 


Also 




The expansion obtained is therefore 
” sU libs,, ) 

b 

X 


b 

f l/Wl/sjyy(aSn)~^y(y«nyM^n)}f(l/) ^V- 


71=1 


A similar result can be obtained from the more general boundary 
conditions. 
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II 

THE SINGULAR CASE: SERIES EXPANSIONS 


2 . 1 . In many of the most interesting examples the function q(x) has 
a singularity at one end or both ends of the interval considered, 
or the interval extends to infinity in one direction or in both. We 
shall now consider these cases, but, in this chapter, with a limitation 
which makes the expansion still a series. The result is a particular 
case of the more general one obtained in the next chapter, but it 
seems worth while to prove it separately on account of its com- 
parative simplicity. 

We begin by considering the case in which the interval is (0, oo), 
q{x) being continuous over any finite interval. The case of a finite 
interval with a singularity at one end is quite similar to this. 

If F{x) and G{x) satisfy the differential equations 

+ ( 2 . 1 . 1 ) 

and the corresponding equation with A' instead of A, respectively, then 

h b 

(A'—A) J F(x)G(x) dx = j [F{x){q(x)G{x)—G"{x)}— 

-G{x){q{x)F(x)~-F"{x)}] dx 

h 

= - J {F(x)G"{x)-G{x)F"{x)}dx 
0 

=^W,{F,G)~W,{F,G). (2.1.2) 

If A = A' — A ™ u—iv, and G = F, this gives 

b 

2v J \F{x)\^ dx = iWo(F,F)-iW,,{F,F). (2.1.3) 

0 

Now let (fix) = ^(x,X), 0(x) = 6(x,X) be the solutions of (2.1.1) such 

^(0)==sina, ^'(0) = — cosa, \ (2.141 

0(0) ™ cosoc, 0'(O) = sin a, ] 

where a: is real. Then 

= 1I^((^,0) = sin^a+cos^a: 1. 

The general solution of (2.1.1) is of the form d{x)-\-l(f>{oc). Consider 
those solutions which satisfy a real boundary condition at a; = 6, say 
{0(6)+/^(6)}cosi3+{0'(6)+/0'(6)}sin^ - 0, 



20 


THE SINGULAR CASE: SERIES EXPANSIONS Chap. II 


where y8 is real. This gives 


I = Z(A) = 


6(b)QOt^+d^(b ) 


(2.1.5) 


For each b, as cot^ varies, I describes a circle in the complex plane, 
say CJ,. Replacing cot^ by a complex variable z, 


Z = Z(A,2:)- 


e{b)z+e'(b) 

<f>{b)z+^(by 


Here I ~ oo corresponds to 2 ; = — <^'(5)/(;4(5). Hence the centre of 
CJ, corresponds to the conjugate, z = ^'(6)/^(6); it is therefore 


Also 


w,{4>,^y 

1 <f,{b)j 'U(6) m "W)p’ 


which has the same sign as v, by (2.1.3) with F — (f>, since 


== 0. 


Hence, if ^ > 0, the exterior of corresponds to the upper half of 
the 2 :-plane. 


Since — 0'(6)/^'(6) is on (for z — 0) the radius of CJ, is 

1 




e’(b) 

w„(e,f) 

_f 

W,(0,<f>) 

<f>'(b) 


“f 



2v 




( 2 . 1 . 6 ) 


dx 


Now I is inside (7^ if l{z) < 0, i.e. if i{z-~z) > 0, i.e. if 
A_if^)+e\b) ^ 

I lHb)+6(b) m)+B(b) I ^ ’ 

i.e. if > 0, 

i.e. if > 0, 

i.e. (by (2.1.3)) if 

6 

2v f dx < iW^{e-\-hl>, 


Since WJ(<^, 0) = 1, WJ,(^,<^) = 0, etc., 

W^{e+l4,5+l^) = l-l = 2il(?). 


Hence I is interior to C* if v > 0, and 
6 


1 ( 1 ) 


f \e-\-i4,\^dx< 

J V 


(2.1.7) 
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The same result is obtained if < 0. In each case the sign of I(Z) 
is opposite to that of v. 

It follows that, if I is inside and 0 < 6' < 6, then 

r 1/3 1 ^ r 1/3 I 7J'I2^^ ^ 1(0 




Hence I is also inside Hence includes if 6' < b. It follows 
that, as 6 -> 00 , the circles Q, converge either to a limit-circle or to 
a limit-point. 

If m = m(A) is the limit-point, or any point on the limit-circle, 

r m I <> I(^^) (ey i q\ 


J \e+m<l>\^ 


( 2 . 1 . 8 ) 


for all values of b. Hence 


J \ e+m(l>\^ 


It follows that, for every value of A other than real values, (2.1.1) has 
a solution ip{x,X) ^ 6{x,X)-\-m{X)(f>{x,X) 

belonging to L'^{0,oo), 

In the limit-circle case, tends to a positive limit as 6 -> oo; hence, 
by (2.1.6), ^(x) is ^^(o^oo); so in fact, in this case, every solution of 
(2,1.1) belongs to L\0,oo), 

2.2. For a given j3, Z “ Z(A) is an analytic function of A; in fact 
it is a meromorphic function, regular except for poles on the real 
axis. For the poles of Z(A) are the zeros of 

(f){b, A)cos j3H-^'(Z>, A)sin /3, 
and this is — a7(A) in the notation of § 1.6. 

Also on the circle (7^ (if > 0) 


I \d-\-i^\^dx = 


Solving for |/|, this gives 




J ^ \e\^dx. 


\l\< -V + "t, + -yi, , - 

2vj \<f>\^ dx J \<f>\^dx J 

n V A 'A ' 
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But the above argument shows that, for a given A, the region of the 
/-plane covered by the circles (7^ decreases as b increases. Hence /(A) 
converges boundedly in any region entirely in the upper (or lower) 
half of the A-plane. Hence its limit m{X) is analytic in either half- 
plane. 

Since the above right-hand side is 0{llv) as ?; -> 0, for any fixed 6, 
it also follows that m(A) = 0(1 jv). Hence, if m(A) has poles on the 
real axis, they are all simple. 

We shall assume in this chapter that the only singularities of m(X) 
are poles. Let them be Aq, A^, Ag,..., and let the residues be rg,... . 

2.3. Lemma 2.3. For any fixed complex X and A' 

hm W{iIj{x,X),iIj{x,W)} -- 0. (2.3.1) 

. C ->-00 

Since 6{x,X)~\-l(X)(f){x,X) satisfies at a: — 6 a boundary condition 
independent of A, 

W^{/9(^r,A)+/(A)^(.r,A),%,A')+/(A')9i(a;,A')} == 0, 

i.e. 

W,[ilj(x^)+{l(X)-m(X)}cl>{x,X),i/j(x,y^^^^ - 0 , 

i.e. 

A), ^{x, \')}+{l(X)-m{X)}W,{4>{x, A), ^(,{x, A'))+ 
+{liy)-m{X')}W,{il>ix, A), ^(x, A')}+ 
+{l(X)-m{X)}{l{y)-m(X’)}W,{,f,(x, A), ^(x, A')} = 0. (2.3.2) 

Now 5 

^,{<|>ix,X),^l,ix,X')} (X'-X) j 4,ix,X)>l,(x,X')dx+Wo{<^ix,X)^ix,X'}} 

0 

= 0i^j\<t>ix,X)\^dx^^+0{l) 
as 6 -> 00 , A and A' being fixed. In the limit-point case 
|;(A)-m(A)| < 2r, = jw J |^(cc, A)|2 
so that \im{l{X)—m(X)}Wf,{(l>{x,X),ili{x,X')} = 0. 

()— >co 

This also holds in the limit-circle case, if /(A) -> m(A), since then 

b 

I \(f>(x,X)\^ dx is bounded. Similar arguments apply to the other 
0 

terms in (2.3.2), and (2.3.1) follows. 
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2.4. The following general lemma on sequences of integrals will 
be used. 


Lemma 2.4. Let f^{x) be a sequence of functions which converges in 
mean square to f{x) over any finite interval y while 

f IfMl^dx^K 
0 

for all V. Then f{x) is ^^(O, oo), and if g(x) belongs to ^^(0, oo), 

00 00 
lim r f^(x)g{x) dx = ( f(x)g(x) dx. 

*'-"*0 0 

We have 

X X 

f \f(^){^d^ = lim f ^ 

J v~*-co d 

0 0 

for every X, so that /(a;) is L2(0,oo). Now 
\(f-f^)gdx < / + / 

n ' n n ' 


+ 


[f\f-L\^dxj\g\^dx^^ 


The second term can be made less than any given e, for all v, by 
choice of X. Having fixed X, the first term tends to zero as -> oo, 
and the result follows. 

2.5. By (2.1.2) 

b 

{X'-X)^4,{x,X)^{x,X')dx = W^{^x^),^x^')}-W,{i|,{x^),^{x^% 

0 

The first term on the right is 

{cos a+m(A)sin cx}{sin a—m(A')cos a}— 

— {cosa+m(A')sina}{sina— m(A)cosa} = m(A)—m(A'), 

and, if A and A' are not real, the second term tends to zero as 6 -> oo, 
by Lemma 2.3. Hence 

m(A)~m(A') 


J ijj(Xy X)il/{x, A') dx 


A'~A 


(2.5.1) 
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In particular, taking A' = A, 

J \4,{x,\Wdx = (2.5.2) 

0 

so that the case of equality holds in (2.1.9). 

Now let be an eigenvalue, and let A' == A^^-l-^^;, «;-> 0. Then 
for any fixed X 

X 

J I V>jj{x, A') + ^(x, A„) 1 2 rfx 0. 

0 

For the left-hand side is 

J |x0(x, A')+{xm(A')+iV„)^(a:, A')— A')— ^(x, A„)}|* dx 

0 

and each of the three terms obviously contributes zero to the limit. 
Also, by (2.5.2), 

00 

J \vilj{x,X)\^ dx < \vm{X)\ = 0(1) 

0 

as V since the pole of m(A') at A,^ is simple. On multiplying 

(2.5.1) by ivjr^, making -> 0, and using Lemma 2.4, it follows that 
(f>{XyXJ is L\0,oo)y and 

J ^(x,A)^(x,AJ dx = (2.5.3) 

0 ^ 

If A tends to a different eigenvalue A^, on multiplying (2.5.3) by 
ivjr^ and making ?; -> 0, we obtain 

00 

J dx == 0. (2.5.4) 

0 

If A tends to the same eigenvalue A^, it follows similarly that 

00 , 

\m^,Xn)fdx=~-, (2.5.5) 

0 

Hence the functions 

-A>) = ri,^(x.A„) (2.5.6) 

form a normal orthogonal set. 

(2.5.3) can now be written as 

\^{x,X)>ls^(x)dx = -^. 

0 ^ \ 


( 2 . 5 . 7 ) 
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2.6. The function 0(a;, A). Let/( 2 /) be i^^OjOo), and 

X CO 

<l)(a:,A) = tp{x,X) J <f>(y,X)f(y) dy +<l>(x,X) j <p(y,X)f(y) dy, (2.6.1) 

0 X 

where ^ and ip are the functions defined above. Then, for every rr, 
A) is an analytic function of A, regular for 1(A) > 0 or 1(A) < 0. 
Also, if f(y) is continuous, 

A) == ifi'{x,X) j 4>{y,X)fiy) dy -\-cj,'{x,X) J if'{y,X)f(y) dy 

0 X 

and 

cD"(a:,A) = rM j <f>(y,m(y) dy +<A"(x,A) J Hy,X)f(y) dy + 

0 X 

+{f{x, X)<f>(x, X)—4>'(x, X)>fi(x, X)}f{x) 

= {?(;«)— A}0(a;,A)+/(a;). (2.6.2) 

Thus 0(a:, A) satisfies the differential equation suggested by §§ 1 .3-1 .4; 
and 

00 

0(0, A) = <f>(0,X) j tf,(y,X)f(y)dy, 

0 

00 

O'(0,A) = cf>'{0,X) J >p(y,X)f{y) dy, 

0 

so that O satisfies the boundary condition 

0(0, A)cosa4-<i)'(0, A)sina — 0. (2.6.3) 

If O Y(i*^, A) is the corresponding function with f{y) — 0 for 2 / > X, 
then 

X X 

<1>JS-(X,A) = e(x,X) J 4>(y,X)f{y)dy +<l>(x,X) j e{y,X)f(y)dy + 

0 X 

X 

+m{X)^{x,X) J <f>(y,X)f{y)dy. 

0 

This is clearly regular everywhere except at A “ Aq, A^,..., where it 
has simple poles with residues 


4967 


X 

j <f>{y,K)fiy)dy. 

0 

E 
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Let C denote a rectangle of the form 

< R(z) < < I(^) < Vy 

excluding all A^/s. If A is in C, 

rh ^ r 




Now O y(^,A) -^<1)(x,A) if A is not real; and as = 1(A) -> 0 


= on ils(y,X)f{y)dy 


J \4^(y^WdiJY 0{v~^) 


by (2.5.2), if A tends to a value not an eigenvalue. Hence making 
X -> 00 we obtain by dominated convergence 


<E)(:r,A) 


_L 

IttI ] Z~ 


‘^iz. 


It follows that 3>(a;, A) is analytic throughout C, i.e. that the functions 
so denoted in the upper and lower half-planes are analytic continua- 
tions of each other. If C includes a point A„, we find similarly that 
<l)(a;, A) has a simple pole at A,,, its residue there being the limit of 
the residue of i.e. 

00 00 

r„H^,K) / <f>{y<K)fiy) ^y = J My)f{y) ^y 

0 0 

2.7. The following two theorems will now be proved. 

Theorem 2.7 (i). Letf(x) and 



L{m} = q{x)f(x)-r{x) 


be L^{0,oo)] let 

/(0)cosa-|~/'(0)sina = 0 

(2.7.1) 

and 

lim W{<f>{x,\)J{x)} = 0 

(2.7.2) 

for every non-real A. : 

Then 


fix) 

(0<a:<oo), 

(2.7.3) 


the series being absolutely and uniformly convergent in any finite 
interval. 


Theorem 2.7 (ii). Let f(x) be L\0,oo)> Then 

J {f{^)Y dx== 2 4- 

J »«0 


(2.7.4) 
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Either theorem can be used as a means of proving the other. We 
shall first prove 2.7 (ii) and deduce 2.7 (i), and then give the alter- 
native procedure. 

We require a number of lemmas. 

2 . 8 . Lemma 2.8. If f(x) is any function of oo), 

00 oo 

J |0(a;,A)|2fi;a: < i j \f(x)\^ dx. (2.8.1) 

0 0 

Suppose first that f{x) ^ 0 for x > X. Then the ‘condition of 
self-adjointness ’ 

00 00 

J 0(a;, X)L<b{x, A') dx — J ^{x, A')LO(a;, A) dx (2.8.2) 
0 0 
is satisfied. For 

J {0(x’, X)m(:x, X')~-ci)(x, X')m(x, A)} dx 

- - J {<I)(a;,A)cI)''(a;,A')---d)(a;,A')0''(a;,A)}dr 

0 

= -[<D(.f, A)<D'(x, A') A')a^'(cc, A)]f. 

The integrated term vanishes at x = 0, since 

00 

^(O, A) = sin a J iPiy, ^)f{y) dy, 

0 

oo 

<1>'(0,A) = —cos a J ip(y,^)f(y) dy. 

0 

The integrated term at a; = ^ tends to 0 as ^ oo, since, if a; > X, 

A) == i/»(x,A) J <f>{y,X)f{y) dy, 

0 

X 

(D'(x,A) === 4,'{x,X) I <i>(y,X)f(y) dy, 

0 

and the result follows from Lemma 2.3. 

Putting A' ~ A in (2.8.2), and substituting from (2.6.2), we obtain 

00 00 

J <^{x,X){X<!>(XjX)~-f{x)}dx ^ J (t>(x,X)[X^{x,X)—f{x)}dXj 
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i.e. (A~A) J 

0 0 

Hence, \i\ ^ u~]-iv, v '> 0, 

GO OO 

2v J A)|2 dx < 2 J |a)(a;,A)/(a:)l dx 
0 0 


. OO GO 

2j J |<D(a:,A)|2da; J \f{x)\'^dx 


Hence J A)|2 dx < J |/(x)|2dx, 

0 0 

the required result in the restricted case. 

If now f(x) is any function of i^(0, oo), and the above functions 
are distinguished by a suffix X, 

A" A^' 

liin r 10 y(^, A)|“ 6 ^ 0 ; — f |0(a;, A)|2 

A'— 00 I ^ J 

0 0 

by uniform convergence, for a fixed X\ Now 


A"' 


GO 


J |0^.(x,A)|2dx < J |0^.(x,A)|2dx 

0 0 


< ^ J lLx(x)l^dx 

0 

X oo 

= ^/ mi^dx. 

0 0 

The result therefore follows on making first X -> oo, then X' -> co. 

2.9. Lemma 2.9. If f(x) satisfies the conditions of Theorem 2.7 (i), 
then 

<1>(^.A) = 1{/(x)+0*(x,A)}, (2.9.1) 

where O* depends on Lf in the same way as O depends on f. 
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We have 

X X 

J ^(y,my) % = ^ J [ym{y)-Viy))S{y) dy 

0 0 

X 

= l[mr(y)-<l>'iy}f(y)ro-\-l J {qiy)fiy)~riy)m) dy 

(I 

on integrating by parts twice; and the integrated term vanishes at 
the lower limit, by (2.7.1). Similarly 

J ’Piy>^)fiy) dy = J-[<A(y)/'(y)-f/«'(y)/(y)]“ + 

X 

oo 

+l J {Q(y)f{y)--riy)my) dy, 

X 

and the integrated term vanishes at the uj)per limit, by (2.7.2). 
Substituting in (2.6.1), the result follows. 

If G{x,y,X) denotes the ‘Green’s function’ 

G{x, y, A) ^(x, X)(j>{y, X) {y < x), 4>(:x, A)0(?/, A) {y > x), 

oo 

then <I)(x, A) = J G(x, y, A)f{y) dy, 

0 

00 

^*{z, A) = J G{x, y, X)L{f(y)} dy, 
and (2.9.1) is oo 

f{x) = J" G(x,y,X)[\f{y)~L{f{yy^ dy. (2.9.2) 

0 


2.10. Lemma 2.10. dff(x) is L^0,oo), 


f 4 < f dx. 

n==0 i 


This is the ‘Bessel’s inequality’. We have 


oo ^ 

J dx 

f ifix)}^ dx + 2 cl— 2 2 C„, I f{x)ifs,,{x) dx 
d n=0 n^O 

? {fix)Y dx — 2 ol 

d n = 0 


for every N, and the result follows. 
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2.11. Lemma 2.11. Let F(X) be an analytic function of X = u-{-iVy 
regular for — r ^ ^ r, — r ^ ^ r, and let 

in this square. Then 

r 


Let G{\) = (A2-r2)J?’(A). 

On the upper and lower sides of the square 

IG'WI < (|A|H»-2) — < 3rif. 

r 


On the left- and right-hand sides 

l<y(A)| < M(|A|+r) < 3ri¥. 


Hence |(t(A)| ^ ZrM throughout the square, 
nary axis 


l-P’(A)| < 


3rJf 

[A^i 


3rif ^ 


Hence on the imagi- 

3^ 

r 


2.12. Proof of Theorem 2.7 (ii). Suppose first that /(a;) satisfies 
the conditions of Theorem 2.7 (i), and also that f(x) — 0 for suffi- 
ciently large values of x. Let 

00 

T(A) =: J f{x)<^{x,X) dx, (2.12.1) 

0 

this being really an integral over a finite range. Then T(A) is regular 
except for simple poles at the points A^, where it has residues 

00 

f ^n(^)/(^) dx = cl- 

By (2.9.1) 

00 00 

'I^(A) = ^ J + X J (2.12.2) 

0 0 

and the last term is 


0 


P 00 00 

J |0*(a;, A)P dx J {/(a;)}^ dx 
*" 0 0 


by Lemma 2.8, applied to Lf, 
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Let C{R) denote the contour formed by the segments of lines 
(R—i, R-\-i) and {~R—i, *— joined by semicircles of radius R 
and centres Then 


( T(A) dX = 2ni 2 cl 

C(R) -ii<K<R 

if none of the is equal to ±^R, 

On the part of the upper semicircle in the first quadrant, we have 

A — i~\-Re^^ (0 ^ ^ ^ i^)* 

Hence the last term in (2.12.2), integrated round this quadrant, gives 

( iTT . . HR V , iTT . 

0 ^ 0 ^ hiR ^ 


A similar argument applies to the other quadrants. Hence the 
integral of 'F(A) round each semicircle tends to 


J {/M}^ 

0 

as iZ -> 00 . To prove the theorem for the class of functions con- 
sidered, it is therefore sufficient to prove that 


lim 

R-^oo 


R\-i 

r T(A) dX^O 

R-i 


and a similar result with — i? in place of R. 


Let 


XW = TCA)- 2 

R~l<Xn<R + l /I — 


Then xW is regular for R~l < A < i?+l, and 




cl < 


4R) 


R~l<Xn<R + l 

where e(i?) -> 0 as i? -> oo. Hence, by Lemma 2.11, 

|X(A)| < Se{R) 

on the segment {R—i,R-\~i), Hence 

JR+i 

lim r ;!^(A) d\ — 0. 
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Also 


R-^i 


I 


dA 

A A^ 


= 0 ( 1 ), 


since the path of integration can be replaced by a semicircle on the 
side opposite to A^, and on this semicircle the integrand is bounded. 
Hence 

R+i 


/ 


R 


I 


ft+1 


A-A. 


dA 


= 0{ 


R~1 


1 c; 


i) =0(1). 


This proves the theorem for the special class of functions. 

Now let f(x) be any function of L^(0,oo), Then we can determine 
a function /o(^} of the special class such that 


Let 

Then 


0 

00 

0 


f (/(^)- = f + 2 y|-2 2 c„y„ 

^ ^ n^-0 ’ J 71 = 0 n=0 

°° N 

> f dx — %cl 

n = 0 


^2 


and also 

00 N 

< 2 J {/(x)-/o(x)}2 dx+ 2 f [Ux)~ 2 y« dx 

< 2€+2[ J {/o(x)}2 dx — ^2 y»] • 

By what has been proved, the last bracket is less than e if N is large 
enough. Hence 

r {fi^)V dx^^ c* + 4e 

J »=o 

if N is large enough. Hence 

f dxi^f cl 

J »=o 

Combining this with Lemma 2.10, the result follows. 
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It also follows that, if g{x) is another function of L^, with ‘Fourier 
coefficients’ then 

00 00 00 

dx= \ ^ {f{x)-^g{x)Y dx — I j {f{x)—g{x)Y dx 

0 0 

n«0 71=0 

= fc^d^- (2.12.3) 

71 = 0 

This formula is also true for complex functions /(x), g{x), by separa- 
tion into real and imaginary parts. 


jf(x) 9 (x) 


2.13. Deduction of Theorem 2.7 (i) from Theorem 2.7 (ii). 
From the formula 

h 

(K—^) j dx 

= Wo{iIj(x, A), 4>(x, A„)} — W^{>li(x, A), <f>{x, A„)}, 
and (2.5.3), it follows that the second term on the right-hand side 
tends to 0 as 6 -> 00 . Hence we may take f{x) — (2.9.2), 

which gives 

J Gix,ij,X)^Jt/) dy = (2.13.1) 

0 

Also 0(x, y, A) is L'^{Q < y < oo), for a fixed x^ and A not real. Hence 


00 



a-a; 


(2.13.2) 


is convergent. 

Now let/(x) satisfy the conditions of Theorem 2.7 (i), and let 
g{x) = Xf(x)-L{f(x)}. 

If A is not real 


00 


00 


00 


J 4'{y> A)I>(/(y)} dy = ^ f{y)L{ilj{y, A)} dy = A J ^(y, A)/(y) dy. 


0 0 0 
Put A = multiply by r, and make 0. Using Lemma 2.4 

as in §2.5, we obtain 


J 'Pn{y)iAM} # == A„ J ifsn{y)f{y) dy 

0 0 
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Hence, if is the ‘Fourier coefficient’ of g{x), 

dn = (A~AJc„ (2.13.3) 

and 2 |A-AJ 2 c 2 (2.13.4) 

is convergent. 

From (2.9.2), and (2.12.3) with f{y) replaced by G(a:,y,A), it fol- 
lows that 

00 

fix) = J G{x,y,X)g(y) dy 
0 


2 


n »=0 


A A^ 


(A-A„)c„ 


n-0 

the required result. The absolute convergence of the series follows 
from that of (2.13.2) and (2.13.4). 


2.14. To prove Theorem 2.7 (i) directly, the following lemma is 
required. 

Lemma 2.14. If f{x) is 1/^(0, oo), x is fixed, and v ^ 0, 

<D(x,A) = 0(|Ali|i;|-i). 

The formula 

i ( 

j {^—y)\y—^)^''(y)dy = (^— x)2a)(x)+ J i&y—2x—^)<biy)dy 

X X 

is at once verified by integration by parts. The left-hand side is 
equal to 

f 

J {^-yf{y-^){<i{y)^{y)~>f^(y)+f{y)) dy. 

X 

Writing C* = J {f{x)Y ^ 

we have 


J (^—y)^{y—x)^{y)dy < (^— xfj \^{y)\dy 


f I 




^ dyj my)?dy^ < C{^-x)il\v\ 
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by Lemma 2.8. Similarly 

{^—yf{y—^My)^{y)dy < C{^—x)in\fix\q{y)\l\v\, 
i 

J {i—y)Hy—x)f{y) dy < C{i-x)^, 

X 

J (Qy~2x—^)<b{y) dy < W(^—xfl\v\. 

X 

Heice lOMl s; 

The result follows on taking ^ = a;+|A|~^. 



2.15. Proof of Theorem 2.7 (i). The absolute convergence of 
the series (2.7.3) is first proved as in § 2.13. This depends on Lemma 
2.10 only, not on Theorem 2.7 (ii). Now 

f 0(0:, A) dX == 2ni 2 (2.15.1) 

din) -R<X„<R 

where C(B) is the same contour as in §2.12. By (2.9.1), and Lemma 
2.14 applied to O*, 




(2.15.2) 


The integral of the last term round the semicircles tends to zero, as 
in the case of the corresponding term in § 2.12. Hence the contribu- 
tions of the semicircles to (2.15.1) together tend to 27rif{x) as iZ oo. 
Hence finally it is a question of proving that 


lim r 0(ir, A) dA “> 0 

r—¥co j . 


(and similarly for —R). Let 
Q(A) 


J?— hi 


<^n'l>n{x) 

A-A„ ■ 


Then fl(A) is regular for iZ— 1 < A < iZ+1, and 








as cx). Hence by Lemma 2.11 

f 2 (iZ+i 2 ;) — o (1) 
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uniformly for — 1 ^ v < 1. Hence 

R+i 

j £i(A) dX -> 0. 


R+i RA 

J 2 2 


R-1<A„<-R+1 


This completes the proof. 


R-l^X„<R^l Rli 

0 \ 2 |c„'An(^^)ll^0- 

ue-i^A„<ie+i ; 


2.16. Deduction of Theorem 2.7 (ii) from Theorem 2.7 (i). 

Let /(a;) satisfy the conditions of Theorem 2.7 (i), and let 

n>=0 

Then //v(x) ->f{x) uniformly over any finite interval, and 
f {/iv(«)}^ dx= '2 cl {f(x)}^ dx. 

J 71 = 0 *' 

0 " 0 

Hence by Lemma 2.4 

00 00 

f fNi^)f(^) dx= { {f(x)Y dx, 

J 

i.e. lim ^ cl = f {/(a:)}" dx. 

N-><x> n==0 ^ 

This is the required result for the special class of functions. The 
general result then follows as in §2.12. 

2.17. Theorem 2.17. If f{x) is ^^(o^oo), and A is not equal to any 
of the A^, 

<D(a;,A) = y 

n=0 ^ 

The series is absolutely convergent. This follows from Lemma 2.10 
and the convergence of (2.13.2), if A is not real; and the result, proved 
for any A, clearly follows for all A. 

Now consider the integral 

r ^(x,z). 
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and proceed as in §2.15. By Lemma 2.14 

^(x,z) _ J 1 \ 

2-a \\z\im\l' 

Hence the integrals round the semicircles tend to zero. The proof 
is then completed as in §2.15. 


2.18. The interval (— 00 , 00 ). Now consider the case where the 
interval extends to infinity, or has a singularity, at each end. We 
actually consider the case where the interval is (— 00 , 00 ), and q(x) 
is continuous for — 00 < a: < 00 . 

Let <f>{XyX), 6{x,X) be the solutions of {L-~~X)y — 0 such that 

^(0) = 0, f (0) = -1, 

0(0) = 1, 0'(O) = 0. 

Then = 1. 

By the previous theory, there are functions ^^(A) and 
regular in the upper half-plane, such that 


^i(a;,A) == e{x,X)-\-m^(\)<l>{x,X) ( 2 . 18 . 1 ) 

is L^—oo, 0), and 

</>^ix,X) = d{x,X)+mSmx,X) ( 2 . 18 . 2 ) 

is L^O, oo). Then 


W"('/'i.>/'2) = mi(A)— W2(A). 

As in §2.5 

J I Wdx= J I ij,(x. A) dx = - 

— CO 0 

(2.18.3) 

Hence l{m^) > 0, < 0, for v > 0. 


Let 

0{x,y, A) 


OTi(A)— » i2(A) 


{y < x), 


^i(a;, X)t]fi (y,X) 
mi{X)—m^(X) 


(y > x), 

(2.18.4) 


and ‘I>(x,A) = J 0(x,y,X)f{y) dy, 

— 00 

where / is the arbitrary function to be expanded. The expansion 
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will reduce to a series if both mi(A) and m^Q^) are meromorphic 
functions. We have 

6(Xy A)+m 2 (A) 0 (a:, A) 


cD(x, A) 


7/^l(A)— ^^(A) 


J {6{y,^)+mS)<f>iy>^}}fiy) dy 


+ 


+ 


d{x, A)+w^l(A)^(a;, A) 

»ii(A)— wiaCA) 


J {6(y,>^)+mS)‘l>iy,Wiy) dy- 

( 2 . 18 . 5 ) 


There are three possibilities: 

(i) At an eigenvalue A^^, mi(A^J ™ '^hi^n) ~ 

mi(A)— 7^2(^) ~ 

Then 0 (:r, A) has the residue 

00 

^Wa;,A,j+^9i(x,A„)} J {d{y,K)+H-<f>iy>K)}fiy) dy. 

— 00 

(ii) m^(X) and m2(A) both have simple zeros, 

mi(A) ~ m^(X) ~ 

Then 0 (a:, A) has the residue 


H'l 


~ — 0{x,K) f ^(y>K)f{y)dy. 

—y-i J 


(iii) ?/2 i(A) and mo{X) both have simple poles, 


mAX) ~ , rriAX) ~ . 

Then <l)(a:, A) has the residue 

00 

J^?^-^(x,XJ f <t>(y,XMy)dy. 

J 

— oo 

The theory of this case is much the same as that already considered. 

Suppose in particular that q(x) is an even function. Then (f>(x,X) 
is an even function of x, and 0 (x,X) is an odd function. It follows 
that, if 6(x,X)-i-m2(X)(f>(XjA) is L^0,oo), 6{x,X)—-m2{X)(f>{x,X) is 
0 ). Hence mj^{X) — m2(A). All eigenvalues occur under 

(ii) or (iii), and each eigenfunction is either odd or even. 


REFERENCES 

Weyl (1), (2), (3), Titchmarsh (4), (5), (7), (9), Kemble (1). 
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3.1. We now consider the same problem as in the last chapter, but 
we no longer make the hypothesis that m(A) is a meromorphic func- 
tion. All that is known so far is that it is an analytic function of A, 
regular in the upper half-plane, and that Im(A) < 0. 

We proceed formally as follows. Defining A) as before, it can 
be proved that 


~~ I* <I>(a;,A)6fA , 


-JiUS ^ 

where S > 0. It is then a question of the behaviour of this integral 
as 8 0. Now 

X CO 

A) = i/r(a;,A) J 4>(y,X)f(y) dy + ^(x, A) J ^(y,X)f(y) dy 

0 X 

X 

-= {6{x,X)+m{X)4>{x,X)} J (f>iy,X)f(y) dy + 


-t- <f>(x. A) J {e{y, X)+tn(X)(f>{y, X)}fiy) dy. 


We shall prove that 


lim r du = ^'(A), 

S— >0 X 


where k{X) is a non-decreasing function of A. Hence we obtain 
formally 


- J a)(x,A)rfA 

-Ti+iS ^ 

-i J {e{x,X)+miX)<l>{x,X)}dX j ^{y,X)f{y)dy + 

0 

^ R+ib CO 

+I -i J <f>{x, A) dX J {e(y, A)+m(A)^(2/. X)}f(y) dy 


00 00 

i J ^x,X)dk{X) ^ <f){y,X)f{y)dy 
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as i2 -> 00 , S -> 0 , since 6(x, A) and ^{x, A) are real for real A. Hence 
the expansion formula involves Stieltjes integrals, and is formally 

00 00 

^ J ^^ 0 ) J dy- (3.1.1) 

— 00 0 
oo 

If we write g{X) = J (f>(y,X)fiy) dy, (3.1.2) 

0 

00 

then (3.1.1) is f(x) = ^ J M%(^) <^'^(A). (3.1.3) 

Also 

00 00 00 

J {f(^)Y “ J /(^) J 

0 0 - 00 

oo 00 00 

= ^ J gW dk{X) J <f>{x,X)f(x) dx = ~ j {g{X)}^ dlc{X). (3.1.4) 
— 00 0 —00 
This is the Tarseval formula’. 

Now consider the interval (— cc,oo) instead of ( 0 , oo). In this case 
0(a;, A) is given by the formulae of § 2.18. We obtain 
A 

K / “ ^ ™ = ^(^)> (3 - 1 -5) 


mi(u-i-iS)~-m2{u+i8) 

mi(u~\-i8) 


lim f -1 , = 

0 

du = £(A), 


(3.1.6) 


(3.1.7) 


where |(A) and ^(A) are non-decreasing functions of A, and 77 (A) is of 
bounded variation. Now 

f R-Yih -I 




— - J (!){x, A) dX 
-n+ih 

f zR-f-iS X 
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Proceeding as before, we obtain formally the expansion formula, 

00 00 

^ J ^y + 

— 00 —00 

OO 00 

+ 1 J e{x,\)dr^(\) J ,l,(y,X)f(y)dy + 

— 00 —00 

OO OO 

+ 1 J ^(a;,A) J e(y,X)f(y) dy + 

— OO — OO 

00 00 

+ 1 J <t>{x, X) dC{X) J <f,{y, X)f{y) dy. (3.1 .8) 

— OO — OO 

Writing 

OO 00 

g,(A) == J e(y,X)f(y)dy, /t(A) = J <i>{y,X)f(y)dy, (3.1.9) 

— 00 —00 

this is equivalent to 

00 OO 

/(X) = i f ff(x,x)y(X) dm +- f dm + 

77- J TT J 


+i r<^(x,A)i7(A)d,,(A)+-- f ,^(x,A)MA)c^^(A). (3.1.10) 

J TT J 

— CO —00 

Multiplying by f(x) and integrating again, we obtain the Parseval 
formula 

OO 00 00 

J {f{x)Ydz = A J {gr(A)}2rf^(A) + ^ J g{X)h{X) dm + 

— CO — OO — OO 

00 

+ A j{A(A)pd^(A). (3.1.11) 

— 00 

In many cases these formulae take a simpler form. Suppose, for 
example, that ^^(A) tends to a real limit as 1(A) -> 0, so that 

0^(a;,A) = 0(a:,A)+mi(A)0(a;,A) 
is L\—QOy 0). We then have formally 

dr){u) ~ ni^{u)d^{u), = {^i{u)Y d^{u), 

G 
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Hence (3.1.8) reduces to 

00 00 

f(x) = i J m) J Uy,m(y) dy. (3.1.12) 

— 00 —00 

Again, consider the case where q{x) is an even function. As in 
§2.17, this involves ^^(A) = — m 2 {A). Hence 

nti(A) __ 1 
mi(A)— m2(A) 2 

and rj(A) == 0. Hence the expansion formula is 

00 

0(x.A)di(A) I 0(y,X)f(y)dt/ + 

— 00 —00 
00 00 

+ 1 j </>(xJ)dm j <i>(y,X)f(y)dy. (3.1.13) 

— 00 —00 

3.2. To justify the above formulae, we first prove two lemmas. 
Lemma 3.2. For any fixed and 

Ux 

J l{m{u-\-iv)) du (3.2.1) 

Ui 

is bounded as v-^0. 

Consider first the case of a finite interval (0, 6), with given boundary 
conditions at the ends. Let A,^ iffn,b(^) be the eigenvalues and eigen- 
functions. Let Z(A) be the function so denoted in §§2.1-2. 2, and let 
5 be the residue of l(X) at A^ Since the formulae of Chapter II 
hold in particular in this case (or may be verified directly), 

(3.2.2) 

by (2.5.7). Hence the Parseval formula gives 

J \d{x, X)+l{X)<l,ix, A)| 2 (to ^ 2 

Taking A = t, the left-hand side is bounded as 6 -> 00 , by the analysis 
of §2.1; e.g. in (2.1.6) I(Z) is bounded as 6 -> 00 , by the property of 
the circles Ct, and V = 1. Hence 

00 

V 




is bounded as 6 -> 00 . 
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Also, by (2.10.2) applied to the finite interval (0,6), and (3.2.3), 




^ n,h 


du 


Kence J " W) ^» = 2/- J 

Ui Ui 

Let —N < U 2 ^ Ny where A > 1. Then if |A^ > 2A, 

\'^-'Kb\ > iKb- Hence 

Ut U2 

r V du ^ C 4: du ^ HN ^ 16A 

J J AS., ^ AS7+ 1 ■ 


[f |A„,,| < 2N, 

Ut 


r V du ^ C vdu _< ^^^+1 

J (U — K.bf+V^ ^ J {u — X„^bf+V^ - W 7T j • 


— 00 
Ut 


[fence 


/ -I{^(A)} 


du 


s less than an upper bound independent of b and v. Since 1(A) -> m(A) 
IS 6 •-> 00 , the result follows. 


3.3. Lemma 3.3. The function 

A 

lc{X) = lim r — \{ni(u-\-ih)] dii 


(3.3.1) 


xists for all k and A except possibly for values belonging to an 
numerable set; k(X) is a non-decreasing function, and 
A A 

lim f —l{i{j{x,u-{-i8)} du — [ ^{x,u) dk(u), (3.3.2) 

8->o ^ ^ 

K K 

t follows from the above lemma and (2.5.2) that as v -> 0 

Ut 00 

jduj \ip{x, A)| * dx == 0/-J (3.3.3) 

or fixed Ui and U 2 . Now by (2.5.1), with A' = i (or any fixed value) 

00 

7 /i(A) = m{i)-\-(i—X) j iff(x,X)ip(x,i) dx 


= o|j |^i(a;,A)|2da;|^. 
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Hence 

. J |m(A)| du < 

Ui 

Hence 

is finite, and so 
exists for almost all u. 
Let 


|(M2— Ml) J |»l(A)|2dM|^ = 0{v-i). 

Ui 

1 Wa 

Jdt; J |m(A)l du 

0 Ux 
1 

J |m(A)| dv 
0 

A 

if (A) = J m(z) dz. 

i 


(3.3.4) 


(3.3.6) 


Then M(X) tends to a limit as for almost all u. Let Mq be 

a value of u for which (3.3.5) exists. Then 


( Wo+i Wo+ iv u+iv V 

j + j + j yn{z)dz, 

i Wo+i Uo^iv' 

and the first term is constant, and the second term tends to a limit, 
as 0. Hence the last term also tends to a limit. Taking imaginary 
parts, the existence of A;(A) for almost all A follows. Since — Im(A) 0, 
i(A) is non-decreasing. 

Actually the limit (3.3.1) exists wherever A;(A) is continuous, i.e. 
except in an enumerable set. For at such a A we can find r] and 
such that 

k{X-\- 7 j)—k(X—rj') < € 


and such that the limit exists at A+t^ and X—r)\ Hence 
A All; 

J —l{m{u-\-ih)] < J —l{m{u~\-ih)] du < € 


A-r 


A — T 


for 8 < Sq. Hence 
A 

J —l{m(u-\-ih)~m{u-\-ih')}du 

A ' 

A — 7] 

= J ~~\{m{u~\-ih)—m(u-\~ih')}du-]- 

Wo 

A A 

+ J —-l{m{u-\-ih)} du — J —l{m{u-\-ih')] du 

X — 7 j' X~ 7 j' 

< 3e 


if 8 and 8' are sufficiently small. Hence the result. 
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Now suppose that k is continuous at k and A. We have 
A A 

J lilt(x,u-\-iv) du = I j {d(x,u-\-iv)~{-m{u-{~iv)(f){x,u-{-iv)} du 

K K 

A A 

= J l{d(x,u-\~iv)] du ^ K{7n(u-\-iv)]l{i:l){x,u-]-iv)] du + 

K K 

A 

+ J lm{u-{-iv)K(f>{x,u-\-iv) du, 

K 

Since 6 (Xj u) and (j>(x, u) are real for real u, \d{x, u-\-iv) and l<j>{x, u-\-iv) 
are 0 {v), uniformly with respect to u over a finite interval. Hence 
the first two terms tend to 0 with v (using ( 3 . 3 . 4 )). The third term is 

A 

[Iiff(i^-f it;)R<^(a;, J \M{u-\~iv)K^J^x,u-\-iv) du 

K 

A 

-> [{k(u)-\-U}<l>(x, m)]^— I {H'u)+ 0 } 4 >u{x, u) du 

K 

A 

~ J (f>(x,u) dk{u), 

K 

The process is justified since IJI/(A) is bounded, by Lemma 3 . 2 . 


3 . 4 . Let 


x(a;,A) =- I* 4 >(x,u) dk(u). 

K 

J I J lifj{x,u-]-iv) du\ dx < K, 


( 3 . 4 . 1 ) 


( 3 . 4 . 2 ) 


where K is independent of v, if k and A are in a fixed interval. For by 
( 3 . 2 . 2 ) 

b X 

J 'f’n.bi^) dx j li/)i,(x,u+iv) du 


= ~ o( 

“’"J Ia^.,+ 1 

K 

Hence the Parseval theorem gives 
b / X \ 2 

I [ J IM^yU+iv) dttj dx=--0^ 


by ( 3 . 2 . 4 ). The result follows on making 6 -> oo. Making t; -> 0 , it 
follows that is L^( 0 ,oo). 
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Let/(a;) be jL2(0,oo). Then 


= / x(y’^)f(y) ^y 


(3.4.3) 


exists, and is a bounded function of A in any finite interval. 
Now consider 


~ J a)(a;,A)dA 

-R+ih 

[ i2+iS X 

-i J ^(x.A)dAj <i>(y,X)f(y)dy + 

-Jf2+iS 0 

p 00 

+1-^ J I 


I{^(a;,A)^(y,A)-.^(x,A)^{^,A)} = I{%,A)^(y, A)-^(x,A)%,A)} = 0(S) 
as S -> 0 (for x and y in fixed intervals). Hence (3.4.4) is equal to 


~ J A) dX J xlj(y, X)f(y) dy +0(8) 


as S -> 0 (i? fixed). Next 


1 mx, A) dX { R(/f(y, X)f{y) dy 

liih 0 

= 0(S) jduj \iff{y,u+i8)f(y)\ dy 


^^(S) J (J \^y,n-]-i8)f(y)\dy\ du 


= 0(8)1 jduj \ilj(y, ^+ 18 ) 1 ^ dyr =- 0(8*) 


by (3.3.3). Similarly 


J {R<f>{x,u+i8)—<f>{x,u)}duj li/i(y,X)f(y) dy = 0(8*). 
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Hence (3.4.4) is equal to 

R 


Now 


J J “+*S)/(y) dy 4- 0(8i). (3.4.6) 

-i? 0 


^ W A 

jduj —my,u+iS)fiy) dy = j f(y) dy j —my,u+i8) du 

K 0 Ok 

S'i(^) 

as S -> 0, by Lemma 2.4, uniformly over a finite A-range. Hence, 
integrating (3.4.5) by parts, we obtain 

p tfc 00 

i <I>{X, u) j du' J —lijjiy, u'-\-i^)f(y) dy 


K 0 

R 


-R 


— i J u) J du' J ~my, u'-\-ih)f(y) dy 

-R K 0 

R 

-^-[<i>{x,u)gi{u)fi ^^{x,u)g^(u) du. 

7T TT J 

-R 

If is of bounded variation, this is equal to the Stieltjes integral 


R 

i J <l>{x, u) dg^iu). 


(3.4.6) 


3.5. We now require the following theorem. 

Let <f>{X) be an analytic function of X = u-]-iv, regular for > 0. 
Let it be bounded on each line v — constant, and let its maximum 
modulus on the line tend toO as v -> oo. Let (f>(X) — p(u, v)-\-iq{u, v), and 

00 

J \p{u,v) \ du ^ M (v > 0). (3.5.1) 

— 00 

Then there is a function p(t), of bounded variation in (—oo, oo), such that 


m - ^ J 


dp{t) 

t—X 


{v > 0). 


(3.6.2) 


Also 


lim I p{u,v) du = p(^ 2 )“P('^i) 

t ^-►0 ^ 

Ui 


for all values of and 


(3.6.3) 
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Integrating (f>(z)l(z—X) along the straight line R~\'iy) and 

round the semicircle above it, where 0 < t/ < v, and making R~>co, 
we obtain 

R-Hy 

lim r dz. (3.5.4) 

' 27 nR-.oo J z-X ^ ^ 

-R+ty 

RHy 

Similarly 0 — lim f dz, (3.5.5) 

27r^i^->oo J z — A 

—R~\-iy 


where A' = u-\-i{2y—v). Subtracting the conjugate of (3.5.5) from 
(3.5.4), 





p(x,y) 

x-\~iy—X 


dx. 


(3.5.6) 


From this and (3.5.1) it follows that 

00 

n^—y) J 7r{v—y) 

— 00 

M 

and making y -> 0 l<^(^)| ^ — • (3.5.7) 

ttv 

Now let X 

<f>iW = f <^(2) dz = pj{u,v)-\-iqi{u,v), 

i 

A 

(^ 2 (A) = J ^i(z) dz = p^(u,v)^iq^{u,v). 

i 

By (3.5.7), ^i(A) = 0{log(l/i;)} as uniformly over a finite 

2 ^-interval. Hence (f) 2 {X) tends to a limit as v~> 0, uniformly over 
a finite i^-interval. This limit ^ 2 ('^) = ^ 

continuous function of u. 


Now 


^i{A) = 


1 u 

—ij <l>(iy) dy + j <f)(x+iv) dx, 

V 0 


1 U 

Pi{u, v) = j q(0,y) dy + j p{x, v) dx 

V 0 
1 

= j3'(0.2/)# + 

V 


and hence 
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where ^ 

x{u,v) = i J {\p{x,v)\^p(x,v)} dx, 
0 


u 

o){u,v) = i J {\p{pi:,v)\—p{x,v)) dx. 

0 

For each v, the functions x co are positive non-decreasing func- 
tions of Uy and x ^ a» < . Let 


Then 


P{n,v,h) = 


P(tt, «;,A) 


W \ h 

\ J Pi{x,v)dx 

U 


where 


Xx{u,v,h) 


1 

= J dy 4- 


uA-h u ]-h 

a>{x, v) dx. 


For given u and v, the functions xi oj^ are non-decreasing func- 
tions of h, and < M, Icu^l < M, Hence if {hy,h^-{-h^) are any 
non-overlapping intervals, 

V 

< 2 xi(«. K)}+ 

V 

+2 


Making 0, it follows that 

J |P(tt,0,A^+SJ— P(w,0,/ij| < 4i¥. 

V 

Hence P{u, 0, h) is of bounded variation, and so tends to limits as 
^->±0. Thus ^ 2 (^ 5 ^) everywhere right-hand and left-hand 
derivatives P 2 ,+(t^, 0) and ^)- ^^o 

\P{UyV,h)~P{UyVyl)\ ^ 4if. 

Making z; -> 0, 


|P('l^, 0, i^)| 4Jf-|-|p2(^+L + ^)|- 

This is bounded in any finite interval. Hence p^iu, 0) is absolutely 

4957 -rr 
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continuous, and so is the integral of its derivative, which exists 
almost everywhere. 

Further, Xi are non-decreasing functions of u, for given 

h and v. Hence, if are non-overlapping intervals, 

2 \P{u^-\-h^,v,h,)—P{u^,v,h)\ < 4if. 

V 

Making ?; -> 0, then A -> ±0, it follows that P 2 .+(w, 0) and 
are of bounded variation in (— 00 , 00 ). Let 
p{u) = 

Integrating (3.5.6) by parts, 

C» 00 

m ~ f = -. f 

TTi J (x-[~iy~Xf TTi J (x+zy—Xf 

—00 —00 

the integrated terms vanishing since Pi(x, t/) = 0(1), 

for fixed y. In the last formula we can make y 0, and obtain 


m = -. 

Ttt j (x— A)® 


(3.5.8) 


To justify this step, we observe that 

/W+i u+iVi 


a -% u-riv^ / , V 

+ J mz)(iz = 0(«)+o|iog-| 

I. u+i ' ' ' 

by (3.5.7). Hence 

Piiu, v)—p^(u, 0 ) = R J i<f>i(u-\-iy) dy = o|r|«+logij|. ( 3 . 6 . 9 ) 

Hence as 1 / -> 0 

0 f gj M+logW 


and also 


00 


-A)® (x- 

—00 

since ^ 2 ( 3 ;, 0) = 0{x), by (3.5.9) with v fixed. 
Integrating ( 3 . 6 . 8 ) by parts, 




I J (a:— A)* TT% J 


■ni J (a:— A)* 


P(^) 

(x-xy 


dx, 


(3.5.10) 


and (3.5.2) follows on integrating by parts again. 
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Since (3.5.10) is uniformly convergent, on integrating over (%,'U 2 ) 
and taking real parts we have 

Ut 00 

/»» 

\p{x) dx. 




'2 (x~Ui)^-\-V^ 

Ui —00 

(3.5.3) follows from this by the theory of Cauchy’s singular integral.f 

3.6. Let f{x), the arbitrary function to be expanded, satisfy the 
same conditions as in Theorem 2.7 (i); viz. /(a:) and Lf{x) are i>^(0, oo), 
/(0)cosa:+/'(0)sina == 0, 
and IV {f(x)y 0(a:, A)) 0 

as a; -> 00 for every complex A. 

Using the notation of §3.2, and denoting by ^^(a;. A) the function 
corresponding to 0(a;, A) for the interval (0,6), we have by (1.9.3) 


where 

Hence 

Hence, as in §3.2, 


n = 0 
h 

^'v,b "= / 

0 




7 

|I<D6(a;,A)| du < ^ |c„,6>A„.6(x)|. 

J n = 0 

Ml 

By (2.13.3), with A = i, 

^n,b ^n,b)^n,bi 

where is the Fourier coefficient of if{x)—Lf{x). Now 

f \ifix)-Lf{x)\^dx ^ J\if{x)-Lf(x)\^dx. 

„.o J J 

Also ^n,b{^)l{'^~~^n,b) Fourier coefficient of G^ft(a:, y,A). Hence 

00 

2 

n=0 




\Gb{^,y,i)\^dy 

0 

X b 

= \'lfb{^,i)\^ \ \<l>{y,i)\^dy +\<i>{x,i)\^ ^ \tPt,{y,i)\^ dy. 

0 X 

f See E. C. Titchmarsh, Fourier Integrals ^ pp. 30-1. 
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Now i/jf,(x, i) == t/j(x, i)+{l{i)—m(i)}^(x, i), 

where tp{Xy i) is L^{0^ oo), and either ^{x, i) is and l{i) -> m(i), or 

l{i)—m(i) — t>| J \(j>{x,i)\^ 

Hence, as 6 oo, 

h 00 

i) -> «). J Wb(y^ i)\ ^ <^2/ -> J \Hy, i)\ ^ dy, 


Hence 


j \Gb{^.y,iV dy ^ \G{x,y)\^dy. 
0 0 
U-i 

J m(x,A)l du < K, 


where K is independent of tio, v, and 6. Making b -> oo, then 
% -> ■— 00, 2^2 00 , it follows that 


J m^A)\ 


du < K» 


( 3 . 6 . 1 ) 


From (3.6.1) and (2.15.2) it follows that the conditions of the theorem 
of § 3.5 are satisfied by the function i(b(x, A). Hence 


<^{x 


CO 


1 f dp{t) 
X-t’ 


(3.6.2) 


where p{t) = p(t,x) is of bounded variation in (— 00 , 00 ). 

We can now prove the following theorem. 

Theokem. If f{x) satisfies the above conditions and 

A CO 

x(x,X) = j <i>{x,u) dk{n), gfiX) = J x(2/.A)/(«/) dy. 


then 


00 

f{x) = ^ J <f>{x,u)dg^ 


Consider the in^pgral J ^>(a;, A) dX 


(u). 


(3.6.3) 


(3.6.4) 


taken round the rectangle with corners at iJJ+iS, As in 

§2.16, 


R+i 


lim f <^{x,X) dX = —TTif{x). 

jR->oo J . . 





3.6, 3.7 THE GENERAL SINGULAR CASE 

Also by (3.6.2) 

c» 
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Hence 


X 

J 


RO(a:, R-{-iv) dv 


00 

;/(■ 

~QO 


arctan 


1 


R—t 


- arctan 


R — ti 


dp{t) 


\dp(t)\ 


oo 

•0 J \dp{t)\ 


= ojMMOI+o/ 

— CO —A A 

which tends to zero as R oo (choosing first A and then R) uni- 
formly with respect to S. 

On taking the imaginary j)art of (3.6.4) and making first R suffi- 
ciently large, then S sufficiently small, and using (3.4.6), we obtain 
(3.6.3). 

It is also easily seen that, under the conditions of this section, 
gi{u) is of bounded variation in (— 00 , 00 ). If sin a 7 ^ 0 , 

00 

<5(0, A) = sina J <P(y,^)f{y)dy, 

0 

and (3.6.1) with x = 0 gives 


J [ my,>^)fiy)(iy 


-CO I 0 

Uv \1 \ CO 


Hence 2 / J '')/(’/) 

I 0 

for any set of Hence 

Wv I 1 CO 

1 J dll ^ li/i(y,X)f{y) dy 


du < K. 
du sC K 




Uv 0 

and making -> 0 

1 il7K+i)-!/K)l < K. 

If sin a ™ 0 , we can argue similarly with 0^(0, A). 

3.7. The Parseval formula. 

Theorem 3.7. Let f{x) belong to L\0,oo), Then the sequence of 
functions ^ 

ffnW = j ^(y.^)/(y) % 



THE GENERAL SINGULAR CASE 


54 


Chap. Ill 


converges in mean, with respect to k{X), over (—00, 00), to a limit g(X); i.c. 

00 

lim f {g'(A)-^„(A))M^;(A) = 0; 

n~>oo •' 

—00 

00 00 

and J {/(x)}^ = {^^(A)}^ dk{X). 

0 —00 

To prove this, we require the mean convergence theorem for 
Stieltjes integrals. 

Consider first a finite interval (a, b) in which k(x) is bounded and 
non-decreasing, and let there be a sequence of continuous functions 
fni^) such that 

b 

J {fm{^)-fn{^)Y d'Hx) 0 (3-7.1) 

a 

as m “> 00, n -> 00. 

The function ^ ~ k(x) maps the interval a x ^ b on the interval 
k{a) ^ X ^ k{b) in the following way. If k{x) is continuous, to 
each X corresponds just one value of |. At a discontinuity, e.g. if 
0) < k{Xi-]-0), the point x^ corresponds to the whole interval 
0), A:(a:i+0)} (open or closed as the case may be). If k{x) is 
constant over an interval {x^,X2), every x in the interval corresponds 
to the same value of Thus an inverse function x ^ x(^), of the 
same type as k(x), is also defined. 

If the function f(x) is measurablef B (in particular, continuous), 
the function f{x{^)} = F(^) is also measurable B, The Stieltjes 
integral of f(x) with respect to k{x) is then 

b kib) 

jf(x)dk{x)= I F{^)d^. 

a k{a) 

Conversely, if F{^) is measurable B, and L{k{a), k(b)}, then 

kib) b 

j F{^)d^^ I F{k{x)}dk(x). 

k{a) a 

We can therefore put (3.7.1) in the form 

kib) 

J {FJO-Km^d^^O. 

k(a) 

t See E. W. Hobson, The Theory of Functions of a Real Variable (2nd ed.), vol. 1, 
§§ 132 and 445-8. 
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Hence there is a function F(^) of L^{k{a),k[b)], such that 

k(b) 

k(a) 

The usual method of defining this function shows that in fact it is 
measurable B, Hence f(x) — F{k(x)} is measurable B, and 

J dk{x) 0. 

a 

The case of an infinite interval can be discussed similarly. 

Now let f{x) be any function of iy2(0,oo). Then we can define a 
sequence of functions satisfying the conditions of the above 

theorem, and zero for x '> n, such that 

oo 

lim r {f(x)-f„(x)Y dx = 0. 

n—>oa J 

0 

n 

Let A„(A) = J <i>{y,X)fn(y) dy. 

0 

A An 

Then J hjn) dk{u) = J dk{u) J (f>(y,u)f„(y) dy 

K K 0 

n A 

= / fn(y) J «) dk(u) 

0 K 

n 

== If n(y)x(y>>^)dy 
0 

say. The above theorem then gives 

00 CO 

IJp^) = “ J ~ J dk{u). 

~~QO —00 

This integral converges uniformly with respect to x over any finite 
interval. Hence the process leading to (3.1.4) is valid, so that 

00 00 

J {fn{^)Y = ^ j dk(\). 

0 —00 

Similarly 

00 00 

J {M-Li^)rdx = i J {hJX)-Ki\)}^dk{X). 

0 —00 
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Hence the functions h^^{X) converge in mean square to a limit ^(A), 
such that 

00 00 

J {/(*)}^ j 

0 —00 

If G^(A) corresponds in the same way to F{x), we shall clearly have 
similarly 

00 00 

J {f(x)-F{x)Ydx = i J {^(A)-(?(A)}2rfA;(A). 

0 —00 

Taking F{x) = f{x) (x < n), 0 (a; > n), we have G(X) — ^„(A). Hence 

J {^^(A)- 9 '„(A )}2 dk(\) = ^ J {/(^)P dx-^0 

— 00 n 

as n -> 00 . Thus ^^(A) is also the mean square limit of ^^^(A). This 
proves the theorem. 


3.8. The interval (—oo, oo). The rigorous discussion of this case 
is similar to that which has been given for the case (0, oo). We shall 
indicate briefly how it proceeds. 

We define functions \fj^{x,X), \Ij^(x,X), G(x,y,X) as in §2.18. Let 
2/? Greeii’s fuiiction, eigenvalues, and 

eigenfunctions for the interval (a,b). Then (cf. (2.13.1)) 

f GaA^,y^)Ka.M dy - 


Hence the Parseval theorem gives 

J \Gafi{x,y,»?dy = 2 

Hence, as in § 3.2, 

U2 b oo , 

J d« j J 

n. OT. = n ' n,a,o/ i 


Ui a 


n = 0 
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as v -> 0, uniformly with respect to a and fc, x, and being fixed. 
Hence /i\ 

j du j \G(x,y,X)\^dy=^oitV 


Now let a; = 0; by (2.18.4) 

Q/Q y \) — (y < 0 ) 

Hence, by (2.18.3), 


iy > 0 ). 


Uz 

f I{ mi(A)-OT2 (A)} , _ ^ 

J |mj(A)-m,(A)|2^ '' ^ ’’ 

Ui 

Uz 

f l( -]d.u = 0(1). 

J lwi(A)-mi,(A)j 


(3.8.2) 


Also it is easily verified from (3.8.1) that 


/ 0'„t,{x,y,X)>f‘n,a,biy)dy 


^n,a,b(^) 


where the dash denotes differentiation with respect to x. Arguing 
as before, we obtain 


°? / 1 \ 

J duj lCr(x,y,X)l-dy=^ OA\. 

Ui — 00 ' ' 


mi(A) — m2(A) mi(A) — m2(A) 

Hence, by (2.18.3), 

Uz 

r |m2jA)|2I{mi(A)}-|mi(A)|2I{TO2(A)} , ^ . 

J |m;(A)-m,(A)12 ' ( 


Uz 

J {mi(X)-m2(X)) 


It is also easily verified that 


^ 1 1 

'' \my—mj \mj — mj‘ 


(3.8.3) 
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Hence 


Us Us 

f Ilf 1 du< fllf ^ )if 

J I Wi(A)-TO2('^)/ J 1 Wi{A)-W2(A)/ 1oti(A)— W2(A)/ 

Ui Ui 

< f f f - 0 ( 1 ). 

J l>Wi(A)-m 2 (A)/ J I \mi(X)—m^(X)\ 

Ui Ui 

The existence of the functions ^(A), r)(X), and ^(A) defined by (3.1.5)- 
(3.1.7) now follows as before. 

Let X A 

Xi(x,X) = j d{x,u)d^(u), = / 0(x,u)d7)(u), 

K K 

A A 

X 3 (a;,A) = J ^x,u) dri{u), Xi(^>^) = / ^(^.«) di(u), 

K K 

CO 

and g'^(A) = J Xviy>^)f(y) (>' = 1> 2, 3, 4). 

0 

Then the expansion formula is 

00 

f(x) = :^ J dg^(u) + 

— 00 

+d(x,u) dggiu) + (f>(x,u) dg^iu). 

This can now be proved in the same way as before. 


3.9. The spectrum. In the case of the a:-interval ( 0 , oo) the 
spectrum may be defined as the A-set which is the complement of 
the set of points in the neighbourhood of which A;(A) is constant. If 
A:(A) is constant in an interval, so are ^^(A) (§3.4). Such 

an interval contributes nothing to the representation formula (3.6.3) 
or the Parseval formula (3.1.4). 

If m(A) is a meromorphic function, the spectrum is the set of its 
poles. This is called a point-spectrum. On the other hand, an 
interval throughout which A:(A) increases steadily belongs to the 
continuous spectrum. 

In the case of the a:-interval (~oo,oo), the spectrum is the com- 
plement of the set of points in the neighbourhood of which all the 
functions ^(A), 77 (A), and ^(A) are constant. 
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IV 

EXAMPLES 


4 . 1 . In this chapter we consider a number of examples of the above 
theory. 

The simplest case is the Fourier case, with q{x) 0. If the interval 
is (0, oo), and the boundary condition at 0 is (2.1.4), then 
6{x,X) = cosacps(a:VA)-l-A"^sinasin(a:VA), 

(^(:r,A) — sinacos(:cVA)— A“^cosQ:sin(WA). 

The function ifj{x,X) ^ 6(x,X)-{-m{X)(f)(x,X) must be a multiple of 
gtxvA jf ^ since is not L^O.oo). Hence 


Hence 


-Im(A) 


w(A) 


VA 


sina— iVAcosa 
cosa+iVAsina’ 


cos^a+Asin^a; 


(A > 0), 0 


(A<0) 


and (3.1.1) reduces to 


7T J cos^a-j-Asin'^a J 

0 0 

This of course can be verified under the usual conditions for Fourier’s 
integral formula by direct consideration of the integrand. The general 
theorem of §3.6 gives the formula in a comparatively indirect form. 
Consider, e.g., the case a ^ {n. Then 

d{x,X) = A“*sin(a:VA), (/>{x,X) cos(:rVA), 

A 

/:(A) — J u-^ du (A > 0), 0 (A < 0), 

0 

;:^(a;, A) == f cos(a;Vi^)'?4~* du = , 

J X 

0 

0 

Thus we obtain Fourier’s cosine formula in the form 


f{x) = " J cos(xV' 2 ^) Ji j f{x) 
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In the case of the interval (- 

-00, oo) we have 


0(XyA) cos(x^/^), 

, .. sin(a;VA) 


nij^(X) = iVA, 



^^(A) = — iVA, 


11 

-i( ' - 1 = 

(mi(A)-m 2 (A)| 



== 

i (A> 0 ). 0 (A< 0 ) 

11 

j( mi(A)m 2 (A) ) _ 
|mi(A)-m 2 (A)J 

-rn 


= 

^ (A> 0 ), 0 (A< 0 ). 


Hence ( 3 . 1 . 13 ) gives 

00 CO 

^ ^ J J cos(t/VA)/(y) dy + 

0 —00 

00 CO 

+“ j* J 

0 —00 

Putting A — 5 ^, this gives the ordinary form of Fourier’s formula. 


4.2. The Hermite expansion. Let q{x) — a;- (— oo < a; < oo). 
We have then to solve the equation 

g-(x=*-A)2/ = 0. (4.2.1) 

Putting y ~ we obtain 


dx^ 


.2x^+(A- 

ax 



A solution of this is 
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where = exp{(— JA— |)logz}, and log 2 is real at the beginning 

of the contour. For this gives 


dx^ 


_ 1 ) 2/1 = J e-*^-i"Vi^-*(22+2a;2+A-l)(i2 

00 

( 0 +) 

= _2 r dz = 0. 

J dz 


A solution of (4.2.1) is therefore 


( 0 +) 




Since (4.2.1) is unaltered if :r is replaced by —a;, another solution 
is <^o(— ^»A). 

If we take the above contour so that K(z) ^ ~ 1 on it, then for 
^ ^ ^o(^) = 0(e-»*’+*) 

as x -> 00. It is also fairly easy to see that ^q{~x) is large for large 
positive x\ for the maximum of xz—\z^ is at 2 = 2a:, which indicates 
that the integral is roughly of the order of e^*. It follows that the 
functions 0i(a:,A), 02 (^j A) defined in the general theory are multiples 
of ^o(~^>A), <f)Q{x,X) respectively. Let 

= co(A). 

Then 

G(x, y,X) = (y^x), 

a>(A) 


Now 


w{A) 


(y > x). 


( 0 +) 




dz 


= _(l+c-’'<A)2-iA-ir(J-JA) 

for R(A) < 1, and so by analytic continuation for all A. Similarly 
( 0 +) 

^;(0,A) = - J e-l*Vi^+l dz = (l+e-’^<^)2-l^+ir(|-JA). 


Hence 

w(A) = 2^o(0,A)<^;(0,A) 


.(l^.e-^tA)2 2i-Ar(J_JA)r(|-iA) 
.( 14 - 2S-* Vr( J- ^A). 
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This has zeros at the points A = 2n-\- 1 ; they are double zeros for 
n < 0, single for n > 0 (on account of the poles of the F-function). 
But (l>o(x,X) vanishes if A = {n < 0). Hence G{x,y,X) has 

simple poles at A = 2n-\-l {n >0). Now 


Putting z — 2{z'—x), this is 


f rl^' 

J (z'-x)«+i " 


„ - AY 

2TTl\-y-] 

2"'.n! \dxj 




e-i-^H„{x), 


where //,j(x) is the Hermite polynomial of degree n. Also if 

A = 271 -\~ 1 

co(A) = -(l-c-’^«)22i-»-JMr(-w~^e) 

-^772e22i-«77iL_iL _!L ^ 

nl n\ 

Also <f)Q{-—x,2n-\-l) = A<f>Q{Xy2n-\-l), since the Wronskian vanishes, 
and A =■ (— 1)'S by considering 0. Hence the function 


<l)(a;,A) = J G{x,y,X)f{y) dy 


has the residues 


r n\ 

J ^ 


!_ 47r^ 

VS 22''(n!)2 


e-i^’HJx)e-iv"HJy)f{y) dy 


2”n! 77* 


oo 

J e-iy‘HJy)f{y)dy. 


Hence the normalized eigenfunctions are 

e-i^HJx) 

4.3. The Legendre expansion. The more general second-order 
differential equation 

a(X)^^+b{X)^+{X-c{X)}Y = 0 


(4.3.1) 
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can be reduced to the standard form as follows. Let 

r dx 

Then J + {A-y(x)}r = 0, 


where 


^(x) = 


6(Z)-K(^) 

{a(Z)}r”’ 


y{x) = c{X). 


Putting Y — r(x)y, where 


we obtain 


r(x) = 

= 0 . 


(4.3.2) 


The general Legendre equation is 




r = 0. (4.3.3) 


Here x = sin“^Z, ^{x) = —tana:, r(x) — Vseca:, and the equation 
reduces to 

~-^ + {n(n+l)+itan2a:4-i— w^sec^a:}?/ = 0. (4.3.4) 


4 . 4 . Lemma. 

r cos^ , { — 2 )^ ( 27 t )^ 

J (cos 2 ;)"*+* J...(m— I) r(f— ^m+i5)r(f— Jm— kj 

c 

if C is a closed contour including — Jtt and in but excluding the other 
zeros of cos z. 

Denoting the integral by fj^s), we have 

r / \ r cos(5— l)2;cos2:— sin(5— l)2;sin2; 

“ J —— (cos 2)”^+^ ^ 


cos(5— 1)2; 
cos^-i^; 


dz-i-^ f 

m-i J 


cos(5— 1)2:. (cos 2:)“"*+^ dz 


(w +a— |)(m+a— f — 2 )...(— w+ 3+^) 


fo(a-m). 
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/.(»-») 4 / 


co 8 (s—m)z 


2^7T^ 

r(H- 15— |TO)f (f — is+^m) 


= (— i— |5+Jm)...(|— |m)r(f— | to). 

Hence the result follows. 

4.5. The ordinary Legendre expansion. Taking m = 0 in 
(4.3.4), the corresponding basic equation is 

^+(A+Jtan2a;+J)y = 0 (-fw < a; < §n). (4.5.1) 


Putting 3 / = yVcosa:, 


d^r . dr , . 


Putting X = i—^TT, 


^+cot^|+(A-i)7 = 0. 


(4.5.2) 


A solution of this is 


COSS2J . 

(cos cos ^)* 


where A ~ s^, and C is a closed contour surrounding f and — but 
excluding all other singularities. For it gives 


= _ f 

J (cos 


\ sin ^ cos sz 

\08Z—COS^)^ 


dFY r f sin^l^— 4cos^(cos2;— cosf) , 

J (cos 2;— cos ^)* 


Hence 


, ,udY r f “4cos^ C0S2;— icos^z , 

^ + cot|— -iF= * cos szdz. 

d^ J (cos 2 — cos ^)^ 

c 

On the other hand, integrating by parts, 

7 = _ dz 


(cos 2 — cos ^^)* s 
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r |cos2{(cos2— cos^)+|sin22: COS52; , 

J (C0S2:— ^cosll^ 

c 

1 r TCOS“;: "-icos2;cos^ , 

J (cos 2 ;— cos 


Hence Y satisfies (4.5.2). A solution of (4.5.1) is therefore 

I 

<f)Q{x) — Y cos^o; = 4cos^.r f dz. 

(cos 2 ; — cos^)* 


Clearly another solution is ^q(— a:). 
If a; — 8, 8 -> 0, 


7T — 

(l)f^(x) = 4 sin^8 J 


(cos 2 :+ cos 8)1 


4sinl8 r 

J (cose 


cos < 5 ( 7 ?— ^) 


(cos 8 — cos ^)^ 


481 cos S7T 


, r 


4V28^cos57r -j-~ — 4 V 2 cos,57r811og- 

J 0 


0 

^ 4sin^8 J 


(cos 2 :— cos 8)1 


dz - 0 81 


d 

r dz 

J 


Hence all solutions of (4.5.1) arc and we are in WeyFs 

limit-circle case. We shall here merely select the point on the limit- 
circle which leads to the ordinary Legendre expansion. 

The function is a point on the limit-circle of the circles 

g(6)cot^+g^(6) 

as 6 -> \7T, where 

9(x.A) = #(X,A) = 
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Taking cotp — oo, we obtain 


Similarly ^i(A) = 
Hence for y < x 

G{x,y,X) = 


^ 2 ( 0 :, A) = 

Moy 


■ UO) 

^o(O) ■ 
i/ii(x,X) 


Moy 


M—^) 4>o{y) 

^o(o) <^o(o) 1 <f>oi—^)^o(y) 




2 ' 


^ J cos *2 " r(|+is)r(|-|7)’ 


Hence 


m = -^ J 


^«(o)^;(o) = 


cos sZy _ 2M 


= 167rsin7r(J— ^5)sin7r(J4-i»^) 

=: StT cos 77.5 = 877 cos 77 VA. 
Hence the eigenvalues are 

K=(^+iy (n = 0,l,...). 


<f>oi^,K) = 4cosix 


cos(»+|) 


^ ^ ^ J (cos 2 ;+ sin o;)^ 

0 

— 2V277COs*a;i^(— sino;) 

= {-~l)^2^j27TCos^xP^{sinx), 

and, as A -> A„, ~ . 

Hence 3>(a;,i/,A) has a pole at A^ with residue 


*77 

J -P»(sin a:)cos*y P„(sin y)f{y) dy 

iTT 

= (»+|)cosia:P„(sina:) J cos^yP^{smy)f(y)dy. 
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Thus the normalized eigenfunctions are 

(^+ cos^o; /^(sin x). 

4.6. Next let 

q(x) = — itan^aj+m^sec^aj—J {—^rr < x < ^tt), 
where m is a positive integer. The basic equation is then 

^-^ + (A+itan2a;— m2sec2a;+i)2/ — 0. 

ctx 


Putting y — Yi Vcosa;, x = 


dW^ 


pcot^ 


dY^ 

di 


-|-(A — wi^cosec^f — = 0. 


A solution of this is 


Y, = sini»»^(^^)”V, 
\d cos 

where Y satisfies (4,5.2); i.e. 

Fi = (— l)^|...(m—|)sin^"*| J 


cos sz 


dz. 


(cos 2 ;— cos 


Hence solutions are <f>Q{x) — cos*a;Ti and Hence 


^o(0) = 


c 


cos 8Z 
(cos 


dz 




and « 0 ) = (-!)-*■ J) J 

c 

2”>+in^ 

Hence 

Momo) 

22m + 4^3 

"" m^m+|57r^Pi+ k)r(f^=l^ 

22m+4^3 ^ 

” Trr(|— m+ajr(^— m— «) 

Stt^ 

~ f(|— m4-s)r(^— TO— A’)' 
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This vanishes if ~ n — i.e. A = (n—m~\~\Y !>•••)• 

Also 

( fj \m 

rf~ C o7 g ) 

This is 0 if < 2m, and if ri ^ 2m it is 

2 V2 TT sin^l P^_^(cos |). 

Putting n = m+r, it follows that the eigenvalues are r ~ 

Jf s = n— m+^+^» 

O 2 

(w— 2m)! 2(/i— m+i) 

Hence the residue is 



(_ 1 )«-18b- 2 cos^a; a;)co8* y JT-^(sin y ) . , 

87t 2(— l)^n!/(7i— 2m)! 2(71— m-j- 2 ) 


iTT 

= (»— m+ Dcos^x P^_^(sin x) [ cos*?/ P™_^(sin y)f{y) dy. 

n\ J 

-JTT 

The expansion may be written in the form 


00 

f{x)= 2 f ^oa^y PTi^i^y)f{y) dy- 


4.7. If we apply the general expansion theorem to the above 
formulae, we obtain a proof of the Legendre expansion under rather 
restricted conditions. We shall now indicate briefly how to justify 
the expansion under the same conditions as an ordinary Fourier 
series. 

We have in the notation of § 4.5 

I* ^-iaz 

where the contour surrounds the points 2 = ±^. This contour can 
be "replaced by loops round these points, each loop coming from 
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— ioo, passing over the point in the positive direction, and returning 
to ~ioo. In the neighbourhood of z ^ 

(cos 2 ;— cos = ^2 sin sin 

(^— 2)*sin*^. 

Hence the loop round this point contributes 

e-ix( 2776 - 1 '’^ 

with error 0 

Treating the other loop similarly, we obtain altogether 

U^)= -^-sin(^^+j77)[l + o|^ijj, (4.7.1) 

provided that |s|| > A. 

If |s^l < Ay take the integral round a circle of fixed radius, and 
we obtain <^o(^) = 0(|sin*^|). 

We have 

iTT ' 

+^o(^) f i>oi-y)fiy)dy . 

X i 

Substituting the leading term in (4.7.1) in the part with y < x, we 
obtain 



^cos 


1 r 

sin(~c9.T+ 2775+^77) sin(5y+-j775+ j7r)/(2/) 

)s tts " J 

~\TT 
X 

[cos{(y— x)s+-;7s4-^7t}— cos{(x+?/)s}]/( 2 /)d 2 /. 

os TTS j 


2s cos 1 
For I(s) > 0 this 


X 


-i-f 

25e-»’7» J 


-iv 


g-i(yT-x+v)s-i£nJ-f^^ 


X 

^ J e»(^-i'>®/(y) dy, 

-in 


and we proceed as in § 1.9. The result is that the expansion is valid 
under the same conditions as an ordinary Fourier series. 
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4.8. Fourier -Bessel series. This is the case where the interval 
is (0,6), and q{x) = Let A == Solutions of 



are x^J^{x8), xW^{xs), Let <f>{XyX), 6{XyX) be the solutions such that 

^(6) = 0, ^'(6)=-l, . 

0(6) = 1, 0'(6) == 0. 

Then 


= ^^^^b^{Jv{^s)Yy(bs)—Y^(xs)J^(bs)}, 


Q{x,X) — -x^bh{Jy{xs)Yl{bs)—Y^{xs)J\ 


’Abs)). j 


(4.8.1) 


If V ^ 1, this is an example of Weyl’s ‘limit-point’ case. The only 
solution of L^(0, 6) is x^J^{xs), so that 


i.e. 


tp{x,X) = e{x,X) — ^^j^S(f>{x,X), 

Jy\0S) 


(4.8.2) 


The eigenvalues A„ are the zeros of J^(6VA). Now 

e/.(6VA) (A-AJ|6A;rV;(6VAJ+.... 

Hence 7n{X) has the residue = — 2A,^/6 (negative because the 
singularity is at the lower end). Also, if A^^^ = 5^, 


Hence the normalized eigenfunctions are 

\ b ) a„biJ'AbsJ bJ'Abs„) 
and the Fourier-Bessel expansion is 

n-l Q 

If 0 V < 1, all solutions of the equation belong to ^^(0, 6), and 
we are in the limit-circle case. Consider first the case 0 < v < 1. 
Since 

Yy{z) = {Jy(;2;)cosi/77— -t7Lj,(2)}cosecv7r 
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we have 


<f>(x. A) = 

eix,x) = 


^^*-{J,(xs)J_,{bs)-J_,{xs)J,(bs)}, 

2 sin V7T 

iS” {J.MJ’_Abs)-J-A^s)Jl(bs)}. 

2 sin V7T 


The limit-circle is the limit of the circles 


6(a, A)cot (x+6'(a, A) 
^(a, A)cot (x+^'(a, A) 


as a ^ 0. Now as a; -> 0 (5 fixed) 


(4.8.4) 


J^{xs) = 


(a:s)*' 

2‘'r(l+v) 


+ 0(x'’+% 


J'^{xs) — 


va;>’-V 

2-'r(i+i/)'^ 


0(x‘'+i). 


Hence 

()(a)cota+r(a) = 

— 2 -vpYi^) “+ (I— 0(ai-''\ cot a | )+ ©(a*--) j . 


Now let 

2-’r(i-.'P "" ^ ^(T+v) 

where c is a constant. Then 


cot a — 


2^r(l+v)(i~v)a-^-^-c2-^r{l-v){i+ix)a~^+^ 
c 2-»'f (1— v)ai+*'"— 2*T(r+ 



and 

ai+^cota+(^+v)a-^+^ — 


-2^+hr{l+v) 


c2-^r{l— v)a^+^ — 2’^r ( 1 + v)a^- 


: a^-K 


Hence the 0-terms in the above expression are negligible if 
1 / > i.e. 1 / < 1. Treating the denominator in (4.8.4) similarly, 
(4.8.4) gives, when a -> 0, 


-^JUb8)-s^JL^{bs) 


Since c may have any value, m describes the circle obtained by 
varying c, and so this is the limit circle. For each value of c, m is 
an even function of s, and so is a one-valued function of A. Its only 
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singularities are poles, viz. the zeros of cJ^(b^lX)—\''J^^{h^lX). Denoting 
these by we have 

TTX^b^ 




2 sin vTT 
TTX^b^ 






2 sin VTT 

The normalized eigenfunctions are \r^\^<j>{x,\j^), where the are the 
residues of m(A) at the points A,^. For c = oo the expansion is the 
ordinary Fourier-Bessel expansion of order v, and for c = 0 it is 
the expansion of order — v. 

In the case — 0, we have 

3^o(^) = “(y+log?-){l + 0(a;% r;(a:s) = 0{x\\ogx\) 

77 y J ttXS 

as -> 0 (5 fixed). Hence 

e{a)cotoc+d'(a) = !I^[jr;(65)-Ji(65)J(y+log|Jj(alcot<.+Ja-i)- 


Trb^s 

T“ 


^ ^ I I ) j 

Y'oibs) — Jo(65)log5 j(a^ cot a+ 

77 } 

2 

TTd^f 




+ 0(a^| log a I cot a) + 0(a^| log a | ) 


Taking 


-(y+log (ffll cot a + |a-i) + A = cot a+ Jos-i) 

77\ 2/ 77a^ 

and proceeding as before, we obtain 

cJ„(fe)-{y„(6«)-(2/77)J„(6«)logi} • ^ ^ 

This is an even function of s (see Watson, Theory of Bessel Functions, 
3.51 (3)), and a result similar to the previous one is obtained. 


4.9. Direct discussion of Fourier-Bessel series. The general 
theory of Chapter II gives the Fourier-Bessel expansion under very 
restricted conditions. A direct discussion of the formulae involved 
gives the following more general result. 
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Theorem 4.9. Let f(y) be integrable over (0,6), and let v > — 
Then if 0 <i x < b, the series (4.8.3) behaves as regards convergence in 
the same way as an ordinary Fourier series. 

Apply (1.7.5) to Bessel’s equation, taking ^(6) — 0, <^'(6) ~ 1, 
instead of the conditions at a: = a. We obtain 


J^{b8)Y,(X8)-Y,{bs)JAxs) = 

775 \ | 5|2 / 

Also, by the well-known asymptotic formula for Bessel functions, 

for 0 < arga: \z\ > 1. 

We have 

0 

b 

+ / yKJAbs)Y,(ys)-Y,{b8)J,(y3))f(y)drj. 


Let 0 < S < a:, and let 

a)(.r) ::=z <l)j(:r)+(l)2(^J::), 


where 


rh 7T ^^J^{bs)Y^(xs)-Y^{bs)J^[xs) 

- 2^ JJbs) 


O 

j* y^Jv{ys)f{y) dy, 


and 02(:r) is the remainder. Now for \s\ > l/S 
^ 1/lsI 8 

J y^Jv(ys)f{y) dy = j + j 


0 1 /lsl 

l/is| 


== j y^lysl^lfiy)] #|+o| J 2/‘|^ \f(y}\ dy 

= J l2/«l*’+M/(2/)l «^2/| + o|j^ J \f(y)\dy 

/ ^ 8 / \ / ^ht \ 


= 0 (| 5 |-i )+0 


4967 
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We now integrate round a quarter-square in the 5-plane (s = 
with sides on the lines 

ba — n7r-{-^v7T-\-^7T, bt — mT-\-\v7T-\-\Tr, 

Here ^i(^) = 

and the integral of this is easily seen to be 0(;n-^) (for fixed S and x). 
In <1>2(^) we can use the asymptotic formulae throughout, and the 
argument proceeds as in § 1.9. We obtain, for example. 


J y>My,)m d,j 

8 

= _ f cos(y g - |^7T -i7r) 

J I S 003(63 — hvTT—l-n) \ | 5 |® jj 


On the contour considered 


1 + o/i-l 1 

C0S(65 — ^VTT — Jtt) [ \kl// 

and the leading term is then the same as in the case of an ordinary 
Fourier sine series. The result therefore follows. 

Similar methods apply to the other Bessel-function formulae. 


4.10. The Weber formula. Consider the case of Bessel functions 
with the interval (a, oo), where a > 0. Taking, e.g., a = 0 in the 
boundary condition at a, we have 

<f>(3:,X) = ^aixi{J^(xs)Y^(a3)—Y^{xs)J^{a3}}, 


6{x,X) = ^a^x^s{J^(xs)Y\,(<is)—Yy(xs)J'y(as)]. 

The only solution which is small as x ^ oo, for 1(5) > 0, is 
Hf''(xs) = J^{xs)+iY^(xs). 

It follows that 


^(x,A) — 6(x,X)-\-m{X)(f>(x,X) 
must be a multiple of this; hence 

—sJ(,(as)—m{X)J^(as) = i{5r'(a5)+m(A)r^(a5)}, 


m(X) 


^J 'Jfls)+iY'Xas) 

Jy{as)-[-iY^(as) 




H,(as) • 


(4.10.1) 
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Also 


I J^(as)+tY^ias)f 
_ J^{as)Y' ^(as)~Y^{as)Jy { aa ) 
^ J*(a^)+r*(as) 

_ ^ 1 

77a J^{as)-j- Yl(as) ' 

TTl 


(4.10.2) 


and K^,{z) is real for real 2;. Hence 

HW’ias) 

U<;^'>{as) 

is real for purely imaginary s, i.e. for negative A. The final result 
is therefore 

0 

00 

X J y'^{Jv(ys)Y^{as)—Y^(ys)J^(as)}f(ij)dy. (4.10.3) 

a 

If, instead of 0 (a) ~ 0, it is assumed that 

0(a)cosa:4-^^'(^)sina: “ 0, 

we obtain the corresponding formula in which Jj,(as) and Vyias) are 
replaced by 

J^(as) ^cos a -f- + J' (a5)sin as, 

n. 

tan a — 


etc. For example , it VCVi 4 . V.^, — — , 

2v — 1 

then Jy(as) and T^(a5) are replaced by J^^^(as) and Yy^^{as). 

4.1 1. The Hankel formula. This is the Bessel-function case with 
interval (0, oo). Each end is now singular. Take x ™ a as basic point 
instead of the a: = 0 of the formulae of § 2 . 18 . Then 

^(x,A) = '^a>x^{J^(xs)Y^(as)—Y^(:xs)J^{as)}, 

Jj 


d{x,X) = -aixh{J^(xs)Yl,{as)—Y^(xs)Jl,{as)}. 
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If V > 1, the solutions of 2i^(0,a) and L'^{a,Qo) are x^J^{xs) and 
x^H^^^{x8), Hence 

Ux^) = ^a^xhJ^{xs)\^'Aas)-YMs)^^^ 

X^ J^{X8) 


As in §4.10 
Hence 


a* J^{as) ' 
m2(A) — 


H<^Kas) • 


(4.11.2) 


J.K(as) 


\J„(as) 

//W(a«)| 


I{ _^J^(as)//(i)(as)| = ‘~Jl(as) {s > 0), 


2i 


rra 


= 0 {8 = it, t > 0). 


Hence (3.1.12) gives 

00 CO 

f{x) = I J xiJ^{x^X) dX I y^Jy(y^X)f(y) dy 

0 0 

00 00 

= J xij^(xs)s ds J yijy(ys)f(y) dy. 

0 0 

Of course it is easily seen directly that 

X 

^{x, s) = ~^^H\^\xs) J yiJ^{ys)fiy) dy + 


(4.11.3) 


+; 


and 


■ j.xU^(xs) J y^H^}\ys)f{y) dy 

X 

00 

-IO(a:, 5) = '^xU^ixs) J y^Jy{ys)fiy) dy (s > 0) 


= 0 (s = it, t > 0). 

This gives the result again. 

Consider next the case 0 < j/ < 1. Then ^^(A) is given by (4.8.5), 
with a instead of b, and m 2 (A) by (4.11.2). This gives 

X^ CJ^(X8) — S^^J_^,(X8) 


’piix, A) 


cJy{as)~s^''J_^,{as)* 
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Also m (X)-~m (X) = 3 ^^ 

mi(A) m,{X) ^[Ha)(as) cJ,(as)-s^-J_,{as)] 

Using the formula 

H^^\z) = i cosec v7T{e“^’'^J^(2:)~JL^(2;)} 

we obtain 

wi( ) m^{ ) - — 

Hence for A > 0, i.e. s real and positive, 

_j 1 _ Wj^^(aa){cJ^(as)—s^''J_^{as)} 


mi(X)—m^{X) 


TOi(A)— TO2(A) 2 




__ rra {cJ^,(as)~s^^J_^(as)Y 
2 c^-~2cs^^ cos VTT * 

For A < 0, 5 = where t is real and positive. Now 
m^iiat) = —.e-i’'^iKJat) 

TTl 

and 

cJ^{iat)—{itf''J_^{iat) = cei'’^%{at)—e'’^H^'’e-i'’”U_^{at) 

HpnPP ^ - (^K,(at){cI,(at)~i'^^I_^{at)} 

mi(A)-m2(A) 

This is real, and if c < 0 it is continuous. In this case the formula 
(3.1.12) gives 

/(X) = i'J- 

« 0 (4.11.4) 

If c > 0, there is also a pole at t = cW^”. Here 

_ _ 1 ^ aK^{at)c{I ^{at)—I^ ^{a t)] 

m■^{X)—m2{^) _gi/2vj 

_ „^i/2.sinv77 Kl{at)^ 

VTT 




_ -A-ci/*’ 
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Hence the residue of l/{^i(A)--m 2 (A)} qua function of A is 


2 acV^^^^^Kl{acVi'’). 

VTT 




xi K^(xc^l^^ 


' aiK^iacV^-y 

Hence we have to add to the above right-hand side tlie term 




00 

(xc^i-'’) J y^Kyiyc}’^‘’)f{y) dy. ( 4 . ] 


In the case r = 0 , m^X) is given by ( 4 . 8 . 0 ), with a instead of b, and 
m^iX) by ( 4 . 11 . 2 ), with r = 0 . Hence 

./. \\ _ cJo(a;s)-ro(^«)-f (2/7r)Jo(a;s)Iogs 


iIji(x, X) 


ai cJ(f((ts)--Y^{as)-\-(2/Tr)jQ(a^)\ogs' 


m a) = c»^ o(«'^)- ^'o(a ^ ) +(2MJ;(as)Iogs l 

^ cJo(a«)— I'o(««)+(2/7r)d'o(««)Iogs / 

__ 2 l+^c+i( 2 / 7 T)log 5 

ttu //^^^(a«){cJo(a 6 ‘)— io(a 5 ) + ( 2 / 77 )Jo(a 5 )logs}* 

Hence for A > 0 , i.e. s > 0, 

_I_^ ^ 1 {cd;,(aa) -y o(a ^)+ ( 2 / 7 T)Jo(as)log 

mi(A)— W 2 (A) ■■ - {c-f ( 2 / 7 T)logsp-f-l 

For A < 0 , i.e. s ~ it, ^ > 0 , we obtain 

1 ^ aAVaO{c/o(a O+ ( 2 / 7 r)/o(a 0 lQg t + ( 2 / 7 r)A^ o(aO } 

mi(A)— m 2 (A) c+( 2 / 77 )log^ 

This has a pole at / “ As ^ it 

aKl{at) _ 2aKl{ae-^'^^') 

Hence the complete formula is 
00 

fi^\- ^ ^H<^Joi^^)-Yo{xs) + {2/7T)Jo(xs)\ogs} ^ 

nx) _ j {c+(2/„-jKg,)i+i X 


X jy^{cJaiys)—ra{iJB)+{2ln)Jo[ys)\ogs}fi-!f)dt/ + 

0 

oo 

•f 2 a:Wo(a:c-J’^<^) J y^Ko(ye-i”<^)f{y) dy. ( 4 . 11 . 6 ) 
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4.12. Further Bessel -function expansions. Let q{x) = x 
(0 < a; < oo). Solutions of 

are 


Ux,X) == (a;-A)i/i{i(.T-A)3}, = (x-A)iArj{f(a:-A)8}, 

where (.r— A)^ is real and positive for A real, x > A. 

(i) Let a — 0 in the boundary condition at a: ~ 0. Then 

Six A1 — 

0(x X) — A ) 

Since large as x -> oo, and small, we must have 


ip{XyX) — 0(a;, A)- 


j/ri(0,A) 


i.e. 


m(A) 


i/fo(0, A) 
0^(0, A) 


^(a;,A), 


i/ro(0,Aj* 

The eigenvalues A^^ are the zeros of i/fo(0,A). Now 


^„(0,A) == (A-Aj[l./.o(0,A)j^ +... 

Hence m(A) has the residue 1 at each pole. Hence the normalized 
eigenfunctions are 

^„(x) == <f>(x,K) = 

To determine this in a real form, consider the value of i(jq(x, A) when 
X < X, As A passes above x from the real axis on the left of it to 
the real axis on the right of it, arg(ir— A) goes from 0 to — tt, and so 
arg(a;— A)^ goes from 0 to — f??. Now 

K^iz) “ 

X10I1CO 


iw Sill "gTr 
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Hence, if x < X, 

U^, A) = :^(A-x)i[Jj{|(A-:r)»}+ J_*{§(A-ar)*}]. 

Hence the are the zeros of 

•^dA^+^^-idA’)- 

Also (from the recurrence formulae) 

== ;^(A-a:)[Jj{|(A-x)8}-jr_5{|(A-:.)t}]. 

Hence the normalized eigenfunctions are 

j. ^ -V3 ./.n (a;,AJ 

^ ^A„{J,(fA|)-J_,(|A*)}- 

(ii) Similarly, if a — ^tt, we obtain 

,/, =, V 3 ^ o(a;, fij 

where the are the zeros of 

The two expansions together make up the expansion corresponding to 
q(x) = \x\ {—CO < X < oo). 


4,13. Let q{x) — —x {0 ^ x < oo), and let ot 
condition at a: — 0. Solutions of 


dhj 


dx^ 


^+(A + X)7/ =: 0 


0 in the boundary 


are (x+X)^J^{Ux+X)^, {x+XnMx+X)^ 

and the Wronskian of these is S/tt. Writing for brevity 
X U^+X)K Y = I (y+A)2, Z - l-A^ 

we have therefore 


cf,(x,X) = ^-X^(x+\nJ,{X)Y^{Z)-Y^{X)J,{Z)}, 


e{x,X) = ^X{x+XnJ^{X)Y\(Z)-Y^{X)J\{Z)}. 

Now H^\X) contains a factor 

= expj^(a;+A)jj = expj^x»|l + |^ + ...jj 
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which is exponentially small as rt: ^ oo, if 1(A) > 0. All other solu- 
tions of the equation are exponentially large, so that \fj{x, A) must be 
a multiple of Hence 

-AVi(Z)-m(A)Jj(Z) = i{^iY’^{Z)+m(X)Y,{Z)}, 
j,{Z)+iY,(Z) H\^\zy 

For A > 0, i.e. Z real, 

-\m(X) = ^MZ)nZ) -Y,{Z )J',(Z) ^ ^ I 

' jf(Z)-\-YI{Zi ttX JI(Z)+Yl{Zy 

For A < 0, let A — Then Z = Now 

7Tl 


TTl 




Hence 


Hence 


Hence 


— 2|^ — r~K^(x)—e~^'^^I^{x"^ 

= -2jlifj(x)+e-*«7i(a;)j. 

K^(x)-\~'rre ^'^'^I^{x) ’ 








2'\ 


1 




Stt 


Also 


ilxKl+TT-^m-^JSnKiiy 


2tt 


i>ix,X) = ~^X^x+XnJ^{X)J_^(Z)-J_^(X)J^(Z)}. 
Hence, if —x < A < 0, 


277 
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and if A < ~x, 

The expansion formula can then be written down from (3.1.1). 

A similar result may be obtained with an interval (— oo, cx)). We 
obtain 

X 

^(x,X) = .^.(x+A)Wi«(X) J (y+X)iHf(Y)f(y)dy + 

— CO 

oo 

+ J(^+A)J^f(Z) J {y+X)iH<^^\Y}f(y)dy, 

X 

etc. 

4.14. Let q{x) —■ (— oo < a: < oo). Solutions of 

g+(A+.--.)i,„o 

are Jy(e^)i J^y{e^), where v ~ iVA. If 1(A) > 0, the former only is 
L\—qo, 0), but both are L^(0, oo). It is therefore the limit-point case 
as a: — - 00 , the limit-circle case as a: oo. 

Since 

TTZ 

the solutions d{x), <j){x) satisfying 6(0) == 1, ^'(0) 0, ^(0) = 0, 

^'(0) = — 1 are 


e{x) = - 


m = - 



2 sm V7T 

Hence 

mi(A) = -j;(i)/j,(i), 

xf)y{x, A) 

= 0(x)+mi(A)^(x) = -J,(e-)/J,(l). 


Now consider the interval (0,6). In the notation of §2.1, we have 
-Z(A) 



4.14,4.15 EXAMPLES 83 

As 6 -> 00 , the denominator is asymptotic to 

1 )cos(e*— |v7 T— Jtt) — J^( l)cos(e*4- ^vtt— i7T)}cot 

— {«/-„( 1 )sin(e*’— |v77 -— Itt) — 1 /^( 1 )sin(e'’ + ^vn—^-n)}, 

and similarly for the numerator. Hence 
-/(A) ~ 

— J^(l)}co8 |v7rcos (e'*+A)4- { J!_„( 1 ) + J'v( 1 )}sin sin (e''+ A ) 

<4(l)}cos |v7TCOs(e''4-A)4-{J_^(l)+J„(l)}sin |v7rsin(e*-4-A) 
where cot A — e~'’cotp. The limit-circle is obtained by giving 
cot(e''-t- A) any constant value C in this formula. Take, for example, 
C 0. Then 


mn{X) == — 




Hence 


+ 


+«^vA(e®) J RvA(e*')+</-iVA(e®')}/(y) dy\ 

X 

and the expansion formula is 

“ J “(W- "" j” "!'• 

0 —00 

4.15. The Laguerre or Sonine polynomials. The Laguerre 
polynomial of order n is defined by 


LniX) 


I d y 

«T\dAj 


(e-^'A"). 


( 4 . 16 . 1 ) 


The generalized Laguerre polynomials are 

^A-“ / d 
n\ 

™ r(».+a4-i) 


L\f{X) 


^(i-r 

r! in — r)\ 


(e-^A'‘+“) 


{n—r)\ r(» — r-f-a+l) 


i—X)'^-'. ( 4 . 16 . 2 ) 


Thus L<i\X) = (-l)>T(w-f«-fl)rs(A), 

where the T"{X) are the Sonine polynomials. 
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Now iiu = then 

x~ = {n+cc-X)u. 

Differentiating n-\- 1 times, 

Hence — {dldX)^(e-^X*^^°^) satisfies the differential equation , 

^2l+(X+l-a)g+(H-lK == 0, 


and «2 = e^% satisfies 

d^u. 


^ 

dX'^' 


(■+ t ) 


dK'^~X 


Uo 


0. 


This is therefore the equation satisfied by X^L^^\X), 
Putting X = ~ we obtain 

dH^ / o . oc^—l A 9 o\ 

^n—2oL—2\u^ 


= + — 4w-2a-2j 


^^2-1- ' ™2 * )-3- 

We should therefore expect to obtain the functions 

a;a+}e-ix»jr^(a)(a.2) 

as the eigenfunctions associated with 

q(x) = (0 < X < oo). 


(4.15.3) 


(4.15.4) 


Consider then the equation 

d^y 


dx^ 


+ 




2 / = 0 . 


(4.15.6) 


(4.15.6) 


This is equivalent to (4.15.4) if A = 47i+2a+2. Hence Y = 
satisfies 


dW 1 

a-l\d7^ 

A+2a— 2 

dX^ [ 

+ X jdX^ 

^ 4X 


0 


(4.15.7) 


corresponding to (4.15.3). 

Assume as a solution 

7 = J g[z)e^^ dz 

taken round a suitable contour. Integrating by parts, 

1 r 
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if the integrated term vanishes at the limits. Similarly 
~ = J zg{z)e^<‘dz = — i J {zg'(z)+g(z)}e^‘‘ dz, 

^2 == J ^ J {zY{z)-\r2zg{z)}e^‘‘ dz. 

Hence (4.15.7) gives 

— J J {2^Sf'(2)+2zgr(z)-2gf'(2)—gi(2)+(a—l)2sr(2;)— 

— \{\-\-2oL—2)g{z)}e^^dz — 0. 

This is true if 


g'(z) _ -(a +I)2+i( A+2ct +2) __ A-2a-2_A^a+2 
g{z) ~ 2(2—1) ““ 4(2—1) 42 

g{z^ ( 2 — 2 ^i(A— 2a— 2)^— l(A+2ot+2)^ 

Since 

e-iXY ^ J g{z)e^i^-^'> dz = J g{z' dz' 

we obtain finally as solutions of (4.15.6) 

(J+) 

~ J ( 2 :— (4.15.8) 

—00 

and 

(— j-f ) 

(f)^{x,X) ~ J ( 2 — -|)KA-2a-2)(2;-[- |)-i(^+2af2)ga;*2;^2;. (4.15.9) 


Here the integrands have the value which is real and positive for 
z real and greater than and A real. 

As X -> 00 , is dominated by the part of the integral near z = ^, 
and so is asymptotic to 

( 1 +) ( 0 +) 
xi-oc j ( 2 — == J dz 

— CO —00 


= sin{|7T(A- 2a- 2)}r{J(A— 2a+ 2)}. 


Similarly 


(— ) 


=: a;i^-le-l^’‘2ie‘‘’^(^-2«-2)sin{i7T(— A-2a-2)}r{i(— A— 2a-l-2)}. 
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<j>[ ~ a:‘-“ J (z— |)«^-2“-2)2a:2e®’® dz x<f>i 


and similarly ~ — a:<^ 2 - 

Hence 

^ 8el*"<^"^“~®sin{|7r(A— 2a— 2)}sin{j7r(— A— 2a— 2)} X 

X r{i(A-2a+2)}r{i(-A+2a-2)}, 

. 8^2gli,7(A-2a;-2) 


W^(<^1.<^2) = 


We have 


a)(x,A) = 


w^(^ir<^ 2 ) 


r{i(2a+2-A)}r{i(2a+2+A)}' 


i(2/.A)/(y)dy + 


+^i(x,A) ^2(y, 


Now llW(<f>i,(f> 2 ) has poles at the points 

A„ = r|:;(4n+2a+2) (w = 0, I,...)* 
Taking the upper sign 


^n) — * jf 


(z+ |)"+“+l 


dz. 


In this case there is no singularity at z — I, so that 


(-t 

<f>i{z:,K) = a;‘-“ J 


(-1 + ) ( 0 +) 

J ( 2 + 1 )"+“+^ J Z» + “+l 


dz 


f .r-n-a-l 

Z7!(n^: J 


dz 


^ ( — 1 sin TTCX X 


r\{n—r)\ 


X T{r~n~oL)x^ 


n 


(— l)ji-r^2»i-2r^! 

r! (n—r ) ! r(ri— r-[-a+ i) 


~-2i7rn! 

T(n+a^ 
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and 

Taking the lower sign 


(-i+) 




e*«* dz = 0, 

(z+i)-" 


SO that 0(a;, A) is regular at these points. 

The residue of 1/ If (<^ 1 ,^ 2 ) at A ~ 4n4-2a+2 is — r(/j+a+l)/27r^n! 
Hence the residue of 0(a;, A) is 


r(n+oc+l) 27 r^n\)^ 






00 

J y-+ie-iy‘L^-\y^)f(y) dy 


= f y'=‘+h-iy' L\f^y^)f{y) dy. 

l(n+a+l) J 

0 


4.16. The ‘hydrogen atom*. Let 

q{x) 


X 


(4.16.1) 


where r is zero or a positive integer. Then 


d^y 

dx^ 




This is the celebrated equation from which physicists have been able 
to derive the theory of the hydrogen atom. 

Putting y ~ x~^Y , A — we obtain 


^ - 0. 

dx^ X dx \x j 

Assuming Y == f g(z)e^^ dz 


(4.16.3) 


and proceeding as in the last section, we find that Y satisfies 
(4.16.3) if 

g\z) __ c— 2 (r+l)i 2 ; _ k—r~l A:4-r+l 
g{z) 2^+52 z-{-is z~is 


where k = \icl8. Hence 


g(z) — 
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Hence solutions of (4.16.2) are 




taken round suitable contours. 


Let 






( 2 — dz, 


(4.16.4) 


where C is a closed contour surrounding both is and —is, and where 
arg( 2 : 4 -^* 5 ) and Q^vg{z-~is) both go from — tt to tt round the contour. 
Thus (2:+i«5)^ is multiplied by and {z--is)~^ by and the 

whole integrand returns to its original value. 


(is+) 


Let 


A) 


J {z-{'isY'''^''\z—is)-'^^'^''^e^^dz, (4.16.5) 


where the loop includes is but not —is. Here 2 iXg{z—is) goes from 
—rr to TT, and, for s in the first quadrant, dbj:g{z-\-is) starts with the 
value 77 and returns to rr again. 

Putting z = is-\-^lx, 


( 0 +) 


where arg(l+^/ 2 i 5 a;) -> 0 as ^ 0 . Hence by Whittaker and Watson, 

Modern Analysis, § 16.12, 


<I>2M == 


The integral in clearly represents a function regular near ir = 0 . 
The coefficients of all powers of x up to x^^ vanish, as is seen on 
expanding the contour to infinity. The coefiicient of is 


1 


(2r+l) 


i J {z+is)'^-’'-^ 


(z—is)~^-^~h^^^^ dz 


27ri 

(^+lj\ 


(again expanding the contour to infinity). Since 
small 2 :, it follows that 


<f>i{x,X) 


(2r4-l)! (— 2i5)^+^ 


^r+l fQp 


Let ca(A) ~ TP(^i,</> 2 )* Denoting the integrals in (4.16.4) and (4.16.5) 
by /i and ^ x-^W{I„ I^). 
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Hence (2r) ! o>(A) = W{I^, Z^) j 

As above we have ^ for v ~ 1 ,..., 2r, and = "Zni, 

Hence 

(2r)!cD(A) = — 277^/2(0) 


(is+) 

= — 27ri J (z-\-is)^-^-^{z~is)~^~^-^ dz 


— CO 

( 0 +) 

= -2TTi j dC 

— 00 


Let arg(i5) ~ 9 (In < 0 < tt). We may turn the contour into a 
position in which ^ comes from infinity with arg^ = 6 — 27r, passes 
round the origin, and returns to infinity with arg^ = 9; and 
SiYg{^-\-2is) ~ 9, Hence the last integral is equal to 


uu 


Hence 

Now 


.(A) = 


- 27Ti(2*s)-2'--i ( 1 — e 2 »E*) 


(r—k)(r—k— l)...{—r—k) 


X 00 

j <f>iiy)f{y) dy J <f> 2 iy)fiy) dy. 

0 X 

This has poles at the zeros of co(A), i.e, at 


k ^ r+Ti+l, 


Now 


^'(AJ = 


A = A„ 


C2 

4 (r+w+l )2 


{n = 0, 1,...). 


— ^TT^{2is dk 

( — n — !)...( — n — 1 — 2r) dX 

477^w! / 7i+r+l\^''‘*'^2(w+r+l)® 

^-|-2r-|-l)!\ c / 


— _ (»+>•+ 

(w-}-2r+i)! c2r+3 

Also 

<l>i(^,K) = <f> 2 i^>K) = j* dz 

= x'’+ie-‘'®'2(«+r4i) f (^-f dC, 


4957 
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where z = is^-\-^jx, and the contour is now a closed curve round 
the pole. From (4d5.1) and Cauchy’s formula for the^th derivative 

£,(„ = f f 


27riJ (z—l)^-^^ 


: f' 

27riJ 


Hence 




Ip+i 

(p!)^ f (C+O^-’^e-i 


dl 


27ri(p—m)l j 


dc, 


where denotes {djdt)”^Lp{t). Hence 


J \ _ ~r+lg-cj/2(m+r+l) 277m! j-2r+i / 


t)- 


Hence 0(a:,A) has the residue 


^2r -f'3 


n\ 


2(n 4-^+1)^^'*'^ {(7i4-2r4-l)!}^ 


,r+lp-cxKn^r+l)r.2r+l / 


t) 


X 


X 


00 

J 


;-ci//(H+r+l)i2r+l^^ 




The sum of these terms is the contribution of the part of the spectrum 
with A < 0. 

We can write 

— CO 

(is-f ) 

4--™ J dz. 


where in the first integral arg(2;4-i^) goes from —tt to tt, arg( 2 ;--i 64 
from — TT to tt; in the second, arg(2;4-^«5) goes from tt to tt, arg(^— i^) 
from —TT to 77. The second term is thus ^ 2 ( 0 ;, A). Now if in 
we increase args by 77, we get an expression apparently the same as 
the above first term, but with ^vg{z—is) going from 77 to 77 instead 
of from —TT to —77. Hence 

Also a>(Ae2^'’^) = — e-2^^*a>(A). 

Hence for A real and positive 

cu(A) cu(Ae^’^^) 
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Hence 


I<D(x,A) = l{cl)(x.A)-«D(x,Ae2”^)} = f UyMiv) dy, 

0 


where 


ico(X) — 


Hence the interval 0 < 


27t(1— 

(45^+6^).. . (4r%^ H- c^)c 
A < 00 contributes a term 


00 


1 

477 ^^ 


I 


(4s2+c2)...(4r%2-j-c2)c 


s^)2s ds 


J <l>i{y,s^)fiy)dy. 

0 


4.17. Hypergeometric functions. Consider the equation 
X{X+l)^^ + {y+{<x+l)X}^ + X*Y^O (0<X<co). (4.17.1) 
Solutions are 

= F(a, b, c, -X), - X^-^F{a~-c+ 1, 6~-c+ 1, 2-c, -X), 

where c ™ y, a-\~b = ocy ab = A*, 

say a = 4A*)^ b ™ -|a— 4A*)^. 

We transform into the standard form as in §4.3. Putting 

X - sinh2 
rlW dV 

we obtain ^ ^ (4.17.2) 

where 

__ y+(a+l)X— (X+l) _ acoshx+2y— 1 — cv 
“(X(X+1)P ' “ sinh:r 


Putting Y := 

r{x)y, where 





\i+ia-y 



’■<*'= Wi, 

1 sinh“*“:r, 


we obtain 


q*{x)}y = 0, 

(4.17.3) 

where 




q*(x) = W^{x) + iP'{x) 

COsh^ird-2(a — l)(2y — 1 — a)cosha:+(2y — 1 — a)^ — 2a 
4 sinh^.r 
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Finally, put A A* — and 
q{x) = q*{x)-lo? 

_ 2(a-])(2y-l-a)coshx+2a2_4yaf(l-2y)2_ 

4 sinli^o; 

we obtain + (4.17.5) 

In this notation a ^ Ja+iVA, h = Ja— iVA. Corresponding to and 
Y2 we obtain solutions of (4.17.5) asymptotic to 

Axy-K A{x-)^~yxy~^ Ax^-y 

SiS X -> 0. Hence the origin is of limit-circle type if 0 < y < 2, and 
otherwise of limit-j3oint tyi3e. To take a definite case, let y > 2. 
Then the only solution of (4.17.5) belonging to L^(0, 1) is 

y ™ r}^{XyX) ~ {r(;r)}-^jF(Ja-biVA, ia—iVA, y, — sinh^l-x’). 

A solution of (4.17.1) for X >1 is 

X"^F(6, l-c+6, l~a+6, -1/X). 

This leads to 

y — r)2{x,X) = {r(a:)}“^ sinh“2^ |a;F(6, 1— c-f6, 1— a+6, — sinh-^l^jr) 

as a solution of (4.17.5). 

As X -> 00, 

rj2(x,X) ~ sinhl“:rsinh~2^ lx ~ 

Hence rjo is L“(l,rx)) if 0 < arg VA < in. 

The other solution, obtained by interchanging a and 6, gives 
and so is not L^. Also 


- {r{x)]~W{F{a, b, c, -sinh^ ^x), 

ginh-2ft lo; jP(6, 1— c+6, 1— a+6, — sinh~2|a:)} 
== {r{x)}~%X{X+l)}iW{F(a,b,c, -X), 

X~f>F(b, l-c+b, l-a+b, -1/X)} 


{X(x+i)}i r{c)r{b~~a) 

{r(x)Y r(b)r(c-a) 


W{X-^F{a, 1— c-l-a, 1—6-fa, —1/X), 


X-^F(b, 1-c-f 6, l-a-f6, -1/X)} 


^„,, r{c)r{b-a) _ _^2“+iiVAr(y)r(-2iVA) 

r(6)r(c-a)^ ’ r(^a-iVA)r(y-i(x-iVA) 
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as X -> 00 . Since W r/g) independent of x, this is its value. Hence 

Ate - r(Ja-WA)r(y- la-WA) 
a; _ 2a+iiVA r(y)r(-2iVA) " ^ 

X 00 . 

X %(a:, A) J r)i(y,X)f{y) dy -\-yi(x,X) J yi(y,X)f(y) dy\. 

1 0 X ' 

When A is real and positive, is real, and since a and b are con- 
jugates 

J ( rjoix, A) ] 1 

X 1-“+''. -i/J)- 

I r(«)r(6)r(c-«)f (c-~Vr„ . „ 

2“+V(cc)i(a-6) r2(c)r(6-rt)r<a-6)^ h ^ 


yi{x,X). 


1 r(.\a -l-iVA)r(y — la-\-i\/X) 

"= -2“|2va ‘ r(y)r(2WA) 

Hence there is a continuous spectrum from 0 to oo, wliich contributes 
to the expansion of f{x) - , 


Qa+2, 


oo 


r(ia+«VA)r(y-J-a+iVA) 

r(y]rT2»VA') 


Vii^. 


X)^^ f 


’?i( 2 />A)/(y) dy. 


There may also be a finite number of poles on the negative real axis. 


__ 1 


-hoc—n 


h or — ^VA ™ |o 


/—n. 


These will be at the points — iVA 
if these are positive for any positive integer values of n. The corre- 
sponding residues are easily calculated. 


4.18. The case y — h In this case 


q{x) 


a(l — a) 

4 cosh^ ^x ’ 


(4.18.1) 


and there is no singularity at a: 0. The ir-interval may be taken 
as (— 00 , oo). We have 

r(x) — 2”^“cosh““|.r, 
and solutions of (4.17.5) are 


0(^x) = cosh^^o; jF(ia+iVA, Ja— iVA, i, — sinh^^a;), 

(^(x) ~ — 2 cosh^ |a;sinh ^x i+ 2 ^ — f, — sinh^ ^x). 
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These are even and odd respectively. The solution which is small 
as a; -> 00 is 


ginh-2b 

2^“r(x) 



1— a+fe, 


sinh^ ^x/ 


r(>|r(i-2WA) _ 

r(l-|a-iVA)r(i+|a-WA) ^ ' 


|r(--i)r(i-2WA) 

r(i-ia-i\/A)r(|a-iVA)’^^ ’ 


Hence, in the notation of § 2.18 or § 3.1 

r(l — ia — iVA)r(.5-f- — fVA) 

- -"r(r- ;„-iVXjr(i.-iVA) • 

nd, since q(x) is even, ^^(A) — — m 2 (A). In the notation of §3.1 
f (A) = P{l/m2(A)}, r(A) - 

There is a continuous spectrum from 0 to oo, which contributes 


1 

77 


J 0(a;,A)r(A)rfAj %,A)/(2/)# + 

0 —00 



77 


00 00 

J <f>(x, A)r(A) dx j (f>(y, X)f(tj) dy 


to the expansion of/(a:). 

There may also be a finite number of poles on the negative real 
axis. Suppose, for example, that a — n-\-\, where is a positive 
integer. Then m 2 (A) has zeros at the points 


— iVA ~ \n—r {r ~ 0, 1,..., [\n\), 
the residue of l/m2(A) being 


A {2r-n)T(n-r+\) 

" r! ^(i-r)^(?^-r+l)* 

This is case (ii) of §2.18, with ^2 == H'l “ — iMr* Hence 0(a;, A) 
has the residue 

00 

—\A^e{x,x;) J d{y,x;)f{y) dy, 

— 00 

where 

0(x,\) = cosh^+^ — sinh^p). 
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Actually the solutions of (4.17.5), where q{x) is (4.18.1), are associated 
Legendre functions with argument tanh \x. In fact 

e{x,\) = (-l)»-2« [1^-^ Pr2'-(tanli lx). 

This is easily verified by showing that the two sides satisfy the same 
differential equation, and are equal, with their first derivatives, at 
a; = 0. 

There are also poles of at the points 

— iVA = (r =- 0, 1,..., [^n—^]), 

the residues being 

B ^ (— I)" (2r-»i+l)r(w+|-r) 

*■ r! r{-l-r)r(w-r) ’ 

This is case (iii) of §2.18, with pg — ~ Hence 0(a:,A) 

has the residue 

00 

lB,<i>{xX) J HyX)f{y) dy, 

— 00 

where 

(f)(x, A') = — 2 coslP^+^ |a:sinh hx F{r-{- 1, n— r+|, f , — sinh^ ^x). 

It is easily verified as before that this is equal to 

Another interesting particular case occurs when the coefficient of 
cosh a; in (4.17.4) vanishes. The solutions are then expressible in 
terms of associated Legendre functions with argument coth },x. 

4 . 19 . Another formula involving hypergeometric functions arises 
as follows. Consider the equation 

riW flY 

X(l+X)|^^+{c+(a+6+l)X)“£ + a67 = 0 (4.19.1) 

satisfied by F — F{a,b,c, —X). Putting 

= x*«(i+x)«“+*'+i-‘^)r 

ilW 

we obtainf j^-\-ZYi = 0, 


f See Forsyth, A Treatise on Differential Equations, (4th ed.), § 116. 
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where 

y 1~(1-— c)2 \~{c—a~hY (1— c)2— (a— 6)2-[-(c— ci~6)2— 1 

== 4:t2~ + “loTf Z)2 ^ ^ 4Z(1+Xr~ ’ 

Putting y ~ -X" — e-^, this gives 

g+(X2Z+J)y==0. (4.19.2) 


Now 

where 


X^Z-\ = A- 


AX 


BX 


X—q{x), 


X+l (X+l)2 
l(a_6)2_|(l_c)^ B = i_i(c-a-6)2, 
A = 


Writing 1 — = (1 — 2a)^, the relations between the parameters are 

a = (x-{-i^/X-J-^^J(X~A), b == a+iVA— i^(A — ^), c = l + 2iVA. 
Solutions of (4.19.2) which are L2(— oo, 0) and L^(0, oo) respectively 
for 1(A) > 0 are 

and 

cf>^{x) - l-c+6, l-a+6, -1/X). 

The resulting expansion formula is of the type (3,1.8), but does not 
seem to reduce to a very simple form. 

If .4 = 0, the q{x) of this section is the same as that of § 4.18. The 
two sets of formulae are connected by the relation! 

(1+^)2^ P(2P,2P+1-Q, Q.z) - F{P,P-Vh QAzKl+zf}, 
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THE NATURE OF THE SPECTRUM 

5.1. The main object of this chapter is to determine how the 
spectrum depends on the function q{x). 

Take the interval to be (0, oo), with no singularity except at 
infinity. Then there are, broadly speaking, four different cases. If 
q(x) -> 00 , there is purely a point-spectrum; the examples q{x) = x^, 
q(x) — a; of § 4.2 and § 4.12 illustrate this. If q(x) 0, there is a con- 
tinuous spectrum in (0, oo), with a point-spectrum (which may be 

null) in (— 00 , 0 ); the formulae of §4.10 and §4.17 are examples. If 

00 

q(x) —CO, but so that J \q{x)\-^ dx is divergent, the spectrum 

extends continuously from — oo to oo, as in the formula of §4.13. 

00 

Lastly, if q{x) — oo, and J |^(x)|~^ dx is convergent, there is a con- 
tinuous spectrum in (— oo, 0), and a point-spectrum (which may be 
null) in (0, oo). This is illustrated by §4.14. 

Actually in each case, except the first, we have to impose other 
conditions, so that classification is by no means complete; but all 
ordinary examples come under one or other of our theorems. 

5.2. We begin with the second of the above cases, but actually 
assume, instead of q{x) 0, that q(x) is L{0, oo). 

We require the following 

Lemma 5.2. Letf{x) ^ 0, g{x) ^ 0, and let f{x) be continuous, g(x) 
integrable, in 0 ^x ^ X. Let 

X 

fix) <c+ j fmt) dt (0 < * < X). (6.2.1) 

0 

Then f(x) ^ Ce^ (0 < ir < X), (5.2.2) 

Let 

X 

y = j dt, 

0 

2 -/(*)»<*)• 


4957 


O 
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Multiplying (5.2.1) by g{x), 


^ < Cg{x)-{-yg{x), 


dx ^ 

Integrating over (0, x). 


d{ 

^ye 0 


- f 0{t)dt 

< Cg(x)e 0 


-)»«)<« / 

1 J oiDdi \ 

t/e “ ^ C 1 — e ® !, 

2/<(7je» -ij 


and (5.2.2) follows. 


5.3. Let (l)(x) = <f)(x,X) be the solution of 
^^+{X-q(x)}y = 0 

with ^(0) = sina, <^'(0) = ~cos«. By (1.7.1), with a =» 0, A = 

j. V . sinsx , 

0(a;) = cos 50; sin a cosa+ 

5 

X 

+ j J sm{s(x-y)}q(y)(l>{y) dy. (6.3.1) 
0 

Let s == o-{-ity t ^ 0, and write temporarily (l>i(x) = (f){x)e-^^. Then 

j / \ ir • ,„sin5o: , 

d)^{x) ~ cos 50; sin a — cosa+ 

5 

X 

0 

Since |cos5o;| < |sin5o;| < it follows that 

X 

0 

Hence, by the lemma, 

\ 

J \i(y)\ dy\. 

Since q(y) is L(0, oo), it follows that <^i(o;) is bounded for all x, 

l5| > p > 0, ^ > 0. 
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Now consider real positive values of s. Then is bounded for 

3 p. Hence ( 5 . 3 . 1 ) gives 

... . sin 50 : 

(P(x) = cossxsmoL cosa4- 


OO .CO , 

+ ^1 sais{x-y)q{y)^y) dy \q(y)\ dy 


= fl{X)Q08SX-\-v(X)8mSX-\-o(l) 


( 5 . 3 . 2 ) 


(as x-> co), where 

OO 

/i(A) = ^(s 2 ) = Sina— i J sin sy q{y)<f>(y) dy, 


v(A)= f cossyqiy)(f>{y)dy. 

s s J 
0 


Since the integrals converge uniformly, p. and v are continuous func- 
tions of 5 . 

Similarly if 0 ( 0 ) == cos a, 0 '(O) == sin a, 


where 


d{x) ™ p^(X)cOSSX-\-Vi(X)Bm8X+o(l)y 


OO 

HiiX) = cosa — J J smsyq{y)d(y)dy, 


( 5 . 3 . 3 ) 


Vj(A) = r cossyq(y)(f>(y) dy. 

s s J 
0 


Also, differentiating ( 5 . 3 . 1 ), 

OO 

(f}'{x) ^ — 5 sin 5 o:sina— cos 5 o;cosa+ J cob s(x—y)q{y)(f)(y) dy. 

0 

Applying this in a similar way, we obtain 

^'(x) = ~ 5 /Lt(A)sin 5 o;+ 5 jA(A)cos 5 o;-f o(l). 


Similarly 

0'(x) — —Spi{X)BmSX+Sv^(X)cOBSX~\-o(l). 

Hence 

W{(f>,0) = ir{/x(A)cos 5 o:-f v(A)sin 5 o:, p^(X)coBSx-\-vi{X)Bin8x}-\'0(l) 

= s{/i(A)i'i(A)— /ii(A)v(A)}+o (1). 
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Since 0) = 1, it follows that 


/.(AK(A)-^i(AMA) = ^. 


(6.3.4) 


Now consider complex values of s. For a fixed positive t, (5.3.1) 
gives 

(f)(x) = I sin O' -J r- cosa+0(e”^^)-- 

Ji%S 

X . X . 

% + o| I* e-^^^-y^\q(y)<f>{y)\ dy\ 

0 ^0 
as a; -> 00 . Since (f){y) == 0{e^^), the last term is 

o[ f e«*i'-®)|gr(2/)| dy\ 


e«a:-2S) J \q(y)\ J \q{y)\ dy\ = o(e'*). 


Also 

J e-^‘^^-y\{y)(f>{y) dy = ( 

Hence 

^x) = 

where 

if (A) = |sma + -— 
2*8 

Similarly 

e{x) = 

where 

Ti/T/w 1 sin od 

= |cosa — 


(5.3.6) 


(5.3.7) 


oo 

j* e^^yq(y)e{y) dy. (6.3.8) 


Since 0(a;)+m(A)^(a:) is L^^q^oo), it follows that 

m(A) - -ilfi(A)/if(A). (5.3.9) 

As s tends to a real limit, the numerator and denominator in (5.3.9) 
tend to 2 /^i(A)+|ii^i(A) and 2 i^(A)+ J^V(A) respectively. By (5.3.4), 
/x(A) and v'(A) cannot both vanish for any positive A. Hence 

and the imaginary part of this is 

m(A)vi(A)— fJ.i( A )v(A) ^ 1 

fM\X)+v^{\) AV(A)+>'2(A)}’ 
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A 

m - J 


du 


(5.3.10) 


Hence the spectrum is continuous in (0, cx)). The contribution of 
this part to (3.1.1) is 


0 0 

0 0 

Next let s -> il, where I is real and positive. The integrals in 
(5.3.6) and (5.3.8) converge uniformly, and hence M(A) and ilfi(A) 
tend to the real limits 

00 

J e->yq(y)^(y)dy, 

0 

CO 

icosa + ^-4-1 J e-‘yq(y)d(y)dy. 


Hence I(m) 0 except possibly at the zeros of ilf(A). Since ilf(A) is 
an analytic function of s, regular for l(s) > 0, the zeros are isolated 
points. Hence there is a point-spectrum in —oo < A < 0, which is 
bounded below, since ilf(A) ~ |sina or cos ocJ2i8 a,s s -^oo. 


5.4. A transformation of the basic equation. In order to deal 
with cases in which q{x} is large at infinity, we make a transformation 
of the equation 

^ + {A-g'(x)}2/ = 0 (5.4.1) 

which will be used frequently in later sections. Suppose first that A 
is real, q(x) < A, and let q'(x) and q''(x) be continuous. Let 

${x) = J {A— dt, 7)(x) = {A-g'(x)}lt/. 

0 
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Then 

drj 


d^rj 

de 


dr) dx 
dx d^ 




dy 1 q'(x) 


{A-g'(x)} 


dx 4 {A— g(a;)}* 


,dy 1 q'ix) 


*'](a 4 


q(x)}^ 




q''{x) , 5 q'\x) 


4({A-g'(x)}* 4 {X-q{x)}i 

Hence (5.4.1) transforms into 


{A-g(x)}i' 


d^r) 

di^ 


+ V+ 


[i 


q"ix) 


4 {A-<j'(x)}2"^ 16 {A 


-q{x)f\ 


0 . 


(6.4.2) 


This is an equation of the same form as (5.4.1); but in (5.4.2) the 
coefficient of r) in the last term is in general small when A is large, 
or when q{x) is large and negative. 

It follows that, if t; == <^(^), cf) satisfies the integral equation 

<!>{$) = ^(0)cos|+^'(0)sin^— 


f8in(|-T)[l H 

5 q'Ht) ] 


^ 16 {A-g(<)FJ 


(t) dr, (5.4.3) 

where r = ^(^). This can be used as in previous sections to obtain 
asymptotic formulae for <^(|). 

If A is not real, or q{x) > A, | is not real, and the above formulae 
would involve integrals along complex paths. It is not necessary 
to introduce such integrals, since we can obtain the corresponding 
integral equation in terms of the real variable x directly, as follows. 

Let g(0) = 0 (this involves at most a change in the A-origin), 
and write 

P(*) = -g. (6.4.4) 

m ,, 1 < 
-g(*)n ‘ 


Then 


P(x) = {A-g(x)}i^ 




d^ 

dx^ 


_ _ ri f 
_ _ _ 


(il/) 

5 f(x)}^ 16 {A-gr(a:)} 2 j 


+ - 
I ' 1 ^ 
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0 

X 


X 

+ J sin{^(a;)— |(0}{A— gWjty dt 

0 

= / 1+/2 

by (5.4.4) and (5.4.1). Now 

X 

h = [sin{^(^)-^(0}{A-3(0}“*^j + J cos{^(a;)-^(0}^ dt 

= — '»j'(0)A-»sin^(a;)4-[cos{^^(x)-^(<)}i?(<)]o — 

X 

— J sin{^(x)— ^(<)}{A— dt 

0 

= — rj'( 0 )X~^ sin f (a;) +17(0;) — 17(0)008^(0:) — 12. 

Hence 17(0:) = 77(0)cos|(o:)+77'(0)A”^sin^(o:)+/, 

i.e. 77(0:) satisfies the integral equation 

rj{x) = ?7(0)cos^(o;)+'^'(0)A"'^sin^(o:) + 

X 

+ I sin{|^(a;)— dt, (5.4.5) 


(5.4.6) 


where R(t) - ^ ® (5 4 6 

where 16 {A= 5 ( 0 }*- 

In these formulae, if 1(A) > 0, we take 0 < argA < tt, and, if e.g. 
q{t) varies from 0 to oo, arg{A--g(<)} varies from argA to tt. 


5.5. The case q{x) -> oo. We impose the following additional con- 
ditions. Let q{x) tend steadily to infinity, so that q'{x) > 0; let 
q'{x) = 0[{g-(x)}«] (0 < c < f), 

and let q"{x) be ultimately of one sign. Then 


{midx=.o\{j4^dx\ 

J {?{»;)}* LJ J 


= 0[{g(a;)}-J]^' = 0(1), 


(5.5.1) 
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and hence also 

X 


Xo 


dx 




Hence 


Xo 

00 

J \R{x) \ dx 


J {^?(^)}* 


dx — 0(1). 


( 6 . 6 . 2 ) 


is convergent, uniformly with respect to A over any region for which 
I ^“2(^)1 ^ 8 > 0 for 0 < a; < 00 . 


5.6. Theorem 5.6. If q{x) satisfies the conditions of § 6.6, the 
spectrum is discrete. 

Since i^rgA < arg{A— g(^)}i < \tt, 

the imaginary part of ^(x) is positive if 1(A) > 0; as a; ^ oo, while 
A remains bounded, ^ 

!(*) i J {?(<)}* dt, 

0 

and so c~w) _> oo. 

We can now argue with (5.4.5), as we did previously with (6.3.1) 
to obtain (6.3.5). Let ri^{x) — e^^^^^'q(x). Then (5.4.5) gives 

r]i{x) = 7^(0)e^^^®^cos^(a;)+')7'(0)A“le^f^^^sin^(a;)-f- 

X 

-f J c«^(®)-f(«>sin{^(a;)— dt. 

0 

X 

Hence Mx)\ < h(0)| + | V(0)A-‘H- J |i2(«hi(<)l dt, 

0 

so that by Lemma 5.2 

hi(a^)l < {h(0)| + h'(0)A-‘|}expj J|i?(01 dij 

< {h(0)| + h'(0)A-l|}expj J \R(t)\ dt^. 

Hence 

< {l’?(0)| + h'(0)A-M}exp| J|B(«)| 

Applying this to (5.4.5), we conclude as before that, for any fixed A 
in the upper half-plane, as a: -> oo 

7){x) ~ j|^(0)+ii7,'(0)A-i+^‘ J ei«-f>R{t)r){t) 
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Now let ri(x) — {X—q{x)y(f){XyX), Then 


7 ]{ 0 ) = A^sina, 


7;'(0) = — A^cosa — i ^^sina. 


Hence 


|<^(a;,A)| < 0(Ai)exp| J|i?(<)| 


(|A-g'(a;)|i 

for all X and 1(A) > 0; and for a fixed A, as a; -> oo 

Jf(A)e-<f(*) 


^{x,X) 


where 


(A— g(x)}i ’ 
Jf(A) = JAisina-|i{^+^— ) + 


(5.6.1) 


(5.6.2) 


Similarly 


e{x,X) 


00 

+ ^i j eim(t){X-q(t)}mt) dt. 
0 

ifi(A)e-^^ (-^> 

{A— g'(a;)}i ’ 


where i)ifi(A) is obtained from M{X) by replacing sin a, —cos a, (f)(t) 
by cos ex, siiicx, 9(t). Hence 


m{X) = ~M^{X)/M(X). 

Now it clearly follows from (5.4.6) and Lemma 5.2 (as in §5.3) 
that 7](x) — uniformly with respect to A as A approaches 

any interval of the negative real axis. Hence A“iilf(A) is continuous 
up to the negative real axis, and it is also obviously real there, ^(t) 
and R(t) being purely imaginary. 

If we do not choose g(0) 0, the argument shows that 

{A-g(0)}-ii^f(A) 

is real and continuous for A < ^(0). But the whole argument could 
be constructed equally well with an interval {X, oo) instead of (0, oo). 
We should then obtain an alternative expression for M{X) involving 
an integral over (X, oo). Hence {A— g'(X)}~*il[f(A) is real and con- 
tinuous for A < q(X); hence in fact e~^^^M{X) is real and continuous 
along the whole real axis. Since ilf(A) is regular in the upper half- 
plane, it follows from the principle of reflection that il[f(A) is an 
integral function. Similarly il!fi(A) is an integral function. Hence 
m(A) is meromorphic, and the result follows. 

It will be shown later that the result holds if q{x) oo, without 
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any other restrictions. This depends on the theory of the zeros of 
eigenfunctions given at the end of this chapter. 


5.7. Theorem 5.7. Let q{x) < 0, q\x) < 0, q{x) —oo, 

q'{x) = 0{\q(x)\<^} (0 < c < I) (6.7.1) 

and let q'\x) be ultimately of one sign. 

00 

Then if J |^(^)|~* dx (5.7.2) 

is divergent, there is a continuous spectrum over (— 00 , 00 ). 

The conditions imply the convergence of (5.5.2). 

For A real and positive, (5.4.5) now gives 

ri{x) = ri{0)co8^{x)-{-rj'(0)X-^8in^(x)+ 

00 

+ J sin{^(a :)— dt + o(l). 

0 

Hence 

^(x,A) = (A— g'(x)}-*(/x(A)cos|(x)+v(A)sin|(x)+o{l)}, (5.7.3) 

where 

00 

^(A) = A^sina— J sin^(<)i?(^){A— g(^)P<^(^, A) dt, (5.7.4) 

0 

v(A) = J cosmum-qitmit^) dt. 

0 (5.7.5) 

Similarly, if d(x, A) is the solution of (5.4.1) such that 0(0, A) — cos a, 
0'(O,A) — sin a, then 

e{x,\) = {A— g(a;)}-i{/xi(A)cos^(x)+Vi(A)sin|(aj)+o(l)}, (5.7.6) 

where /x^ and are obtained from /x and v by replacing sin a, —cos a, 
and <f> by cos a, sin a, and 0. 

The argument also shows that the integrals in (5.7.4), etc., con- 
verge uniformly with respect to A, and hence that jtx(A), etc., are 
continuous functions of A. 

Again, differentiating (5.4.5), 

r)'{x) = {A— g(a:)pj^A-^7^'(0)cos|(x)— 7^(0)sin^(x)-f 

X 

+ J COs{^{x)—i{t)}R{t)r)(t) c?<j. 
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Applying this in a similar way, we obtain 

^[{X-q(x)y(f>(x,X)] ~ {A-^(x)}»{v(A)cos^(a;)-/i(A)sin^(a;)}, 


whence by (5.7.3) and (5.7.1) 

(f>'{Xy A) ~ {A— g(x)P{i^(A)cos |(a:)— /x(A)sin ^(x)}. (5.7.7) 

Similarly 

d'(x,X) ~ {X—q(x)y{vi(X)cos ^(x)— fi^(X)^in ^{x)}. (5.7.8) 

Hence lim W(<f>, 6) fji(X)v^(X)-ii^(X)p{X), 

x-*-co 

Since W{(f>,d) — 1, it follows that 

fM(X)v,{X)~fi,(X)v(X) = 1. (5.7.9) 


Hence /x(A) and v(A) do not both vanish for any value of A. 

The argument requires modification if A < 0; but then we can 
choose X so that X—q{x) > Oforx ^ X, and we can apply the whole 
argument to the interval (X, oo) instead of to (0, oo). The same 
conclusion then follows. 

Now let A = u-\-iv, where ?; > 0. Taking Xq so that u—q(t) > v 
in a; > we have 


X 

i(x) = J dt = 

0 


X 

^0+ J {u+iv—q(t)}i dt 

Xo 


^ 1^0+ J {«— ffWF di +¥'>’j {u—q{t)}-^ dt +oiv^j |5r(0l~* dt 

Xo Xo Xo 


Hence as a; -> oo 

X 

I(^) (u—q(t)}-i dt, 

and so is large for large x. Hence (5.4.5) gives 


7]{x) = ie-»fw|Tj(0)+»A-*7j'(0)+i J e^(<‘>R{t)y}{t) dt +o(l)|. 

Proceeding as before, we obtain 

7](x) ~ dt^. 

4>{x,X) ~ Jf(A){A— 


Hence 
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where 

Jlf(A) = iA.A.-lip'g^ + 5|2) + 

00 

+ Ji J e^<f>R{t){\-q{t)}mt,X) dt, 

0 

and e{x, A) ~ M^(X){X~q(x)}-~h~^^^^\ 

where Mi{X) is formed from Jf(A) by the same interchange as before. 

According to the general theory of the solutions of (5.4.1), there 
is, for ?; > 0, a solution 

ilf(x,X) = 6{x,X)-{-m(X)(l>{x,X) 

which is L2(0, oo). It follows from the above formulae that 


m(A) = 


]\m 

Jf{A)- 


Also, as V'->0, M(X) -> liJL{u)-\-liv{u) and M^iX) 

since it is easily seen that the integrals involved converge uniformly 

with respect to A. Hence 


lim m{X) = 

v-^Q 


lii(u)+i vi(u) _ + ivi(u)}{n(u) — tV(M)} 


lim{-s-Im(A)} 


fx(u)v^{u) —fXi{u) v(u) 


1 

Ia,^(u)+v^(uY 


Since A;(A) is the integral of this over (0,A), the result stated follows. 


5.8. Theorem 5.8. If the conditions of Theorem 5.7 are satisfied 
except that (5.7.2) is convergent, there is a continuous spectrum in 
(— 00 , 0 ) and a point-spectrum in (0, oo). 


As a: 00 

^{x,X)—i{x, 0) = J dt 

0 


X 

’I 

0 

00 

■/ 


Adi 


{A-q(t)}i+{-q(t)}i 
A dt 

{A— ^(O}* 


X(A) (5.8.1) 


say. Hence I|(a;, A) is bounded, and so cos^(a:) and sin^(a:) are 
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bounded, whether A is real or complex. It is then clear that the 
argument of § 4 holds, in this case, for all values of A, real or complex. 
Thus (5.7.3), (5.7.6), (5.7.7), and (5.7.8) hold for all values of A. In 
particular, all solutions of (5.4.1) are L^{0,ao), and we are in WeyFs 
‘limit-circle’ case, in which the expansion formula involves an 
arbitrary parameter. 

Further /4(A), v(A), /^^(A), and i/i(A) are integral functions of A. For 
consider e.g. (5.7.4). Here {A— g(^)}l is an analytic function of A, 
regular except on the negative real axis; similarly for ^(^, A) and R(t), 
and the integral in (5.7,4) converges uniformly with respect to A in 
any finite region, as the argument clearly shows. Hence /x(A) is an 
analytic function, regular except possibly on the negative real axis. 
However, we could apply the argument leading to (5.7.3) equally 
well with any interval (X, 00 ), where X > 0, instead of the interval 
(0, 00 ). Hence we must also have 

00 

/it(A) - {X~q{X)}mX,X)- j sin^(0i?(<){A-g(<)}¥(^A) dt. 

X 

Hence fi(X) is regular except possibly on the real axis between —00 
and ^(X); and so in fact it is an integral function. Similarly so are 
v(A), jai(A), and i/i(A). 

Now consider the solution 0(XyX)-\-l(X)(f>(XyX) of (5.4.1) which 
satisfies the boundary condition 


{0(6)-fZ0(6)}cosj8-|-{0'(6)-f Z0'(6)}sinj8 — 0 

Sit X — b. This gives 

j/\\ d(b,X)cotp-\-d'(byX) 

^ ^ ” 0(6,A)cot>+f(6,A)‘ 


(5.8.2) 


As cotjS varies, 1{X) describes a circle. According to Weyl’s theory, 
this circle tends to a limit-circle as 6 ^ 00 . Now by (5.7.3) and (5.7.7) 
the denominator in (5.8.2) is of the form 


{A— g(6)}“i{/[4(A)cos ^(b) -f- i/(A)sin ^{b)-\-o(l )}cot 

+ {A- g(6)}*{i/(A)cos ^(6) -/4(A)sin ^^(6 ) + 0 ( 1 )}, 
or, since {A— = {— 

of the form 

{A---g(6)P[{-~g(6)}“^{/4(A)cos^(6)+v(A)sin^(6)-fo(l)}cot^+ 

-fv(A)cos^(6)-/a(A)sin^(6)+o(l)]. 
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Putting {— g'(6)}“^cotj3 = coty, (5.8.2) takes the form 

KX) = /^i( ^)cos{^(^)+y}+vi(A)8in{^(6)+y}+o(l) 

/i"(A)eos{^^(6)+y}+Vi(A)sin{^(6)+y}+o (1) ' 


Putting y = 0), keeping c fixed and making b -^oo, this 

tends to 

m(X) = +c }+ >'i(A) sin {x( A) + c} 

/x(A)cos(x(A)+c}+ v(A)sin{x(A) H-cj ' 


As c varies this describes a circle, which is therefore WeyFs limit- 
circle. 

For A > 0, y(A) is real, and Iw(A) -> 0 as A tends to a real value, 
unless the above denominator vanishes. Hence the part of the 
spectrum on the positive real axis consists of the zeros of the deno- 
minator, which are isolated points since it is an analytic function. 

On the other hand, if A < 0, y ~ Xi+^X 2 > where X 2 ^ 
q{h) = A, (5.8.1) gives 

h 

X 2 = J {?{a:)-A}* dx. 

0 

Hence, using (5.7.9), as A tends to a real negative value, — I{7n(A)} 
tends to 

_____ |8i nh{2x 2(A)} 

I A‘(^)cos{x(A) + c}+ i/(A)sin{x(A)-fc} |2’ 

which is finite and positive. Hence the spectrum is continuous from 
—00 to 0. 


5.9. The zeros of eigenfunctions. This theory depends on the 
following fundamental theorem, due to Sturm. 

Let u be a solution of 

flht 

^+Sr(x)M = 0, (5.9.1) 

and V a solution of ~ (5.9.2) 

aoc 

where g{x) < h{x) throughout the interval {a, b). Then between any two 
consecutive zeros of u there is at least one zero of v. 

Multiplying by v, u respectively, and subtracting, 

u''v—uv*' — {h{x)—g{x)]uv. 


(5.9.3) 
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Let Xi, X 2 be consecutive zeros of u. Integrating from to X 2 j 

-Ml)']®* — J {h{x)—g{x)]uv dx. (5.9.4) 

Xi 

Suppose that v has no zero in (x^,x^. We may suppose without 
loss of generality that > 0, v ^ 0 in {x^.x^)- Then the right-hand 
side is positive (assuming that u and v do not vanish identically). 
Now the left-hand side is 

u\x2)v(x2)—u\Xt)v(x^) 

and 

u\x^) > 0 , u\x2) < 0 , v(xj) > 0 , v{x2) > 0 . 

Hence the left-hand side is < 0, giving a contradiction. 

5 . 10 . We deduce the following theorem. 

Let u be the solution of (5.9.1) such that 

u{a) — sin a, u'{a) — —cos a, 
and V the solution of (5.9.2) such that 

v{a) = sin a, ^;'(a) = —cos a. 

Then if u(x) has m zeros in the interval a < x ^ b, v{x) has at least 
m zeros in the same interval, and the vth zero of v(x) is less than the 
vth zero of u{x). 

Suppose first that sin a: ^ 0, say sin a > 0, so that the left-hand 
end-point is not a zero of either function. 

In view of the previous theorem, it is only necessary to prove that 
v{x) has at least one zero in the interval (a,x^), where x^ is the 
smallest zero of u. Now on integrating (5.9.3) from a to x^, we obtain 

Xi 

u\xf)v{xf) = J {h{x)~g{x)]uv dx. 

a 

If v{x) has no zero in {a,xf^, the right-hand side is positive; but 

u'{xf) < 0 , v{Xi) > 0 . 

Hence we obtain a contradiction. 

The result follows similarly if sin a = 0. 
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5.11. Now consider the eigenfunctions of Chapter I, defined as 
the solutions of 

^+{A-5'(a:)}y = 0, (6.11.1) 

y{a)GOsoL-\-y'(a)mioL == 0, (5.11.2) 

2/(6)cosj8+2/'(ft)sinjS = 0. (5.11.3) 

Consider first the solutions of (5.11.1) and (5.11.2), and suppose 
that sin a 0. Then y(a) ^ 0. The number of zeros of y(x) in (a, b) 
is a non-decreasing function of A, by the above theorem. Let 
\q{x)\ < c in (a, 5). First compare (5.11.1) with 

SH(A+c)2/i = 0. (6.11.4) 


The solution of this satisfying (5.11.2) is 

y^ = cosh{(— A— c)l(a;— a)} — ™^^^^j8inh{(— A— c)l(a;— u)}. 

This has no zeros if A is negative and large enough. Hence the solu- 
tion of (5.11.1) has no zeros in (a, b) if A is negative and large enough. 
Similarly, by considering the solution of 


dx^ 


+ (^-c)2/2 = 0 


with A large and positive, it follows that the number of zeros of 
(5.11.1) in (a,b) tends to infinity with A. 

Also, by the above theorems, a zero of y travels steadily to the 
left as A increases. Hence there is an increasing sequence of numbers 
/x-Q, /xj,... such that y{b) = 0 for A = and y{x) has just m zeros 
in a < a: < 6. 

If sin j8 = 0, the are the eigenvalues. 

Otherwise, we have 





Hence ju—v^ jv) is steadily increasing. 

Suppose that u and v have the same number of zeros in (a, 6). 
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The zero Xj, next before c is a zero of u and not of v, for between a 
and lie at least v, and so exactly v, zeros of v. Hence 


u^{b) 


ju'ib) 

ym 

U(6) 

v(b)i 


> u^x^) 


\u{x^) 






v'(b) 

u(b) v{b) ‘ 


It follows that y\b)jy{b) is a steadily decreasing function of A in each 
interval must decrease from oo to — oo, since 

y{b) ~ 0 at each end, and y'{b) ^ 0. Hence there is just one value 
of A in the interval such that 


y\b) 

y(b) 


— cot^. 


Hence there is an increasing sequence of eigenvalues Aq, A^,... such that 
the eigenfunction associated with A,^^ hOnS just m zeros in a <, x <. b. 


5.12. We now pass to the case in which the interval is (0, oo). 

Let q{x) -> 00 as a; 00 . We shall show that this is quite similar 
to the Sturm-Lioiiville case; there are discrete eigenvalues, and the 
eigenfunction associated with A^, has n zeros. 

We consider the solution y = y{x, A) of the equation 

^^ + {X-q(x}}y = 0 

such that y{0) and y'(0) have given values. Since q{x) -> oo, 
q(x)—X > 0 for > Xx, say. Suppose that for some x^ > Xx 

yiXi,X)>0, y'(xj^,X)>0. 

Then y"{x, A) > 0 (a: ^ xf), y\x, A) increases steadily, and y{x,\)-^co. 

If y{xi, A) > 0, y\x^. A) < 0 

there are two possibilities. If y{x, A) remains positive for x ^ x^, then 
y"(XyX) remains positive, y'(x,X) increases steadily, and so tends to 
a limit (finite or infinite). This limit cannot be negative, or y(x) 
would tend to —oo; if it is positive, y{x) -> oo. If y'(x) 0, then 

y'{x, X) < 0 for a; > x^. Hence y{x) is steadily decreasing, and so 
tends to a limit. Also 

Xi Xz 

I {q{x)—^}y(.x) dx = j y"{x) dx 

Xx Xx 

== < -y'ixj). 

Q 
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Hence {q{x)—)^y{x) is L{x^^ oo), and a fortiori y{x) is L{x^, ad). Hence 
y{x) -> 0 . Hence also y(x) is L\x^,ad), 

On the other hand, if y(x) changes sign, it can do so at most once, 
and then tends to ~oo, by the previous argument. 

Similar arguments hold if y{x^y A) < 0. Summing up, y{x, A) has 
at most one zero for x > Xx, and either y{x, A) -> ±oo, or y{x, A) -> 0, 
y\x,X) -> 0, and y(x,X) is L\x^,ad). 

It follows that each y{x,X) has a finite number of zeros in (0, oo), 
say n(A). By the theorem of § 5.10, n(A) is a non-decreasing function 
of A. The argument of § 5.11 shows that r 2 ,(A) = 0 if A is negative and 
large enough, and that n(A) -> oo as A -> oo. 

Since 7i(A) changes by at least 1 at a discontinuity, it has only 
a finite number of discontinuities in any finite interval. Hence they 
are discrete points. We shall show that they are the eigenvalues of 
the differential equation. 

In addition to the above results, we can show that, if y{x, A) -> 0, 
and /X > A, then y{x, y) has a zero greater than the greatest zero of 
y(x,X), Suppose e.g. that y{x,X) > 0 , y\x,X) < 0 for a: large enough. 
Suppose that y{x, /x) > 0 for a: > a, where a is the greatest zero of 
y{x,X), If y{x,ix) -> 0, (5.9.4) gives 


ao 

(fj,—X)jy{x,\)y{x,n) dx = —y{a,ix)y'(a,X), 

a 

which gives a contradiction as before. Otherwise, if y(x,ix) > 0 for 
large x, then y'{x, /x) > 0 for large x. Hence 


for large x, and we obtain a contradiction again. 

In particular it follows that to each interval where n{X) = constant 
corresponds at most one value of A for which y[x, X) is and such 
a value must be the right-hand end-point of the interval. 

Now consider a decreasing sequence A^, A 2 ,... tending to a limit A. 
Since ^(A^) — 7i(A+0) for jp large enough, we can suppose that 

72,(Aj) ™ ^(A2) — ... — k. 

Let the zeros of y{x, X^fj be Then it follows from the above 

argument about the interlacing of zeros that for every m ^ k 

^(1) <;^ <C' cT 
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Let c be a constant such that q(x)~-\ ^Ofovx^c, Then 

remain less than c for all p, since 2 /(^, A^) has at most one zero 

greater than c. Hence 

lim ^ c 

p—^co 

exists for m = 1, 2,..., k—ly and the are zeros of y(x,X). They 
are all different, since if two coincided we should have y((i^, A) = 0, 
A) = 0, which is impossible. Conversely, any zero of y(x,X) is 
a limit-point of the Hence y{x,X) has k or k~l zeros according 
to whether tends to a finite limit or to infinity. 

Now clearly y\ofjf^,Xp) has the same sign for all p, say +• If 
-> 00 , then > c for sufficiently large p. Since none of the 
y{x, Ap) are we must have 

y'{x,Xp)>0, y{x,\p)<Q (c < a: < ojf)). 

Making p -> oo, it follows that 

2 /'(a;,A)> 0 , ?/(a;,A) < 0 {x > c). 

Hence y{Xy X) is L^, 

On the other hand, if tends to a finite limit then 
ij'{x,X) > 0, y(:x,X) >0 {x> aj. 

Hence y(x, X) is not 

Summing up, either 7i(A) = ri(A-f-O) or /^(A) = 7^(A+0)— 1. In the 
latter case, y{x, X) is L^, in the former case not. 

If we argue similarly with an increasing sequence X^, A 2 ,..., the 
are decreasing, and therefore all tend to finite limits. Hence 
72,(A) = n(A— 0). 

It follows that to every n corresponds an interval of values of A 
in which /i(A) — n, and that the right-hand end-point of each such 
interval is such that y{x, X) is L^, 

Now let A' and A" be interior points of an interval where 71(A) is 
constant, and let Xq be greater than the greatest zero of y(x, X) or 
y'(x,X) for A' A ^ A". Then y(x,X) is positive increasing (or 
negative decreasing) for x'^ Xq, A' ^ A ^ A". Let 7n be the lower 
bound of \y{x^, X) | for A' < A < A", so that m > 0. Then 

|y(a;,A)| ^ m {x x^, A' ^ A ^ A"|. 
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We can now take y{Xy A) to be the function (f){x, A) of Chapter III. 
By (3.4.1) ^ 

x(x,A")— x(a:,A') == J <f>{x,u) dk(u) 

A' 

A' 

^ m J dk(u) 

A' 

> m{i;(A''-0)-A;(A'-l-0)}. 

But is for every A. Hence 

A;(A'+0) k(X"--0). 

Hence &(A) is constant in the open intervals where 7^(A) is constant. 
Let the points of discontinuity be A^, the saltus at A^, k^^. Then by 
(3.4.1) 

x(x,A)= 2 K4>i^>K)- 
k<A»<A 


By (3.4.3) 
and by (3.6.3) 


2 Kf <!>(!/> K)f(p)<ip 

K<An<A g 


00 

/(^) = j ^y,X„)f{y) dy. 


Hence the expansion is a series. 

Also, comparing (3.5.1) with (3.4.4)-(3.4.6), we have 


p(m2)-p(— %) = J w) dg-Su). 

Ui 


Hence in the present case p{x) is constant except at the i)oints A^. 
Hence by (3.6.2) 0(0:, A) is regular except for simple poles at the 
points Ay^. 

It follows also that m(A) is meromorphic; for if/(:r) == 0 for x ^ X, 
(2.3.2) gives 

X X 

<l)(a;,A) = 0{x,\) j <l>{y,X)fiy) dy +^(x,A) J d{y,X)f{y) dy + 

0 X 

X 

+m(A)^(a:, A) | <f>{y,X)f{y) dy, 

0 


and the result follows. 
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5.13. The interval (—00,00). Theorems on the nature of the 
spectrum when the interval is — 00 < a; < 00 can be deduced from 
those already obtained for the interval ( 0 , 00). 

Let g(a:)->oo as —00 and as a:->oo. Then both ^^(A) and 
m2(A) are meromorphic. Hence the functions 

1 mi(A)mg(A) fsioix 

mi(A)— m2(A)’ mi(A)— m2(A)’ mi(A)— m2(A) 
occurring in the (—00, 00) formula are also meromorphic, so that the 
corresponding expansion is a series and the spectrum is discrete. 
The argument of § 5.12 also extends without difficulty to this case. 

Next let q{x)">co as x->oo, and q{x) ~oo as a; ->—00, 
/ \^{^)\~^ being divergent. Then I{mj(A)} tends to a finite non- 

—00 

zero limit along the whole real axis, while I{m2(A)} tends to zero in 
general but to infinity at certain discrete points. Hence the imaginary 
parts of the functions ( 5 . 13 . 1 ) tend to finite limits, which can vanish 
at most at discrete points. The spectrum therefore extends con- 
tinuously from —00 to 00. In fact it is clear that, if q{x) -> 00 at one 
end of the interval (—00,00), the result is the same as if we had the 
other half-interval only with a given boundary condition at a; = 0 . 
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VI 

A SPECIAL CONVERGENCE THEOREM 


6.1. Ik Chapter III it was proved that the general expansion formula 
represents an ‘arbitrary’ function f{x), provided that f(x) is twice 
differentiable, that /(a:) Sbnd f"{x)—q{x)f{x) are ^^(O, oo), and that 
the ‘boundary condition at infinity’ (2.7.2) is satisfied. In this and 
later chapters we shall show that, by imposing further conditions on 
q(x), we can relax those bearing onf(x); and in fact it will be shown 
that in certain cases the expansion converges to f{x) if f(x) satisfies 
conditions similar to those for the convergence of an ordinary Fourier 
series. 

In this chapter it will be assumed that the integral 

00 

^\q{x)\dx (6.1.1) 

0 

is convergent. This is perhaps not the most interesting case, but it 
is the simplest, and the analysis suggests the method which will be 
used later in considering cases in which q{x) -> oo. The advantage 
of this and other special assumptions is that they enable us to prove 
asymptotic formulae for the functions (l>(x,X) and 0(a;, A). In the 
general case these formulae are not known, and we have to rely on 
a direct discussion of the function <I)(x, A). 

6.2. What we require primarily is a solution of the equation 

g + {A-g(x)}y = 0 (6.2.1) 

which is small when 1(A) is large and positive. Consider the integral 
equation 

X 00 

x{x) = + ^ j * dy J e»»(*'-*>^(y)x(?/) 

» ® ( 6 . 2 . 2 ) 

in which A = On differentiating this twice, it is at once verified 
(formally) that y ^ xi^) satisfies (6.2.1). Now a solution of (6.2.2) can 
be obtained as follows. 

Let 


(6.2.3) 
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and, for ^ 1 , a- 

J Ay + 


+ 


00 

^ J ei^y-^'^q{y)xJ.y) dy. (6.2.4) 


Then 


/iisoc If* f* I 

X 2 i^)-Xii^) = J ^^y'> + J dy\. 


Let s = a-\-it {t ^ 0 ), and let J denote the value of the integral 
( 6 . 1 . 1 ). Then 

p — tx T 

IX 2 (^)-Xi(^)l < {6-2.5) 

Next 3 . 

X3(*)— X2(^) == 2 ^^ J c^^*""^«(2/){x2(2/)-Xi(2/)}% + 

0 

00 

+ ^- J e**<"-*>^( 2 /){x 2 (?/)-Xi( 2 /)} dy, 

whence by (6.2.3) 

( X CO V 

J \y{y)\ dy + J e‘^^^-y'>\q(y) \ dy\ 

Q X 



and so on generally. It follows that, if |5| > ^Jy the series 


( 6 . 2 . 6 ) 


2 {x«-n(*)-x»(^)} 

n=l 

is convergent, i.e. that Xni^) f^^ids to a limit, say x(^)- 
For every n 


lxn(^)l < Ixi(*)l + lx2(^)— Xi(*)l + -" + lx«(*)— X«-l(^)l 



g-te 

r=|7//H’ 


(6.2.7) 


and xi^) satisfies the same inequality. It follows by ‘dominated 
convergence’ that we can make n->co under the integral sign in 
(6.2.4). Hence x(^) satisfies ( 6 . 2 . 2 ), and so also ( 6 . 2 . 1 ). 
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For a fixed Sy or 5 in a bounded part of the region t 0, l^l > 
(6.2.2) now gives 


^{x) = 


sx j* 

J- I e-i^Vq(y)^(y) dy - 
0 

00 00 

pisx (* pisx (* 

— J e-^^«q{y)x{y)dy +-^ j e-^^^y-^\{y)x{y) dy. 


The last two integrals tend to 0 as a: -> 00 ; hence 

Xix) = ei^^{K{X)+o(l)}, 

00 

where A"(A) = 1 + L J e-^^«q{y)x(y) dy. 


( 6 . 2 . 8 ) 


(6.2.9) 


6.3. Let i/j{XyX) = d{XyX)-{-m(\)(l>(XyX) denote as usual the solution 
belonging to L^(0y oo). Then 

0(0;, A) == KS)xM+K,(X)cf>(XyX), 

Now L^(^yCX))y but 0(o;,A) is not, by (5.3.5), at any rate 

if |A| is large enough ; for (5.3.6) gives 

If (A) ^ I sin a (sin a 9^0)y 


Hence iL 2 (A) = 0, i.e. 


cos a 
2i8 


(sin a — 0). 


i,{x,X) = K,(X)x(x,X). 


As in §5.3, the asymptotic formulae (5.3.5) and (6.2.8) can be 
differentiated, i.e. we can apply similar arguments to the differen- 
tiated functions, and obtain (as x->oo) 


0'(a;) ~ —isM(X)e~^^^y 
x(^) ^ isK(X)e^^^, 

Hence 

Pf(0,x) ~ M{X)isK{X)+isM(X)K(X) = 2i8K(X)M(X), 
Hence If (0, 0) ~ 2i8K(X)M(X)K^(X). 


But 1F(0, 0) — 1. Hence 

2i8K(X)M(X)K^(X) = 1. 

x(a;,A) 


Hence 


^x,X) 


2isK{X)M{X) 
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6.4. Theorem 6.4. Let q{x) be L(0, oo), andf(x) oo) and also 
L(0, co). Then the eigenfunction expansicm is 


0 0 

00 

+ 2 j <f>iy,K)fiy) 

0 

where /ut(A) and v(A) are defined in § 5.3, and, in the sum, r.^^ {n = 1, 2,...) 
runs through the residues of m(A) at its poles on the negative real axis, 
and irr^ is the saltus of k{X) at A — 0. The expansion converges under 
the same conditions as an ordinary Fourier series) e.g. if f is of bounded 
variation in the neighbourhood of x, its value is 

i{/(x+0)+/(a;-0)). 

Let 

X 00 

•JCa;, A) = <A(a:,A) J </>(y,X)f{y)dy +(l>{x,X) j ili{y,X)fiy) dy 
0 0 


a h X \ / it*+S 00 V 

+ J j+.^(a;,A)| J + J j 

x—h' ' X ic+S' 


R+U 

Consider J Q>(x,X) dX, 

— R-\-i€ 


This has the same value whether it is taken along the straight line, 
or round the semicircle above it. Suppose, for example, that sin a 
in the above formulae. Then for \s\ ^ J 


Also 




2e~^^ 


2|«| I Jf(A)ir(A)| 
<^(x,A) = 0(e‘% 


< 


Ae-*v 


Hence (I>4(x,A) = ofe J e-ty\f(y)\ dy\ = 

^ a;+5 ^ 

The integral of this round the semicircle tends to 0 as i? oo, for 
any positive S. A similar argument clearly applies to 
Now consider 4)3. For x fixed, or in a finite interval, 


(xi 


(J 


Ae~^^ 
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Also K(X) = l + 0 (\s\-^), M(X) = isina+0(|5|-i). 

Hence + 

^s^Ina( \l«|/j 

Also — cossa:sina|l + oLj 

by ( 1 . 7 . 3 ). Hence 

x+S f aj+S 

fii8x\p~i8x C I r 

<1)3(0:, A) = - J e-<'/(y)# + 0 1^3 J e-‘yf{y) dy 


The last term is 


0 


x+8 

ixj '^<*'>1'*^ 


and the integral of this round the semicircle is 

/ic+S 

0 


J \f(y)\ 


dy 


which can be made as small as we please, by choice of S. The term 
involving also gives a zero limit. The other term is the same as 
in the case of an ordinary Fourier series; and similarly for Og. 
Altogether it follows that, e.g. in the bounded variation case, 
it+u 

lim f A) dX -> U7T{f{x~\-0)-{-f(x—0)}. 

R—^^co J 

~R+i€ 

This is true uniformly for 0 < e < 1. 

Since f(x) is L^( 0 , 00), the analysis of § 3.4 aj^plies, so that 

R+i€ X R 


lim I 


where 


- f <I>(o;,A)rfA U=i f <f>{x,\)dgi{X), 

7T J 77 J 

-R+U ^ -R 


9 iQ) = J x( 2 /> ^)fiy) ^y^ x( 2 /> ^) = / “) 


By ( 6 . 3 . 10 ) 


A 

xlS'.A) = J 


4 >kih 




du, 


taking /c > 0. Since (j){y, u) is bounded for 0 < ac ^ ^ A, 7/ -> 00, 

it follows that 

A CO 

du 


9 'i(A) = J 




J Hy , '“)/(«/) dy- 



6.4 A SPECIAL CONVERGENCE THEOREM 123 

Hence, if p > 0 , 

R R CO 

J A) = J / Hv. my) dy. 

p p 0 

In the interval — iZ < A < — p' (p' > 0 ), 0 r^(A) is constant except 
for a finite number of discontinuities at the poles of m(A). Hence 

—p' 00 

f (l>{x,X) dg^{X) = TT 2 'rn’i>{^^K) f Hy>K)J{y) dy, 

-R -K<A»<-P' J 

each (^(y, A^J being L^(0,co), as in §2.5. Finally 

p, 

lim I A) ~ 0), 

p,p'— >0 

-p 

where IZq is the saltus of (/^(A) at A = 0 . It easily follows from the 
analysis of § 3.4 that 

00 

■^0 = J 4>(y> ^)f{y) dy, 

0 

where is defined above. This completes the proof. 

The form of the result makes it seem possible that the condition 
that f{x) is may be unnecessary; it has been introduced here so 
that previous theorems which have been proved with this condition 
may be used. 


Stone (1). 


REFERENCE 
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7.1. In this chapter we suppose that q{x) -> oo, so that the differential 
equation 


g + {A-5(*))!/ 


(7.1.1) 


has discrete eigenvalues Aq, A^, — We then ask how the distribution 
of these eigenvalues is determined by the function q{x). 

If A ~ A„, then y = has just n real zeros (§5.12). On the 
other hand, by comparing the distribution of the zeros of (7.1.1) 
with that of the zeros of known functions, we can obtain approximate 
formulae for the number of zeros in terms of A,^. Approximate rela- 
tions between n and A,^ are thus obtained. 

To carry out the analysis considerable restrictions have to be 
imposed on the function q{x). We begin by assuming that q\x) and 
q'\x) exist, and that q'{x) > 0, q"{x) > 0 for a; > 0. It is also con- 
venient to assume that ^(0) = 0. This is no restriction, since it 
merely involves a choice of the A-origin. 

Let p = ^?(A) be defined by the equation q{p) — A, and let 

X 

^ ^ ^ {X-q{t)f dt (0<a:<p). (7.1.2) 

0 


Let denote (temporarily) a twice-differentiable function of 
and let 


Y = Y{x) = {X-q{x)}-m^). (7.1.3) 


Then 


Y' I q'{x) 

Y - -^4X-q{xy 


I! 

7 


Y^ 1 m r + m ^ ^4 X-q{x)'^4 {A-«(a:)}2- 


Hence 
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In particular, if <^(^) = cos^, then and Y satisfies 

the differential equation 


dW 


jA-g'(a;) — i 


X—q(x) 


A \y 

16 {A-3(x)}2/ 


0. 


(7.1.5) 


7 . 2 . Theorem 7.2. If q{x) satisfies the above conditions , 

Pn 

^ ^ J l^n-q(^)}* dx -I, (7.2.1) 

0 

tvhere = p(X^). 

For 0 <i X < p, the coefficient of ^ in (7.1.1) is greater than that 
of Y in (7.1.5). Then by the theorem of §5.9, between any two zeros 
of Y there is at least one zero of y. Hence if Y has I zeros in (0,p^), 
then n ^ 1. 

The zeros of F, i.e. of cos^, are at the points ^ = (w+i)?^, for 
values of m for which this lies between 0 and 

Pn 

j R-3'(<)F dt. 

0 

Pn 

Hence < J {A^— ^(0}^ dt < 

0 

and the theorem follows. 

Since q{x) is convex downwards, q{\Pq^ ^ lA/^j and 

^ “ J dx—l 

0 

jLTT 

Hence A„ < STr\n-{-lflpl. (7.2.2) 

In particular, A^ ~ o{7i^). 

Actually we can prove this if q(x) is any function which is twice 
differentiable and tends to infinity; for compare the solution of 
(7.1.1) with 

Y^ ^ {A— 1— ^(a;)}”icosj^ J {\~l—q{t)}^ 

The coefficient of y in (7.1.1) exceeds that of F^, in the differential 
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equation for 1^, over any jfinite interval (0,X), if X is fixed and A is 
large enough. Henpe 

” > " J — f 

0 

>AXXi,. 

Since X may be arbitrarily large, the result follows. 

In the Sturm-Liouville (finite interval) case we have of course 

K ~ 

In the case in which the interval is (— oo, oo), and q{x) is an even 
function, and g'(0) ~ 0, there is a slightly more precise result, which 
will be of use later. In this case, if n is even, is an even func- 
tion, with positive zeros and \n negative zeros. Let 

J {X^-q{t)}i dt\. 

*- 0 -* 

Then 7(0) ~ == 0 = ^^(0). Hence, by the theorem of 

§5.10, 0„(a:) has at least as many zeros in (0,p,J as Y{x) has. If 
7 (x) has I zeros in this interval, then 

Pn 

/ {A, -?(<)}* < {l+\> < (ln-\~l)7T. 

0 

Pn 

Hence n ^ - J “!• (7.2.3) 

0 

If n is odd, has \n—\ positive zeros. Let 

•-0 

Then 7(0) = 0 = ^^'^(O), 7'(0) = ^^^(0). Hence, if Y{x) has I zeros 
for 0 < a; < p^, 

Pn 

j {K-qm^dt < (l+l)n < (i«+i)77, 

0 

and (7.2.3) follows again, 

7.3. In the problem of the upper bound of n, we have to assume 
rather more about q{x). 

Theorem 7.3. Let q{x) satisfy the above conditions, and also 

q'\x) = o[{ 2 '(;c)P]. (7.3.1) 
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Pn 


We observe that, if q'(x) is non-decreasing and 
q"(x) = 0[{q’{x)}y\ (1< y < 2), 


(7.3.2) 


(7.3.3) 


X 

q'{x) = J q"{t) dt ^-^'(O) 

0 

= o[j{g'(<)}vd<] 

= 0j^{g'(a:)}y-ij5'(<)rf«j 

= 0[{q'{x))y-^q(x)l 

{q'{x)Y-y = 0{q{x)}, q'{x) = 0{g'(x)}y<2-r); (7.3.4) 

and similarly with o instead of O. Thus (7.3.1) gives 

q’{x) = o[{g(x)}8]. (7.3.5) 

All such conditions imply a certain restriction on the regularity of 
the growth of q{x), but not on its rate of growth. 

Now let be the function obtained by replacing A by (ijl > 0) 
in y, with = cos^. Then the coefficient of Y^, in the differential 
equation for Y^ corresponding to (7.1.5), exceeds \—q(x) if 

5q'^x) , q"(x) ^ , .. 


5'(a;)}^^4{A+/i— g'(a:)} ^ 


(7.3.6) 


This is true for 0 ^ a: if 


59'HPn) I q''{Pn) ^ „ 

Let /X = Using (7.3.1), this is seen to be true if is large 

enough. 

If the number of zeros of Y^ in the interval (0,p,J is m, then 
m ^ n~l, and 

Pn 

(m— 1)77 < I {A„+jw.— g'(x)}i dx < (m+^)TT. 


Hence 


Pn 

“ J +1* 
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To prove (7.3.2) we now require 

Vn r -I 

/ = o J dx\. 

n . L « J 


Now 




(7.3.7) 

< 


for q{x) ^ Also, since q{x) is convex downwards, q{\Pn) 


Hence 


Pn iPn 

J dx^ j (|-A„)i dx = IPrXlKf- 


0 0 
Hence (7.3.7) is true if p, “ o(A,J, which is true by (7.3.4). This 
proves the theorem. 

Suppose, for example, that q{x) ~ x^ {k > 0 ). Then^^ = A^/^, and 

J dx= j (A„-x*)i dx 

0 0 

1 

= iAji+V*J di 

_ nmm 

kv\i+iik) ” • 

Hence it follows from the theorem that 

tatr(5+i/t) i«‘«) 

~ rf(i)r(W”) ■ 

This is easily verified in particular cases; e.g. if A; = 1 , the formulae 
are those of § 4.12. As A -> oo 

~ {^:^^H^Os(l^i-ln-ln)+COs(l\i+^7r-hr)} 

Hence Pi — ^77 ~ (w+|) 7 r, 

agreeing with (7.3.8). 


7.4. The following well-known property of Bessel functions is now 
required. 

Lemma 7.4. In the interval (0, A), J^{x) has A/7r+0(l) zeros. 
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See Watson, Theory of Bessel Functions^ § 16.4; alternatively, consider 

/(r.) =, 

It is clear from this formula that f{x) has at least one zero in each 
interval 

rmT-\-\v7T — Jtt — 8 ^ a; ^ — Jtt-I-S 

if 8 is a fixed positive number less than and m is large enough; 
for fix') has opposite signs at the ends of such an interval. Also 

Hence f'(x) is of constant sign throughout each such interval, so 
that /(a:) has at most one zero in the interval. The result now follows. 


7.5. Application of Bessel functions. In the problem of the 
upper bound of n, the comparison function used above is not the 
best that can be found. It was observed by Langer that a better 
approximation to a solution of (7.1.1) is obtained by using Bessel 
functions. 

Let us now take 


By (1.11.2) 

Hence (7.1.4) gives 


no _ _i 

n) Wp)-i(or ■ 


(7.5.1) 


dW 



i){A-g(x)} 


4 X~q{x) 


5 q'HO ' 

16 {A— g(x)P_ 


F = 0. 

(7.5.2) 


If V =z ±1, we return to the previous formulae. If v = J, the infinity 
of the new term a>t x = p just cancels those of the following terms. 
It is clear that this must be so; for as ^ ^(p) 


and p p 

Op)-^ = J {A-g'(<)}‘ dt ~ ^ J {A-gW}* dt 

X X 

= l{A— g(a:)}*/?'(a;)- 

Hence F ~ .d{ 5 ''(x)}-'’{A— g(x)} 


and no negative powers of A— g(a:) arise on differentiating. 

4967 • a 
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For our purpose, however, it is better to take v == 0, so as to have 
something to spare in comparing (7.1.1) with (7.5.2). 

Let iV(ic)4- 1 be the number of eigenvalues not exceeding x, so that 

^(A)-^ (A,,_x<A<AJ. 

Let X = 'p{y) be the inverse function of ?/ = q(x), so that ^(A) = p, 
PiK) = Pn- Then 

V 

N{X) > i J dx +0(1). (7.6.3) 

0 

For this is true for each A^, by Theorem 7.2; and so, if A^_i < A ^ A^ 
(denoting the integral by /(A)), 

iV(A) > iV(AJ-l 
>i/(A„)+0(l) 

7T 

> i/(A)+0(l), 

TT 

since /(A) is steadily increasing. 

We shall now prove 


Theorem 7.5. Lei q\x) -> oo, q'\x) > 0, and 
q'\x) < {q’(x)}y (K y < I) 

for X > Xq. Then 


N{X) = i J {A-g(a;)}i dx +0(1). 


(7.6.4) 

(7.6.6) 


Let 


..(.r \) - ^ ^ ^7 5 6) 


4{^(i?)-|(^)P 4 X-q(x) 16 {A--^(:r)}2 

so that (7.5.2), with v = 0, is 
dW 


dx"^ 


+ {A-?(x)+x(a;,A)}r = 0. 


(7.5.7) 


Now 


p 

X 

P 

< ^ I dt 


%{X—q(x)yiq'(x). 


(7.6.8) 



7.5 


THE DISTRIBUTION OP THE EIGENVALUES 


131 


Hence 




16 (A— g(a!;)p 4 X—q{x) 

^1 q'^x) 1 

4 {A— 4 X—q(x) ’ 

Hence A) > 0 provided that 

X-q(x) < q'\x)lq"{x), 
or, by (7.5.4), provided that 


5 

16 


(7.5.9) 


(7.5.10) 

(7.5.11) 


If (7.5.11) holds for ^ x ^ p, it follows from Lemma 7.4, with 
= 0, that the number of zeros of y in the interval x^^x does 
not exceed 

V 

i J {A-^(a;)}i ix +0(1). (7.5.12) 

Xx 

It remains to consider the interval 0 < a: < and here the most 
natural source of information seeiias to be (5.4.5). Suppose, for 
example, that rj'(0) = 0. Then, writing rj-^(x) = r]{x)/rj{0)y 


X 

rjlix) = cos|(a:)+ J sin{f(a:)— ^(«)}i?(0T?i(0 dt. 
0 


Suppose that 
Then 


by Lemma 5.2, and 


Xi 

I |i?(<)| dt < S. 

0 

hi(^)l <i+Jl^(<hi{<)|d<. 

0 

|7ji{x)| < ®^p|/ i-®(OI d^. 


(7.5.13) 


X ^ f 

J 8in{^(x)-^(<)}i?(<)’?i(0 < / \R{t)\eV^^''^'^' dt 

0 0 

= expj J \B(t)\ diS-l 

' 0 ' 

< e^~l < x^-g < 28 > 


ifS< + 
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Consider the interval nir ^ (n-{~l)TT, and suppose, for ex- 

ample, that n is even. If S < 1/(2 V2), it follows that r]i(x) is positive 
for niT ^ ^ ^ negative for (n+f)7r ^ ^ ( 72 . 4 - 1 ) 77 . 

Hence there is at least one zero in the interval (ti+Dtt < ^ < (7i-|-f )7r. 
Also 

X 

Tj[{x) = — f{a:)sm^{a;)+f(a;) J cos{^^(a:)— |(0}iZ(<)T/i(<) 

0 


and a similar argument shows that r^'i(x) is negative throughout the 
interval ( 72 . -j- 4)77 < ^ < (72.+ Hence r]^{x) decreases steadily, and 
so has only one zero in the interval. The number of zeros of 7]{x)y 
and so of 7 /, in (0, x^) is therefore 

Xx 

i J {A—g(x)}i dx +0(1). 

0 

The theorem then follows from this, (7.5.12), and (7.5.3). 

It remains to prove (7.5.13) with 8 < 1/(2V2) and A sufficiently 
large. Now 


f ^'^(0 ..2 q\x,) 2^'(0) 2f q\t) 

J {A-?(0F 3{A-g(a;,)}» 3 A» 3 J {A-?(0}* ' 

0 0 

The first integrated term tends to zero, since 

X~q(x) > {q\x)Y~y (x < x^). 

It is therefore sufficient to prove that 


q"{t) 

{\-q{t)Y 


dt — >■ 0. 


Let qix^) ~ |A. Then this integral is less than 


Xq Xx 

J (iA)» ^ J 


dt 


Xq 

g^K)-g'(O) 

(|A)» 


1 r 1 1 ] 

-3y/2[{3'(a:„)}2-3r/2 ' 


This all tends to zero (the first term since q'ix^) ^ q'(v) ~ o (Al) by 
(7.3.5)). The theorem therefore follows. 
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7 . 6 . The following consequences of Theorem 7.5 will be required 
later. 


Theorem 7.6. Under the conditions of Theorem 7.5, 

iV(A+VA)-iV(A) = 0{i>{A)}, (7.6.1) 

N{X+ix)~N{X) ^ 0{fiX-ip(X)} (VA < /X < A), (7.6.2) 

N(X+fji) = N(X){l + 0{p/X)} (VA < ^ < A). (7.6.3) 

By Theorem 7.5 

p(lA) 

^{N{X+^JX)-N{X)} j [{A+VA-g(a;)}i-{A-ff(x)}»](ix + 

0 

p(X) 

+ / [{'^+VA— ^(x)}*— (A— g(x)P]c?x + 

D(iA) 

2)(A4*VA) 

+ r {A+VA— g'(x)}idx +0(1) 

PiX) 

= /i+/2+/3-f-0(l) 

say. Now 

{A+VA-5W)<-{A-5(»))t = jX +VA-;Wp+(A-i(» -)}i- 


For X < i?(|A), q{x) < |A, so that this expression is bounded. Hence 

= 0{i.(P)} = 0{j>(A)}. 

PiX) 


Also 


I 2-0 


PiX) 

S(lA) 


VA dx 


= 0 


VA 


Now 

Hence 


(A+VA— g'(x)}‘ 

X 

q[x) == J q\t) dt < xq'(x). 


[ 


q'{x) dx 


J {A+VA-^(a:)P 

p(iA) 


b(iA)}]- 


<imm^ iXP) -pm' 


Hence 4 = 0{j)(iA)} = 0{i?(A)}. 

Finally 

h = 0[Xi{p{X+^X)-p(X)}] = 0{Aip'(A)} 



= 0{p{X)}, 
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and (7.6.1) follows. (7.6.2) follows on applying (7.6.1) 0(pA-*) times, 
and observing that, since p{x) is concave downwards, 

p{x) < 2p{\x). 

Also (7.6.3) follows from (7.6.2), since by (7.6.3) 

N{\) > AAip(|A) > AA»p(A). 

REFERENCES 

Milne (1), Langer (1), (2), (3), Titchmarsh (8), (10). 



VIII 

FURTHER APPROXIMATIONS TO N{X) 

8 . 1 . The following argument is used by physicists. In quantum 
mechanics the equation 

g + = 0 (8.1.1) 


is considered, where V — V{x)y and E, m, and h are constants, h 
being small. Let be defined by 



f = expj?JJ 

Xdx^. 

(8.1.2) 

Then 

# 27rix , 

dx h 

(8.1.3) 


d^tfj 27ri(/r dx , 

dx^ h dx 



Substituting in 

(8.1.1), and dividing by 0, we obtain 





or (8.1.4) 

If h were zero, it would follow that 

X = {2m(J?-F))i = x„ (8.1.5) 

say. Now suppose that y is expanded in a series of the form 

= ■■ <*•’•** 

Substituting in (8.1.3), and equating to zero coefficients of different 
powers of y, we obtain equations which determine the Xn succes- 
sion. They are (8.1.5), and 

= - i Xn-».Xm {n = 1, 2....). (8.1.7) 

U'X m=0 


Xi = — 


iXo 

2 Xo 


1 v 

4E-V’ 


Thus 
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V = ^ I , F- I 

2xo 2xo \l^E~ vr HE-Vr ^4(E~-V)i 
_ 5F'2 V" 

and so on. 

Now apply these formulae to the problem of the relation between 
n and A^. Replace a; by a complex variable z, and suppose that V{z) 
is an analytic function of z, regular in any relevant region. Then 
any solution </r = iff{z) of (8.1.1) is an analytic function of z. Suppose 
that there is a discrete spectrum, and let E — E^ be the (rt4-l)th 
eigenvalue. Then iff ~ ^ zeros. Suppose that there is 

a contour C which includes these zeros but no others. Then 


by (8.1.3). 


n — 
Hence 


j_ rm 

2TTi] ijtjz) 



dz 


o c c 

Let V{z) = E at z — a and z b, and suppose that these are 
simple zeros of V{z)-~E, Then Xi(^) simple poles at these points 
with residue — i at each. Hence 


c 

Also, integrating by parts. 


Hence 


/ 


n^.dz 


{E-V{z)}i 

J Xzi^) dz = 


1 


32(2m)* 


J r Z: 

U (K- 

c 

r JL' 

J {E- 


^(z)_ 

V{z)f 


dz. 


v(z)y 


dz. 


Hence we obtain the formula 





dz- 


h 


128772(2m)* 


J 


F'2(2;) 

{E-V{z)Y 


dz +... . 


(8.1.8) 


All the above analysis is, of course, purely formal. The validity of 
an expansion of the form (8.1.6) remains to be investigated; and h 
is not arbitrarily small, but merely a numerical constant. Even the 
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formal appearance of successive approximations with h small, upon 
which (8.1.6) depends, vanishes when we come to (8.1.8); for if E is 
to be an eigenvalue of (8.1.1), it must depend on h. Suppose, for 
example, that V(x) — ax^. On making the substitutions 



(8.1.1) takes the standard form 


dX^ 


f (A-Z^).^ 0. 


(8.1.9) 


If Xq is the smallest eigenvalue of this, that of (8.1.1) is 


E,= 


Xo/ h^a U 
4 \7T^my 


It is easily verified, e.g. by making the substitution z = that 
each term on the right-hand side of (8.1.8) is independent of Ji. 

The above analysis is known as the B.W.K. (or W.B.K. or W.K.B.) 
method. 


8.2. The equation (8.1.1) reduces to the standard form on putting 

STT^mE . SirhUj.. . , . /o o 

We have therefore to consider the equation 

g + {A-^(2)}^=0 (8.2.2) 

in which the usual real variable x is replaced by a complex variable 
z, and q(z) is an analytic function of z, regular in a certain domain 
containing 2 — 0. The argument of § 1.5 is easily adapted to show 
that, if iIj(z) satisfies given boundary conditions at 2 ; “ 0, it is an 
analytic function of z, for each A, regular in the same domain as q(z). 
Since the solutions to be considered are eigenfunctions, A may be 
supposed real. 

The situation seems to depend a good deal on the nature of the 
function q(z). In order to obtain precise results, we consider the 
special case q(z) = z^y where k is an even positive integer; but it will 
be seen that it is possible to extend the method to any analytic q(z) 
which mimics this special case sufficiently closely. 

4957 m 
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8.3. It is a question of obtaining approximations to solutions of 
(8.2.2) in certain regions of the 2 :-plane. 

Let A be real and positive, and let p denote the real positive 
value of A^/*. Let 

^(z) = J (A— m;*)* dw, 17 ( 2 ) = (A— z*)ii/r„(z), (8.3.1) 

0 

where (X—w^)^ and {X—z^)^ reduce to the real positive values of A* 
and A^ as 0, 2 ; -> 0. As w passes above p from a real value less 
than p to a real value greater than p, arg(A— decreases from 0 to 
—77. Hence (A— becomes --i(w^~X)^, where the square root is 
real and positive. Hence becomes exponentially large as 2 ; -> 00 
along the real axis. 

The analysis of §5.4 obviously extends to complex variables, so 
that r](z) satisfies the integral equation 

r)(z) =. 7y(0)cos|(2:) + 77'(0)A-*sin^(25) + 


where 


R(w) — 


z 

+ J sin{|(z)— dw. 


0 

4 {A— 16 {A— g'(w)p 

k(k — 6 
4(A— ~I6 (A-m;*)* 


kw^~’^{4c{k — 1 )A+ 
16(A-w;*)« 


(8.3.2) 


Consider first the region 77/^* ^ argsj ^ In, \z\ ^ p(sinn/k)^^^» 
Then \z\^ < Asin77/i;, and 


< arg(A-z'') < ^ 

for all z in this region. Taking the integral defining ^( 2 :) along a 
straight line, so that njk 5 ^ argdw; ^ -|77, it follows that 

^ < arg{(A-M;*)ldw} < 

Z 

Hence I{^(25)} ^ sin^ J |(A— tp*)^ dw\ > KX^\z\, 

0 

where K denotes a positive number depending on k only. 
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Putting r]{z) — in (8.3.2), 

H(z) — •7(0)e^(*'coa^(z)+A~*ij'(0)e»'f<*'sin|(2)-f- 

Z 

J Qm{^{z)-—^{w)}R(w)H{w) dw, 

0 

If n is even, ipni^) is an even function, >^'(0) = 0, and this gives 


\m\ < \vm+ j\R{w)H(w)dw\. 
0 

Hence by Lemma 5.2 

\H(z)\ < h(0)|expjj|i?(M;)rfii;|j. 
It then follows from (8.3.2) that 

rj{z) = r}(0)coQ^{z)+e-^^^^'>x{z), 

where 


lx( 2 )l < h( 0 )||^ J |i 2 (w)|exp| J \R{w')dw’^\dw\ 

= hWll^expj J |i?(w)()fw|| — ij. 

Hence r)(z) — 

where 

\<o{z)\ < |e^'^<*>| + 2|x(3)|/|^(0)| 


< le^'' 


exp J \R{w) d^v\ — l . 


r . 2*-2 


= ol I dp] = 0(A-‘-i"'). 


A« 


a , J) 

Now J |i?(2^) dw| = o| J 
0 ' 0 
If we also assume that \z\ > Kp, it follows that 

co{z) = 0(A-^-i/^). 

The result is that, in the region considered. 


(8.3.3) 


(8.3.4) 

(8.3.6) 


8.4. The above result holds in particular on the segment of the 
straight line I(z) = ^ptamrlk between the imaginary axis and 
argz = -rrjk. We next require a similar result on the same line as 
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far as R( 2 :) = p, and then round the circle \z—p\ = ^ptoumrlk as far 
as the real axis. Let C denote this curve, together with its reflections 
in the real and imaginary axes. 

Since 


0 

0 


and 


j 




J e}ii^^E{w)7]{w) dw + 

0 

0 

J dw 


00 

J e~'^^^'^^R{w)ri{w) dw — 

0 

00 


J 




(8.3.2) is formally (i.e. apart from the question of the convergence 
of the integrals at infinity) equivalent to 


0 

ri{z) — J R{w)ri{w) dw + 

zi 

00 


+|i J e'^^^^^^'‘^^^^^R{w)rj{w) dw, 


(8.4.1) 


where A, B, and are independent of z. 

Let r be the region defined by 

0 < 1(2) ^ Ij^taiiTr/fc, \^—p\ ^ l^tanTr/A;, ^(z) > 0, 


together with the segment of straight line 1 ( 2 ;) = ^ R( 2 ;) p. 

Let Zi be the left-hand end-point of this segment, z any point of F. 
Take ^ = 0, i? — 1 in (8.4.1), and let the path from z^ to z consist 
of the curve G as far as the straight line through z parallel to the 
real axis, and then this straight line. Let the path from 2 ; to infinity 
be a straight line parallel to the real axis. 

If w is in r, \X~w^\ > K\w\^, where K depends on k only. Hence 


\R{w)\ < 


Kp^ 


since \w\ > Kp. Hence 


K K 


0 ^Kp ' 
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A similar argument applies to the rectilinear part of the path 2 ), 
and the circular part is of length 0(p), so that the integral of 
\R(w)\ round this is also Hence, if 8 denotes the maxi- 

mum of 

z 00 

I J J \R(w)dw\ 

Zl z 

formed in the above manner, we have 
Thus 8 < 1 if A is large enough. 

Now take A — 0, 5 = 1 in (8.4.1), and form a solution by itera- 
tion. Let 

y}x(z) — 

and 

z 

~ j* dw + 

2^1 

CO 

for V > 1. Then 

y 2 00 

rj 2 (z)—r]i(z) == J R{w) dw j dw 

' 2l Z 

W 

Now J (^—'^'^)^dw\ 

z 

where X—w'^ lies in the lower half-plane, and so (X—w'^)^ lies in the 
fourth quadrant (since it reduces to A* as w' ~>0); and dw' is real 
and positive. Hence 

RR^(2)— ^(««)}] < 0 , 

|g2i{f(2)-f(w;)>| ^ 

Hence \V 2 (^)~Vii^)\ ^ 

Similarly 

z 

Zl 

00 

J »ji{w)} dw, 

Z 

\Vs(^)-V2i^)\ < 

and so on generally. Hence ^ convergent, i.e. r]^{z) 
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tends to a limit 7]{z), which (as in the case of the similar formulae 
of § 6.2) satisfies (8.4.1) with ^ = 0, J5 = 1. Hence the corresponding 
ip(z) satisfies (8.2.2). If A is an eigenvalue A^, this i/j{z) must be a 
multiple of the eigenfunction since any other solution of the 

equation tends to infinity as 2 -> oo along the real axis. 

Now 00 

7}(z) = 2 + 

v-1 


say, where 

lxi(^)l < = IZTs == 0(X-i-V^). 

(8.4.2) 

The result is 

therefore 




(8.4.3) 


where is a constant. 


8.5, Let Zq and Z 2 denote the points of C on the positive real and 
imaginary axes respectively, and z^ the point defined above. 

The function ip^iz) is either odd or even — suppose, for example, 
that it is even. Let N be the number of zeros of i/^^( 2 :) inside C, 
Then 27riN is equal to the variation of log round C, and so to 
twice the variation of log ip^^z) round the upper half of (7. Hence 

N == il{log^„(-zo)-log^„(2o)}- 

77 

The values taken by between z^ and —Zq are the conjugates of 
those taken between Z 2 and Zq, Hence logi/f^(--2:o)— logi/f,^( 2 ; 2 ) is the 
conjugate oilog^^{zQ)—logi(j,,(z 2 ), and 

I{log<A„{-2o)-log<A«(2'2)} = I{log'/'«(22)-log’An(2o)}- 

Hence N = h{logil^Jz^)-log^Jz,)}. (8.5.1) 

77 


Hence by (8.3.5) and (8.4.3) 




+logi±^ + log 


(8.5.2) 


—ii(z^)+ii(Zj) = — J (X—v^)^dw 
c 

p p 

— j dw i j (X—w^)^ dw, 

-p 0 


Now 
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Next, as z varies from Zq to sjg, z^ describes a spiral, starting to 
the right of A and ending to the left of it, and not encircling it 
again. Hence arg( 2 :*^— A) increases by tt. Lastly, if A is so large that 
|w(z)| < 1, lxi(z)l < 1 on C, log{l+w(3)} and log{l+xi(2)} “ay be 
assigned their principal values; hence 

log{l+co(2)} = 0(A-i-y*), log{l+xi(z)} = 

Hence (8.5.2) gives 

p 

iV = - J (X-iv^)idw (8.5.3) 

0 

Now has n real zeros, and by (7.2.3) 

9 ^ 

n \ {X—w^)^dw —1. 

^ J 

0 

Also N = 72,+ 2m, where m is zero or a positive integer, since complex 
zeros occur in conjugate jjairs. Hence 

2m < 

If A is so large that the last term is less than |, then m = 0, N = n, 
so that there are no complex zeros inside the curve (7. Replacing 
N by n in (8.5.3), we obtain a formula equivalent to (8.1.8), with 
the remainder after the first term on the right-hand side replaced 
by an 0-term. 


8.6. Further approximations. The next approximation to the 
function r]{z) (0 < arga; < Trjk) is 


= ■^2(2;){1+X2(2)}. 

where 

Z 00 

y^^{z) = J R{w) dw J e--^ii<«’>R(w) dw, 

Zi Z 


and X 2 ( 2 ) = 0(82) = 0 (A-i- 2 /fc). 


Now 


00 

z 


f d I R(w) ) 

2i(X — 2*)^ J 2i ^ 


= 0{A-i-2/*|e-2<^(*)|}+o||e-2<f(*)| J 

==-- 6>{A-i-2/*|e-2<^(*)j}. 
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Hence rj{z) = J B{w) dw 

A similar formula may be obtained in the angle Trjk < arga: < \tt. 
It follows that the term in (8.5.1) may be replaced by 

— J R{w) dw 

c 

This corresponds to (8.1.8), with the remainder after the second term 
on the right-hand side replaced by an 0-term. We can plainly pro- 
ceed indefinitely in the same way. The formula (8.1.8) is therefore 
justified for the functions q{z) = 2 ;*, or for any analytic functions 
with sufficiently similar properties, provided that the eigenvalue is 
large enough. 


8.7. The case k = 4, Here 


B{w) = 


wHZp^+2w^) 


Consider first the formulae of § 8.4. The curve C now consists of the 
right-hand semicircle of \z—p\ — ^p, and the left-hand semicircle of 
\z~\-p\ — \p, joined by segments of the lines l{z) = On the 

upper half of the right-hand semicircle. 


\p—w\ = \p, \p+w\ > 2p, 2 p 2 


6-2 


and \w\ < ^p. Hence 

- 3^)^(9j,2/4) + 2(3^>/2)<> 27V2/ 27\ 1 

Integrating over a length ^np, we obtain approximately 4/p^. 
On the straight line w — u-[-^ip, \p ^ u ^ p, have 

5p^l4, 


\p+w\^ = {n-\-pY+\p^ > 


\ 2 / ' 2 




and 
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Hence 


< 5i.2i.34 = 4 

approximately. Integrating, we get again 4/p^; hence 


2 !l 

J 


On the straight line from z to infinity, where z is on or to the right 
of the semicircle, we have 

\p+w\ ^ \w\, \p^-\-w^\ > \w\^, \'^—p\ > 


Hence 


Hence this part is 


I RlwW < _ 3.2ij>i 2i 


p 

= (3.2l.i+2i.2)l = (1+16)^ = ^ 
17 

approximately. Hence ^ “3 • 

For 8 < I we therefore want 

A > 17* 44 

roughly. 

In §8.3, we have |^^;| ^ 2-lp, so that 

< w+v 


Hence 


Z 

J \B{w)dw\ < (AjiJ (Spy+2p^)dp, 


4957 


IT 
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which is found to be less than l/p^, and so it is small if A satisfies the 
above condition. A closer examination of the above integrals would 
no doubt show that the method was valid for much smaller values 
of A. 

8.8. General character of the functions If 

integral function of z, so is for each n. From (8.3.2) it follows 
that 

ri{z) = 0 [exp|I{|( 2 ;)}|], 

and hence that 

= 0[|?(z)|-iexp|I{^(z)}|]. 

For q(z) — z^, 

^(z) = ~i j dw 

0 

z 

= J Aw 

2o 

= — |iz2-|-|A„logz+0(l) (k = 2) 

= -‘p+i+o'‘» 

Hence 

>Pni^) == 0{|2|‘^”-iexp(i|2|2)} (k = 2), 

f«(2) = 0 j 1 2 1 ^ 

Hence iffj^(z) is an integral function of order Ji+l at most. Since 
by (8.4.3) 



for real x, the order is exactly \lc-\-l. Hence or \/zip,^X"^z) is 

an integral function of order \k-\-\, according to whether n is even 
or odd. 

In the case Ic ^ 2, (8.4.3) holds for all sufficiently large values 
of z in the first quadrant. The curve C can be replaced by any circle 
\z\ ~ R with R > A^. Consequently ^n{z) has no complex zeros, if 
n is sufficiently large. Actually this is well known to be true for all 
values of n. 

In the case Z; “ 4, iffjX^z) or ylzifj^^{^lz) is an integral function of 
order f, and so has an infinity of zeros, all but a finite number 
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of which must be comi:)lex. An infinity of these zeros are purely 
imaginary; for if 2 : — iy, (8.2.2) becomes 

g^+(2/^~A).A = 0. (8.8.1) 

The coefficient of 0 is positive for y > A^, indicating that the solution 
is oscillatory. For y ^ (2A)^, y^ — A ^ Hence (8.8.1) has a zero 
between any two consecutive zeros of a solution of 

«%■ of 

It is possible that in this case all the zeros of are either purely 
real or purely imaginary. For k ,> 4, the situation is probably still 
more complicated. 


REFERENCES 
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Kemble (2), Langer (4), (6). 

The method employed here is perhaps new. 



IX 

CONVERGENCE OF THE SERIES EXPANSION UNDER 
FOURIER CONDITIONS 


9.1. The problem of this chapter is similar to that of Chapter VI, 
but it is now assumed that q{x) -> oo. The expansion thus takes 
the series form 

/(^) = I ( 9 . 1 . 1 ) 

«=0 

It will now be shown that, by imposing further conditions on q{x), 
we can relax those which were imposed on f{x) in Chapter II; and 
in fact that, for a wide class of functions q{x), the expansion con- 
verges if f{x) satisfies conditions similar to those which are sufficient 
for the convergence of an ordinary Fourier series. 

As in Chapter VI, it is mainly a question of obtaining asymptotic 
formulae for the functions (f>(x, A) and ijj{x^ A) for large complex values 
of A; these then give what is required for the function 

€>(x,A) = J <l>(y,X)f{y) dy +<f>ix,X) J >fi(y,X)f{y) dy. 

® (9.1.2) 

We obtain satisfactory results provided that A is not too near to the 
positive real axis. A restriction of this kind is to be expected, since 
it is here that 0(a:,A) must have an infinity of poles. However, in 
the neighbourhood of the poles, the formula 


cI)(x,A)= 2 
w=0 


A A^ 


( 9 . 1 . 3 ) 


can be used. For this purpose we have to know something about the 
distribution of the eigenvalues A^, and the analysis of Chapter VII 
is needed. 


9.2. It will be assumed that q(x) is twice difierentiable. 


q\x) > 0, q"{x) > 0, and q\x) < {q'{x)]y 

for sufficiently large x, where 1 < y < |. Then Theorems 7.5 and 
7.6 hold. 

The analysis of §5.4 is required again here. With the notation 
there used and the conditions assumed above, we have 
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Lemma 9.2. J | -> 0 

0 

05 A -> 00 , it being supposed that A = u-\-iv, ^ 0, and A does not enter 

the region between the positive real axis and the curve v = 

x(u) being any continuous function which tends steadily to infinity with u. 

The conditions assumed imply that 

q'{x) = 0[{q(x)y], 

where c = 1/(2— y), by (7.3.4). Hence 


If u 


J \] \X-qm 

0 0 

^ 0, put q ~ u-\-vr. We obtain 


dt 


^ 0 


dq 


o 



{u-\-vrYv dr 


If « < 0, the integral is 



0 


q^dq 

(|Ap+?F, 


= 0(|AM). 


Since q^t) == 0[{q'{t)}v] = 0[{q{t)Y-^q' {% 


r 

J rA-?(oi* 



and similar results hold for this integral. The lemma clearly follows 
from these results. 


9.3. Let (l>(x. A) denote again the solution of 

g + {A-g(x)}2/ = 0 (9.3.1) 

such that ^(0,A) = sin a, <^'(0,A) = —cos a. For a fixed Xy or x in 
a fixed finite range, Lemma 2 of § 1.7 gives 

cl>{x,X) = cos(a:VA){sina+0(|A|“^)} (9.3.2) 

if sin od 5 ?^ 0, and 

^(a:,A)= -?iH^{cos«+0(|A|-»)} (9.3.3) 


if sin a == 0. 
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For x->oOy the relevant formulae are (5.6.1), (5.6.2). Also by (5.6.1) 
and Lemma 9.2, the integral on the right of (5.3.6) is 


0 


ou 

\X\ij |i?(0| 


dt 


o(|A|i) 


as A -> 00 in the region allowed by the lemma. Hence 


and 


ilf(A) ~ lA^sina 


(9.3.4) 


(9.3.5) 


as A 00 in the above region, provided that sin oc ^ 0. Similarly, if 
sin a — 0, 

, . cos oc 


M{X)^ 


Ai ‘ 


(9.3.6) 


The formula obtained by diiferentiating (5.6.2), 
d 


dx 


[{X~-q(x)}^(l>(Xy A)] iMiX){X--q{x)}h-^^^^\ (9.3.7) 


is also easily proved by first differentiating (5.4.5). 


9.4. It is now a question of obtaining a solution of (9.3.1) which 
is small for large x. We proceed as in § 6.2 or § 8.3, but is now 
small for large x, not large as in §8.3. Hence (5.4.5) is formally 
equivalent to 

7]{x) =:= |i7y'(0)A-^{e"f(^>— 


+ 


or 


^ J e^i((x)-m)lt{t)7}(t) dt +~j dt - 

0 X 

00 

2i J 
0 

X 

7){x) = f dt + 

2 % J 
0 

00 

+ ( eim-((x))m)-q{t)dt. (9.4.1) 

21 J 
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Take ^ = 1, i? = 0, rji = and for n > 1 

X 

0 


+ 1 




Let J now denote J |i3?(^)| dt\ by the above lemma, J < 2 if A is 
0 

sufficiently large and in the above region. Arguing as in §6.2, it 
follows that r]^^^{x) tends to a limit r]{x), which satisfies (9.4.1) with 
A ~ 0, B \ and 

\ri{x)\ < (9.4.2) 

|^(a;)_e»-a^)| (9.4.3) 

In particular, y{x) = is L^^q^oo). 


9.5. We have to apply this result to the function ij}{x,X). Now 
\jj{x,X) ~ A{X)(j>{x,X)A-B{X)y{x,X) 

and iIj(x,X) and y{x,X) are /v^(0, oo), while (f>{x,X) is not, at any rate 
if M(X) ^ 0. It follows that, if A is large enough and in the region 
of the lemma, 

iIj(x,X) == B{X)y(x,X), 

Since W((f),^) 1, B{X)W(^,y) — 1. 

Now for a fixed X, x oo, (9.4.1) (with J. = 1, ^ = 0) gives 

riix) ~ 1 +1 J e~»^Wi?(<)7?(<) dt \ = if(A)e»f(*> 

say. As in previous cases, we also have 

7]' {x) ~ iK(X){X — q{x)ye'^^^^K 

Hence 

W{<l>,y) = {X-q(x)}-^W[{X-q(x)}ii,(x,X),r](x,X)] 

r^{X-q(x)}-i[M{X)e-ii<^HK(X){X-q(x)}h^i^-^^+ 

+iM(A){A-g(a;)}ie-^^(*)is:(A)e<^(*)] 


= 2iM(X)K(X). 
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Hence 

i.e. 


2iB(X)M{X)K{X) = 1, 

2/(a;,A) 


A) 


2iM(X)K(X) ' 


(9.6.1) 


Clearly K(X) -> I as A -> oo in the above region. 

For a fixed Xy and A -> oo in the above region, this gives (if 
sin a 0) 

Since ^ 

^(x) = J {X-q(t)}i dt 
0 


X 

= J {Ai+0{A-i)}d« 
0 


= a;A*+0(A“*), 

this also gives 4f{^iX) ~ r — • 

^A^ sin OL 


(9.5.2) 


fiix'^lX 

Similarly j/f'(a;, A) — . 

sin a 


(9.5.3) 


If sin a = 0, the corresponding results are 


0(^, A) 


pix^X 


cos a 


f(x,A) 




COS a: 


(9.5.4) 


9 . 6 . All the above results apply to regions which do not come too 
close to the positive real axis in the A-plane. In the neighbourhood 
of the positive real axis the expression (9.1.3) has to be used, and 
for this purpose some further information about the functions ipj^{x) 
is needed. 

We have ^Ux)+{X^—q{x)}^ni^) = 0- 

Let qiPn) = K- Then for 0 < a: < the coefiicient of ^^(a;) is 
positive. Hence «/^^(a;) is concave downwards where it is positive, 
and upwards where it is negative; it therefore has just one maximum 
or minimum between consecutive zeros. 

Lemma 9.6(a). The successive maxima of \^jSx)\ in 0 < x < 
form a non-decreasing sequence. 

m 

A„-2(a;)‘ 


Let 





9.6. 9.7 


UNDER FOURIER CONDITIONS 


163 


Then F'(x) = > 0. 

Hence F{x) is non-decreasing. At a critical value of = 0, 

so that F(x) ~ The lemma therefore follows. 

Lemma 9.6(6). Let he the zeros of xfs^fx). Then 

WnMl 

{ ^iVv) 

is a nan-decreasing function of v. 

Since F{x) = 

at o; “ 7^1, 7^2,..., this also follows from the above argument. 

9.7. Lemma 9.7. For any fixed x {or x in a fixed interval) 

^;(x) == oMp-i). 

Let the successive maxima (or minima) and zeros of in the 
interval 0 < a; ^ be 

^l<Vl<h< - < Vn- 
Consider a particular Then by Lemma 9.6 (6) 

l <A»(^ft) l < \ ]K(Vv )\ 

Since q(x) is convex downwards, q(^x) ^ lq{^)- Hence 
^iiPn) < kiPn) = 

Hence \<l'n{v^)\ < 2 

^ V^-fl 

Suppose e.g. that is positive between and rj^,. Then 

/*•<")'*“= I 

== <An{’?J) 


4957 


X 
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Hence 


Hence (9.7.1) gives 




Wn{'n^)\ < J dx 


- iPn iPn 

j J {'Pn(x)Y dx J dx 


< A, 
n' — /x+^ 


^ (9.7.2) 

n' — /x+1 

Now let a: be a given positive number. Let be the nearest 

zeros of i/fniv) below and above x, the zero of ip^iy) between 
and Then \ipn(x)\ < \^n{v^)\- if is clear from the argu- 
ment of § 7.2 that ip 


i'-M {K~Qiy)}^ dj 

X 

iPn 


dy -1 


for n large enough. Hence 
and the lemma follows. 


(9.7.3) 


(9.7.4) 


9 . 8 . Lemma 9.8. For a fixed x 

Ux) = 0{p-^). 

Consider again an interval (^^, ri^) where > 0, «/f^(a:) < 0, 

x ^ Then 

'l>'n{xWn{x) = {K-<l{x)}xl)Sx)[-^'n{x)} 


> ¥a4>n{x){—'P'n{x)). 


J 'PUxWnix) dx ^ j <f’nix){—<Ph(x)} dx, 

i i. 


Hence 
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Hence by (9.7.4), if < x < 
and the lemma follows. 

It is easily verified that, in particular cases, these results are 
the best possible. For example, in the Hermite case, q(x) = x^y 
= 2n-\-l,p^ = (2?i+l)^, so that Lemma 9.8 gives — 0{n~^). 


Actually 


>Pni^) = 




(n\ 2'^7T^)^ 

Now (Szeg5, p. 94) for a fixed x 

= 0{p-j 


n\ 




Hence 


•I'M - 


»r(|n+l)i 


= O(w-i). 


9.9. We can now prove 

Theorem 9.9. Let f{y) belong to L^(0, oo), and let f{y) ~ 0 in an 
interval of which x is an interior point. Then 

n -0 


This shows that, as in the case of ordinary Fourier series, the 
convergence of the series ^ depends only on the behaviour 

of the function f(y) in the immediate neighbourhood of the point x. 
Of course we have also assumed that f(y) is U, a condition not required 
in the theory of Fourier series. 

Lct/(^) — 0 for x—h < 2/ < a;-fS. 

We have to prove that, under the conditions of the theorem, 


lim J <!){x, X) dX = 0 
imIr 

as i2 -> 00 through values not equal to any of the eigenvalues. We 
have 

r a>(x, A) = ^ Cn'PrM) I == 2 






n=0 


(9.9.1) 
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say. Now 2 is convergent, since /(x) is L^(0,(X)). Let 


Then 

J A)cJA 




|c™</'»(*)4|-h I K*/>n(^)I„ 

Xn>hR 




Now is bounded for all values of R and n (e.g. since we can take 
it along an arc of a circle with centre at A = A^). Since 

|A— A^J > R—X^, 

and the length of the contour is 2 R 7 ], it is also 

o( V 

From this and Lemma 9.8, it follows that 

say, and 

Y \>Pn(^)In? - of y 


+ 


^R^Xn^R—Efl 

+o( 2 -l+®[ 2 ,T 

0{S,)+0(S,)+0(S,), 


R^T)^ 


-R?P, 


say. 

In S^, p{R — Rr]), and the number of terms in the sum is 

N{R-\-R'q)-N{R-Ry^). Hence 

by (7.6.2), provided that 

^{R — Ry]) ^ 2Rrj ^ R — Rrj, 


which is true if 


1 
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00 

".-2 2 




(K-R)^P„ 




fe=l 

2 2 isV. 

1 N{R+2'‘Rr,)-N(R+2'‘-^Rr)) 


2i 


fc-1 


22k-2 


P(R-{-2^-'^R7]) 


k=l ' ' 


Similarly S 2 ~ 0(R^r]). 

Lastly a similar argument gives 

2-'‘-^R 




0 


\2ik)’ 


2 

and summing terms of this type, 

Si = 0{Riri^), 

Altogether f ^{x,X) dX = 0(R^rj)+0{€(^R)R^ri^, 

This tends to zero if is a suitable function of R, e.g. if 

_ 1 

“ is:*{6a:B)+iF»} 

(or of course for smaller rj, but later we want R^rj oo). 


(9.9.2) 


(9.9.3) 


9.10. Next consider the order of 0(a:, A) on the parabola ~ u. 
We have 


n=0 ^ 

-«(2<x4^.)*- 

3 not exceed 

2 U^Pu'^ 2 {u—XnfPn 


(9.10.1) 


The last sum does not exceed 

1 

^Pn 






T — 

jL^ 


U-V<Xn^U-\-V 

+ 2 . {Xn-ufPv 


Pn 


w+v<A« 

We consider these sums in the same way as those of the previous 
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section, with B and Br) replaced by u and v respectively. The third 
sum, for example, is 

p(u—v) I {v^u^j \u/’ 
and similarly for the others. Hence on the parabola 

<I>(x,X) = 0(u-i) = 0(|A|-i). (9.10.2) 

This also holds throughout the part of the first quadrant above the 
parabola, since (9.10.1) is a decreasing function of v, for a fixed u. 
In the second quadrant 



0(a:,A) - 0 






the convergence of the series following from the above argument. 
Now we also have 

x —8 00 

0(x,A) = ip(x,X) j <^(i/,A)f(i/)dp +<f(x,X) j tjj(y,X)f{y)dy. 

0 a;+S 

Combining (9.3.2), (9.3.3), (9.3.5), (9.3.6), 

<l>{x,X) = 0{|A~W(A)e-^^>'^|} 

for a fixed x, or a; in a finite range. By (9.5.1) and (9.4.2) 
i/r(a;,A) = 0{|A-Mf<^V^(A)|} 


for all X, and A in the region allowed by § 9.2. Let A be real and 
negative, say A = —[i. Then 


Hence 


X 

I{f(a;)} = J {/a+g-COP dt > x^y. 
0 

( X S \ 

^^■^>^\f{y)\dy\-\- 


+0 




00 

■ f e-y^>^\f{y)\ dy 

xU 


The first term is clearly and so is the second, since the 

integral does not exceed 



<I)(ar, -y) == (9.10.3) 

F{X) = Aic-^%(x,A), 


Thus 

Let 
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where is real on the positive real axis. By (9.10.3), F{\) is bounded 
on the negative real axis. On the parabola = u 

^ exp| — ^ + ...| == 0(1). 

Hence, by (9.10.2), F{y) is bounded on the parabola. It is therefore 
bounded throughout the region between the parabola and the nega- 
tive real axis, by the Phragmen-LindelOf theorem. Hence, putting 

'X id 

’ 0(a;,A) = 

Hence, if 0 < jBsin^i^ ^ cosrj, 


[ <D(a:, A) dA = 0 ( f dd) == 


\X\=^R,7]<0<TT 

lame result h 

Theorem 9.9 now follows from (9.9.2) and (9.10.4). 


(9.10.4) 

The same result holds, of course, for the integral over — tt < 0 < —rj. 


9.11. In view of the above theorem, we can now confine our 
attention to a finite interval. We first prove 

Lemma 9.11. Letf{x) be of bounded variation over a finite interval 
(a, 6), and zero elsewhere. Then 

Let f{x) ~ f fix)— f^{x), where and /2 are bounded and non- 
increasing in (a, 6). Then 

b h 

= J 'l>n{y)U(y) ^y - j 'l>n(y)k{y) dy- 

a a 

By the second mean-value theorem 

6 j8 

J 'J>n(y)fi{y) dy = /i(a-f 0) J dy (a < ^ < b) 

a a 


/i(<^+0) 


J 4'niy} dy {a <(X<^) 


= 0{X-ip-i) 
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by Lemma 9.7. Similarly for the other integral. Hence by Lemma 9.8 

cM^) = 

By Theorem (7.3) 

n — 

K^p-^ - 0(n-^) 

and the lemma follows. 


9.12. Theorem 9.12. Let q{x) satisfy the conditions stated in § 9.2, 
and let f{y) be of bounded variation over an interval including the point 
and zero outside this interval. Then 

= - 2 -{/(^+ 0 )+/(a;- 0 )}. 

n ~0 

We start again from (9.9.1), but now use Cj^ijj^fx) — 0{lln) instead 
of 1^2 < 00 . We obtain 

'jV,A)dA = o( y 2 )+o( y ^„%1+ 

L ^ .n<X^n~Rr, 






say. Now 


0(Ix)+0(I.)4-0(23)4-0(2,), 
® 1 N{B-Rri) I - 


by (7.6.3), provided that rj > Next 

I. =2 2 -- 


k=l B+%’‘-^BT]<X„^R+^Rr] 




OO 

2 


N{R+2'‘Rrj)-N(R+2’‘-^R-q) 1 


N(R+2k-^R7}) 


2k-l 


21 

k=l ' 


/ 1 




R+2'^~^R7j 2*- 

2 2 

(,logi). 


: 0 


+ 0(r]) == 0( 
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Similarly for ^ ; and 




iR<Xn<iR 


by (7.6.3), and adding O(logi?) terms of this form, 

2i = O(ijIogJK). 

We may take for example t] = i2”“, where 0 < a < ^. Then 


/ ‘I>KA) 


dX = 0. 


Now A = i?exp(ijR““) ~ lies in the region specified in 

Lemma 9.2, if a is in the above interval. Hence on the remaining arc 
of the circle |A| = jR we may apply the asymptotic formulae obtained 
in §§ 9.3-9.5. 

Now 

0(x,A) = il>{x,X) J <f>(y,X)f(y) dy +^(a:,A) J if‘{y>^)f{y) dy- 

x—S X 

The second term can be written in the form 
x+8 aj+5 

f(x+0)<f>(x,X) j il,{y,X)dy +^x,X) j i>{y,X){f(y)-f{x+0)} dy 

X X 

= <I)i+02. 


x+d x+ 

J 4i(y,X)dy = - J 


f'(y.A) , 

X-q{y) 


f(x,A) ifi'{x+8,X) 
X—q{x) X—q{x-{-S) 


X-t-0 

I -A'(2/.A) 


g'( y) dy 


Suppose, for example, that sinoc ^ 0. Using (9.3.2) and (9.5.3), 
= /(x+0)cos(a;VA){sin «+ 0{ lA|"i)} X 

^ |{A— ?(x)}sina \ [Al \ lAl^ // 

. . , n^cos(xVA)e‘■*^^ _ 1 , 


■'/(*+0)- 


-/(»+0)g 


H- 


Hence 


lim f Oi(x,A) dX = if{x-\-0). 

R—^oo . . . 


largAj^T; 
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In ^> 2 , we can write 

' fiy)—f{x+o) = gi(y)—g 2 iy), 

where and are non-decreasing and tend to 0 as y -> x. Now 
ar+S aj+S 

f R{'/'(y>^)}gi(y) # = o) J Ri/>(y>^) dy 

X X' 

by the second mean- value theorem; and as before this integral is 
0{|e^^^^/A|}. Similarly for the imaginary part, and for the corre- 
sponding integrals involving ^ 2 * Altogether we obtain 

* — /o/ gi(^+^)+g2(a^+8) 

^ \ |A| 

Hence J C>(a;, A) = 0{sri(a;+S)+fii'2(a;+8)}. 

|argrAl>i7 

Since 8 can be taken arbitrarily small, this contributes 0 to the final 
result. 

Similarly for the other part of A), and the theorem follows. 

9 . 13 . Combining Theorem 9.9 and Theorem 9.12, the final result 
of the analysis of this chapter is 

Theorem 9.13. Letfiy) belong to oo), and letf(y) be of bounded 
variation over an interval including the point x. Let q{x) satisfy the 
conditions stated m § 9.2. Then 

n = 0 

REFERENCE 

Titchmarsh (10). 



X 

SUMMABILITY OF THE SERIES EXPANSION 


10 . 1 . It will be assumed here merely that q{x) is continuous and 
tends steadily to infinity, but not that it satisfies the other conditions 
of Chapter IX. The asymptotic formulae previously used are no 
longer available, or at any rate they have not been proved under 
such general conditions. Consequently it is no longer possible to 
prove the convergence of the eigenfunction expansion under Fourier 
conditions. We can, however, prove the following summability 
theorem. 

Theorem 10.1. Let q{x) be continuous and tend steadily to infinity, 
and letf{x) be L^{0,co), Then 

CO 

=fi^) ( 10 . 1 . 1 ) 

for every value of x for which 

xiv) = J dy == o{ri) ( 10 . 1 . 2 ) 

0 

as 7] In ^particular (i) it holds almost everywhere, and (ii) it holds 
wherever f{x) is continuous. 

In view of (9.1.3), this is a question of the behaviour of 0(a;,A) 
on the negative real axis. Putting A = —t^, it becomes a question 
of the solutions of 

^ = {t^+y{^)}y 

for real t, and q{x) tending steadily to infinity. This is dealt with 
in the following lemmas. 

10 . 2 . Lemma 10.2. Let Q{x) > 0/or x > x^, and let y{x) be a solu- 
tion of 

g = «(»)!' 

belonging to L'^{xq,co). Then y(x) and y\x) have opposite signs, and 
tend steadily to zero as x-^ oo. 

Suppose, for example, that y(xQ) > 0. Then y'{xQ) < 0, or y{x) 
would tend steadily to infinity, since it is convex downwards where 
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it is positive. Thus y'{x) < 0 as long as y{x) > 0. But y{x) cannot 
change sign; for if x^ were its first zero, and y{x) < 0 for a; just 
greater than x^, there would be a point such that 2 /(^i+S) < 0, 
Hence y{x) would decrease steadily past this point, 
which is impossible if it belongs to L\Xq, oo). 

It follows that, for x > Xq, y{x) > 0 and y'{x) < 0. Hence y{x) 
decreases steadily, and so tends to a limit, which must be zero if 
y{x) is L^. Since y\x) > 0, y\x) increases steadily, and so tends to 
a limit, which clearly must also be zero. 

10.3. Lemma 10.3. In addition to the above conditions, let Q{x) be 
non-decreasing. Then, if yix^) > 0, 

y{x) < 2 /(a;o)expr— J du\ (x > x^). (10.3.1) 

We have for x ^ Xq 

—y'Wix) = Qix)y{x){-y'{x}} > Q{Xo)y(x){—y'{x)}. 
Integrating from Xq to x^, 

W^i^o)—y'^(xj}} ^ ^Q{xo){y^(xo)—y%x^)}. 

Making x^ oo, and replacing x^ by x, it follows that 

y'Hx) > Q{x)i/(x). 

Hence —y'{x)/y{x) > {Q(x)}K (10.3.2) 

Integrating from Xq to x, 

X 

^ogy{x^)-\ogy{x) ^ J {Q(u)f du, 

Xq 

and the lemma follows. 


10.4. Lemma 10.4. Let Q{x) satisfy all the above conditions, take 
Xq — 0, and let 

y(0)cosa+y'(0)sinQ' = 1, (10.4.1) 

where sin a: ^ 0 {instead ofy{0) being given). Then, if Q{0) > cosec^a. 


\y{^)\ < 


r a? 

exp — J {Q(M)}t du 

L n 


{Q(0)}i|sino£| — 1 

By (10.4.1) and (10.3.2) 

|y'(0)sina| < l+|y(0)| < l+{Q(0))-»|y'(0)|, 

1 


(10.4.2) 


l2/'(0)| < 


sma 


-{<2(0)} 


-i' 
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Also by (10.3.1) and (10.3.2) 
\y{x)\ < l2/(0)|exp 


WlexpT- J{Q(tt)}idJ 

L 0 

'(0)|{e(0)}-^exp[- j{Q{u)}idu 


and (10.4.2) follows. 


10.5. Lemma 10.5. Under the same conditions 

U,M\ 2y(a:o)exp[-(x-a;o){Q(a:+l)}i] 


(10.5.1) 


Suppose, for example, that y{x) is positive. Since Q{x) is non- 
decreasing, 

—y’{x)y”{^) = Q{^)y{x){—y'(x)} < Q{^i)y{^){—y'{x)] < *i)- 

Integrating from x to x^ 

\{y'\x)-y’^{x^)) < \Q(x^){y^{x)-y^{x^)} 

< \Q{xx)y\x)- 

Writing temporarily a = b == — 2/'(^i)» l^his gives 


Hence 




X — Xq, 


where yQ = y(xQ), Hence 

y+ /fy^+y >(2/0+ /fe+yie-“<"^-*'’> > 2y(,e-' 


Also, since y is positive and convex downwards, 

ii,-(»)i < y& =yj ^ . (0 <*■<»). 

X X 

Taking x’ — x— I, this gives 

\y'{x)\ < 2/(a:— 1). 


(10.5.2) 
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Hence, taking = x+l in (10.5.2), 

y(a^)+ ^ 2yoexp[-{e(a:+l)}n*-a:o)]. 

and (10.5.1) follows. 


10.6. Lemma 10.6. Let Q{x) = i^-\-q{x), where q{x) is non- 
decreasing and independent oft, and let y{x) satisfy the above conditions. 
Then as t-^oo 


y{x) ~ y{(i)e-^ 


( 10 . 6 , 1 ) 


uniformly in any fixed range of values of x, if y{0) is given) while if 
(10.4.1) holds, 


y{x)r^ 


g-X/ 

t sin a ’ 


( 10 . 6 . 2 ) 


The first result follows at once from (10.3.1) and (10,5.1). 

It follows from (10.3.2), with a; == 0, that l2/'(0)/2/(0)| tends to 
infinity with t, and hence, if (10.4.1) holds, that y\0) has the same 
sign as sin a, if t is large enough. Hence y{x) has the opposite sign 
to sin a. Thus the sign in (10.6.2) is correct. It will then be sufficient 
to consider one case, e.g. that in which y{x) is positive. 

The required upper bound for y{x) follows at once from (10.4.2). 
Also (10.4.1) gives 

l2/(0)l + l«/'(0)sina| < {Q(0)}-i|t/'(0)|+|2/'(0)sina| 


by (10.3.2) with a; = 0. Hence 
\y'm > 

Also, if Xi > X, 

y{xi)-y(x) = I y'(0 d$ 


1 


{Q(0)}-*+sina 


Xi 


(10.6.3) 


X 

= y'(x)(x,-x)+ d$, 


y'(a;) = L_ r Qi^)y(^){x,-0dt 

X^ — X Xi — X J 

X 

and hence W{^)\ ^ + 


x^—x 
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Now, by (10,5.1), 

y{^i) ^ 0 ( 1 / 0 } 

if X and x-^ lie in a fixed interval. Hence 

T akin g (xi—x)t small and large, e.g. x^ = x-j-f-^, we obtain 

\y’{x)\ < ty{x){l+h) [t > <o(S)). (10.6.4) 

The lower bound for y{x) implied by (10.6.2) now follows from 
(10.5.1), (10.6.4), and (10.6.3), 


10.7. Proofof Theorem 10.1. Since/(a;) is L2(0,oo), (9.1.3) holds, 
and the theorem is true if 

<^{x,-t‘^) f{x)lfi (10.7.1) 

as i ->■ 00 . Now 

X CO 

<l)(x, — <2) = ifi{x, -<2) J cf,{y, —t^)f{y) dy ■\-<t>{x, —t^) J 4>(y, —t^)f{y) dy. 

0 X 

Since ^(0) = sin a, ^'(0) = —cos a, 6(0) = cos a, 0'(O) = sino^, and 
i/j(x,X) = 6{x,X)+m(X)<f>(x,X), 

we obtain i/f(0, — ^^)cosa+^'(0, — ^^jsina = 1. (10.7.2) 

Suppose first that sincx ^ 0. By (1.7.3) 

^(o:) = cosh sin a+O(e^70 (10.7.3) 

as t->oo, uniformly in any finite a;-interval; and y(x) == ^(x,—t^) 
satisfies (10.6.2), uniformly in any finite a;-interval. Also by (10.4.2) 

ifi(Xi —fi) = 0{e~^^) 

for all X and sufficiently large t. Hence, if S > 0, 


^(x, -fi) J <f>(y, -my) dy = 0 ^ J e*"|/(y)| dy\ = 0(e-^) 

0 'o' 

d 

1 Ipiy, -my) 1 e-y‘\f{y)\dy\ 

xU ^ ds ’ 

= J e-^y*dy | \f(y)?d^ = 0{e-^). 



168 SUMMABILITY OF THE SERIES EXPANSION Chap. X 

The remainder of <!)(*, — t^) can be written 


>I>(X, — <2) J <l,{y, dy +<l>(x, —t^) J ijiiy, —t^) dy\ + 

X—h X 

■\-^{x, — «2) r ^y, —fi){f{y)—f{x)} dy + 

+ ^(X, -<2) J 0(y, —ti)[f{y)-f{x)} dy. 


X—Z 


The coefficient of f{x) is asymptotic to 

a+8 


^-xt r r ^-yl Y 

: — ley^^iYiocdy — : — dy = — 

^sina J ^ ^ ^ J ^sina ^ ^ ^ 

x-S X 

If (10.1.2) holds, the last term in the previous expression is 

x+8 


(XTC 

grf J 


t 


■ \f(y)-f{^)\ dy 


= 0 I J e-^\f{x+u)—f{x)\ du 
' 0 

( ® 

= oj|[e-«<x(«)] 5 + J du 

^ 0 
= 0{e-^/t)+0^ j e-«‘udu'j = o{t-% 


Similarly for the second term; altogether 

a)(x,— <2) = — /(x)/<24-0(i-2), 

the required result. A similar argument holds in the case sin a = 0, 
using (1.7.5) and (10.6.1). 


10.8. An alternative proof of Theorem 9.12. There is a well- 
known argument in the ordinary theory of Fourier series by which 
we can deduce the convergence of the Fourier series of a function 
of bounded variation from its summability (0, 1). If f{x) is of 
bounded variation over (0, 2tt), and 6^^ are its Fourier coefficients, 
then 
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The series 2 

n=l 

is summable ((7,1) to the sum f{x) for each x; hence, by Hardy’s 
‘Tauberian’ theorem, it also converges to f(x). 

We are now in possession of the materials for a corresponding 
proof of Theorem 9.12, viz. Lemma 9.11 and Theorem 10.1; but, 
instead of Hardy’s Tauberian theorem, we have to use a case of 
Wiener’s general Tauberian theorem. 

Let = c^ils^(x), so that a,^ ~ 0{lln) by Lemma 9.11. Let 

S = K/(a^+o)+/(a;— 0)}. 


Then by Theorem 10.1 


“2 Art 

n=0 ^ 


as A -> 00 . Let 
Then 


An^X 


A+A, 


i-i 




Now 


>S(A)+A<D(x, -A) = y ajl 

T y l==l y 


-^1 

w(A+A„)/ 


A<A„<2A 


A<An<2A 


iV(A)<n<V(2A) 


_ ^jl i\^(2A)-iy(A) ) _ ^11 


A N{X) i 
by (7.6.3). Replacing A by 2A, 4A,... and adding, we obtain 


Similarly 

Hence 


4957 



2 ^ 

V ^ 

A„<A 


0(A).. 


;Sf(A)+A<D(x,-A) = 0(1), 

^(A) = 0(1). 

z 
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Also 

Hence 


where 

Now 

and 

Hence, 

where 

Thus 


00 00 

- J "W = J 


;Sf = Urn 


00 

J 


A->ooJ {X+uf 


S{u) du 


00 

= lim f 

1-00 J 


e( 


(e^+evf 
e^~v 


S(e‘^)e‘n dr) 


—00 

00 

= lim f Xi(i—7))S(eV) d-q, 

A-^ao J 


K,{x) 


(e*+l)2‘ 

00 00 00 

(* C px-iux C f-iu 

—00 — 00 0 

= r(i+m)r(i-w) 0, 

00 

J Kj^{x) dx = 1. 

— 00 

by Wiener’s Theorem 4,f 

00 00 

lim f K^{^-r))8{€n) dr) ^ 8 \ K^{x) dx, 

A—yco J *' 

^ —00 — oo 

m)=^0 (^<0), e-^ (^>0). 

i 

8 = lim r eV-i8{e‘^) dr) 


X 

= lim - r 8{^J) dy. 

a;— ►CO X J 


t N. Wiener, The Fourier Integral, pp. 73-4. 
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Let S be a number between 0 and 1, independent of x. Then 

,(14-S)a; 

(1+8)S-/Sf l' ‘ 

8 = = hm — 

0 iC->00 OX 


-I J s(y) % — J ^{y) dy 

' 0 0 
(l+8)a; 

f S(y)dy 

X-*oo OX J 
X 

p (1+8)® 


Now 


s(s,)-«w = 2 '•• = 2 

X<Xn<V X<Xn<V ' ' 

= °{ I s) = o{ 2 s) 


^X<An<(l + 8 )X 

'[ 


^iV(a:)<n^M(l+8)x} 

■1 


= = 0(8) 


by (7.6.3), provided that x-^ ^ 8. Hence 

a+8)a; (1+8)0; 

^ I {'Sf( 2 /)-'Sf(x)} dy = i J 0(8) dy 

X X 

= 0 ( 8 ). 

Since 8 is arbitrary, it follows that 

lira 8{x) = S. 

x - y < x > 


Titchmarsh (9), (10). 
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